DOE/NASA/0615-79/1 
NASA CR-159465 


DEVELOPMENT OF A PHASE CHANGE 
-THERMAL STORAGE SYSTEM USING 
MODIFIED ANHYDROUS SODIUM 
HYDROXIDE FOR SOLAR ELECTRIC POWER 
GENERATION 


Barry M. Cchen 
Richard E. Rice 
Peter E. Rowny 
Comstock & Wescott, Inc. 



December 1978 


Prepared for 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
Lewis Research Center 
Under Contract NAS 3-20615 

(R1S4-CR-159465) DEVEIOPHEHT OP A 
PHASE-CHAHGE '^HERHAl STOBA6E STSTEH OSING 
HODIFIED AHHTDBOOS SODIDN BTDBOZIDE PCS 
SOLAB ELECTRIC CO0EB GENEBATICN <COBStOCk 
for and Bescott, Inc.) 252 p HC A12/HF A01 

U.^DEPARTMENT OF ENERGY 
Office of Solar, Geothermal, Electric and Storage Systems 
Division of Energy Storage Systems 


B79-19454 

Onclas 

G3/44 16441 



TABLE OF CONTENTS 


SUMMARY • 1 

INTRODUCTION 3 

Objective 3 

Concept for Phase-Change Thermal Storage 3 

Heat Storing Medium ^ 

Background 5 

Scope of Work 6 

TECHNICAL DISCUSSION 9 

MATERIAL PROPERTY STUDIES 9 

REFERENCE DESIGN STUDY 11 

Computer Model 11 

General Description 12 

Reference Design Analysis 15 

SUBSCALE TEST UNIT 37 

Scaling Procedure . 37 

Subscale Unit Design 

Construction of Subscale TES Unit ^6 

SUBSCALE TEST SYSTEM 6l 

SUBSCALE TEST RESULTS AND EVALUATION 71 

Preliminary Discussion 71 

Description of Test Procedures 74 

Comparison of Computed and Experimental 

System Response 76 

One-Pass Charge/Discharge Test Results .... 76 

Idle Tests 92 

Square Wave Cyclic Testing 93 

Solar Cycle Simulation Tests 104 

PRELIMINARY DESIGN 121 

Cost Estimates 135 

CONCLUDING REMARKS 137 

HRECEDINQ PAGE BLANK NOT FOKO 

ill 



APPENDIX A 

DYNATECH R/D COMPANY REPORT 

"The Thermophysical Properties 
of a Thermal Energy Storage Material” 


APPENDIX B 

COMPUTER MODEL USED FOR REFERENCE DESIGN STUDY 

AND PRELIMINARY DESIGN D-1 

Main Program B-1 

Subroutine ELANAL B-5 

Initialization Section B-6 

Heat Transfer Section B-10 

Secondary Subroutines and Functions B-19 

APPENDIX C 

ORIGINAL COMPUTER ANALYSIS C-1 

Introduction C-1 

Description C-2 

Main Program C-8 

Conclusion C-15 

FORTRAN CODE C-1 7 

APPENDIX D 
LIST OF SYMBOLS FOR COMPUTER ANALYSIS 


APPENDIX E 
LIST OF SYMBOLS FOR SCALING PROCEDURE 


APPENDIX 

REFERENCES 

iv 



LIST OF TABLES 


TABLE I 

Cost Calculation Input Data 

Table II 

Cost Calculation Output Data 



LIST OF FIGURES 


Figure 1 • 13 

Typical Solar daily cycle at Sandla STETP 

Figure 2 19 

Overall efficiency as a function of Thermkeep mass, 
showing effect of heat exchanger surface area. 

Figure 3 20 

Overall efficiency for constant Thermkeep mass 
at various heat exchanger surface areas. 

Figure 4 21 

Cost of storage for constant heat exchanger 
surface area at various Thermkeep masses . 

Figure 5 22 

Cost of storage for constant Thermkeep mass at 
various heat exchanger surface areas . 

Figure 6 23 

Specific cost of storage for various TES designs 
at fixed overall efficiencies. 

Figure 7 2^ 

Reference design with 10 element grid (preliminary). 

Figure 8 25 

Reference design with 25 element grid (preliminary). 

Figure 9 26 

Reference duty cycle- 

Figure 10 29 

Effect of tube size on cost and performance 
(25 element grid). 

Figure 11 30 

Effect of vessel aspect ratio on cost 
and performance. 

Figure 12 31 

Effect of fluid velocity on cost and performance 
(10 element grid). 

Figure 13 32 

Effect of Insulation thickness on erst and 
performance (10 element grid). 

Figure lil 33 

Effect of peak allowable storage flow on cost 
and performance. 

Figure 15 35 

Final full-scale reference design — movement of 
thermal gradient . 

vi 



Figure l6 36 

Final full-scale reference design — variation of 
Therralnol-66 outlet temperature during discharge. 

Figure 17 ^7 

1/10 model design Thermkeep thermal energy stor- 
age unit. 

Figure l8 ^8 

1/10 model design — Thermkeep thermal energy 
storage unit. Variation of Therminol-66 outlet 
temperature during discharge. 

Figure 19 ^9 

Photograph - Subscale vessel and partially 
Inserted thermocouple wells. 

Figure 20 50 

Probe thermocouple. 

Figure 21 53 

Photograph - Upper (left) and lower (right) heat 
exchanger manifolds. 

Figure 22 5^ 

Photograph - Welding of heat exchanger colls to 
lower manifold. 

Figure 23 55 

Photograph - Bottom view of heat exchanger. 

Figure 24 56 

Photograph - Heat exchanger positioned in vessel 
( top view) . 

Figure 25 57 

Photograph - Vessel positioned on shroud assembly. 

Figure 26 58 

Photograph - Vessel (with side heaters attached) 
showing insulation and vertical heated shroud 
partially complete. 

Figure 27 59 

Photograph - Subscale unit enclosed in inner 
(heated) shroud. 

Figure 28 60 

Photograph - Subscale unit enclosed in outer shroud. 

Figure 29 63 

Test circuit schematic. 

Figure 30 65 

Photograph - Major components and piping for 
test system. 


vil 



Figure 31 ^7 

Main T-66 pump control circuit schematic. 

Figure 32 73 

Thermkeep enthalpy from C&W cooldown data 

Figure 33 78 

T-66 cutlet temp, discharge test 0l4 (1 gpm) . 

Figure 3^ 78 

T-66 outlet temp, discharge test 012 (2 gpm). 

Figure 35 79 

t- 66 outlet temp, discharge test 003 (2 gpm). 

Figure 36 79 

T-b6 outlet temp, discharge test 007 (2.5 gpm). 

Figure 37 

T-66 outlet temp, discharge test 009-1 (2.8 gpm). 

Figure 38 80 

T -06 outlet temp, discharge test 001 (2.83 gpm). 

Figure 39 

T -66 outlet temp, discharge test 005 (3.5 gpm). 

Figure ^0 8 l 

t -66 outlet temp, discharge test 019 ( 6.3 gpm) . 

Figure 4l 82 

Thermkeep temp, profile, discharge test 01^< (1 gpm), 
at 3 hours . 

Figure 42 82 

Thermkeep temp, profile, discharge test 012 (2 gpm), 
at 2 hours . 

Figure 43 83 

Thermkeep temp, profile, discharge test 003 (2 gpm), 
at 3«5 hours. 

Figure 44 . 83 

Thermkeep temp, profile, discharge test 007 (2.5 
gpm, at 2 hours. 

Figure 45 84 

Thermkeep temp, profile, discharge test 009 -I 
( 2.79 gpm), at 2 hours. 

Figure 46 84 

Thermkeep temp, profile, discharge test 001 
( 2.83 gpm), at 2.33 hours. 

Figure 47 85 

Thermkeep temp, profile, discharge test 001 
( 2.83 gpm), at 5.33 hours. 


vlli 



Figure 85 

Thermkeep temp, profile, discharge test 005 
(3.5 gpm) , at 2 hours. 

Figure 49 86 

Thermkeep temp, profile, discharge test 019 
( 6.3 gpm), at 2 hours. 

Figure 50 86 

T-66 outlet temperature, charge test 015 (1 gpm). 

Figure 51 87 

T-66 outlet temperature, charge test 013 (2 gym). 

Figure 52 87 

T-66 outlet temperature, charge test OO 6 (2.5 gpm). 

Figure 53 • 88 

T-66 outlet temp., charge test 002 (2.83 gpm). 

Figure 5^ 88 

T-66 outlet temp., charge test 00^ (3.5 gpm). 

Figure 55 89 

Thermkeep temp, profile, charge test 015 
(1 gpm), at 2 hours. 

Figure 56 89 

Thermkeep temp, profile, charge test 013 
(2 gpm) at 4 hours. 

Figure 57 90 

Thermkeep temp, profile, charge test 6 
( 2.5 gpm) at 2 hours. 

Figure 58 90 

Thermkeep temp, profile, charge test 002 
( 2.83 gpm) at 3 hours. 

Figure 59 91 

Thermkeep temp, profile, charge test 004 
( 3.5 gpm) , at 6.5 hrs . 

Figure 60 94 

Initial Thermkeep temp, profile, idle test 010. 

Figure 6l 94 

Thermkeep temr . profile, idle test 010, at 25 hrs. 

Figure 62 95 

Thermkeep temp, profile, idle test 010, at 45 hrs. 

Figure 63 95 

Thermkeep temp, profile, idle test 010, at 65 hrs. 

Figure 64 96 

Thermkeep temp, profile, idle test 010, at 115 hrs. 


1 X 



Figure 65 96 

T-66 outlet temp., cyclic test 017 (2 gpm), cycle 1 

Figure 66 97 

T-66 outlet temp., cyclic test 017 (2 gpm), cycle 2 

Figure 67 97 

Thermkeep temp, profile, cyclic test 017 (2 gpm), 
at end of first discharge phase. 

Figure 68 98 

Thermkeep temp, profile, cyclic test 017 (2 gpm), 
at end of first charge phase. 

Figure 69 98 

Thermkeep temp, profile, cyclic test 017 (2 gpm), 
at end of second discharge phase. 

Figure 70 99 

Thermkeep temp, profile, cyclic test 017 (2 gpm), 
at end of second charge phase. 

Figure 71 99 

T-66 outlet temp., cyclic test 018 (2 gpm), cycle 1 

Figure 72 100 

T-d 6 outlet temp., cyclic test 018 (2 gpm), cycle 2 

Figure 73 100 

T-66 outlet temp., cyclic test OI8 (2 gpm), cycle 3 

Figure 7^ 101 

Thermkeep temo. profile, cyclic test 0l3 (2 gpm), 
at end of first discharge phase. 

Figure 75 101 

Thermkeep temp, profile, cyclic test OI8 (2 gpm), 
at end of first charge phase. 

Figure 76 102 

Thermkeep temp, profile, cyclic test 018 (2 gpm), 
at end of second discharge phase. 

Figure 77 102 

Thermkeep temp, profile, cyclic test 018 (2 gpm), 
at end of second charge phase. 

Figure 78 103 

Thermkeep temp, profile, cyclic test 018 (2 gpm), 
at f .id of third discharge phase. 

Figure 79 103 

Thermkeep temp, profile, cyclic test 018 (2 gpm), 
at end of third charge phase. 

Figure 80 107 

T-66 outlet temp, for 6/10 scale solar cycle 
simulation test. 


X 



Figure 8l 108 

Thermkeep starting temperature, cycle 1. 

Figure 82 108 

Thermkeep temp, at end of phase 1 discharge, 
cycle 1. 

Figure 83 109 

Thermkeep temp, at end of charge, cycle 1. 

Figure 8^1 109 

Thermkeep temp, at end of phase 2 discharge, 
cycle 1. 

Figure 85 110 

T-66 flow rate during charge. 

Figure 86 Ill 

Thermkeep temp, at starting profile, cycle 2. 

Figure 8? Ill 

Thermkeep temp, at end of phase 1 discharge, 
cycle 2. 

Figure 88 112 

Thermkeep temp, at end of charge, cycle 2. 

Figure 89 112 

Thermkeep temp, at end of phase 2 discharge, 
cycle 2. 

Figure 90 113 

T-66 flow rate variation during second charge cycle. 

Figure 91 11^ 

T-66 outlet temp, for full-scale solar cycle 
simulation test. 

Figure 92 115 

Thermkeep temp., test 023 starting profile, cycle 1 

Figure 93 115 

Thermkeep cemp. at end of phase 1 discharge, 
cycle 1. 

Figure 9^ 116 

Thermkeep temp, at end of charge, cycle 1. 

Figure 95 116 

Thermkeep temp, at end of phase 2 discharge, 
cycle 1. 

Figure 96 117 

T-66 flow rate variation during first charge cycle. 

Figure 97 II 8 

Thermkeep temp, at starting profile, cycle 2. 


xi 



Figure 98 118 

Thermkeep temp, at end of phase 1 discharge, 
cycle 2. 

Figure 99 119 

Thermkeep temp, at end of charge, cycle 2. 

Figure 100 119 

Thermkeep temp, at end of phase 2 discharge, 
cycle 2. 

Figure 101 120 

T-66 flow rate variation during second charge cycle 

Figure 102 12^1 

Specific cost of storage for various heat exchanger 
surface areas at constant Thermkeep masses. 

Figure 103 125 

Figures of merit for various Thermkeep masses at 
constant heat exchanger surface areas. 

Figure 104 126 

Figures of merit for various heat exchanger 
surfaces at constant Thermkeep masses. 

Figure 105 127 

Specific cost of storage for various Thermkeep 
masses at constant heat exchanger surfaces. 

Figure 106 123 

Specific cost of storage for various Thermkeep 
masses at constant figures of merit. 

Figure 107 129 

Specific cost of storage and figures of merit for 
various tube diameters. 

Figure 108 13f 

Specific cost of storage, peak pump power, and 
figures of merit at various heat exchanger pres- 
sure drops. 

Figure 109 131 

Specific cost of storage and figures of merit at 
various aspect ratios. 

Figure 110 132 

Effect of insulation thickness on cost and 
performance. 

Figure 111 133 

Specific cost of storage and figures of merit 
at various peak pump flows. 

Figure 112 134 

Specific cost versus deficit heating. 

Figure Cl C-3 

Assumed phase diagram NaOH-NaNO^ 


xll 



SUMMARY 


Anhydrous sodium hydroxide modified by addli;lon of 85? 
sodium nitrate and 0.25? manganese dioxide has been under 
development as a phase-change heat storage medium for a 
number of years for application to electrically powered, 
low temperature space and water heating. This medium 
Is low In cost, has a high heat storing capacity, and Is 
compatible with steel containers and heat exchangers. 
Therefore, It Is of Interest for use In solar powered 
electricity generating systems. 

The purpose of this program was to make an analytical 
and experimental study of a phase-change thermal stor- 
age system utilizing this medium, suitable for use In a 
solar total energy system, and to develop a preliminary 
design of such a storage system consistent with the fol- 
lowing specifications: 

Storage capacity is 3-1 x 10^ kJ operating 
over the range of 516 K to 584 K. The heat 
transfer fluid for charging and extracting 
heat is Thermlnol-66 . The maximum charging 
rate Is 1.8 x 10° kJ/hr with the Thermlnol-66 
temperature at 584 K ± 2 K,,and the maximum 
extraction rate is 1.0 x 10° kJ/hr with the 
Thermlnol-66 Inlet temperature 516 K ± 11 K 
and outlet temperature 582 K ± 2 K. 

The analysis showed that supplemental heat of about 10? 
of the heat extracted from storage, added to the Thermi- 
nol-66 during discharge to help maintain the outlet tem- 
perature within the specified range, had a large oeneflclal 
effect on the size and cost of the storage system. There- 
fore such supplemental ("deficit”) heating was employed 
in the design of the system. 

A mathematical simulation of the system and computer pro- 
gram were available as a result of prior work. These were 
amended to facilitate use in the present work. 

The program comprised the following tasks: 

. . . Measurement of physical and thermophysical 
properties of the storage medium. 

... Adaptation of the previously developed simu- 
lation and computer program to present re- 
quirements. 



. . . Development of a reference design of r. sys- 
tem meeting the operating requirements. 

. . . Deslgr. and construction of an experimental 
model of 1/10 scale, and a test-bed In 
which to test It. 

. . . Operation of the 1/10 scale model under test 
conditions consistent with the specifica- 
tions, and collection of experimental data. 

. . . Correlation of data with computer-predicted 
results . 

... Modification of the simulation and computer 
program until correspondence was acceptable. 

. . . Development of a preliminary design meeting 
the operating specifications and estimation 
of the system cost. 


Satisfactory correspondence between experimental data and 
computer predicted results has been obtained, and the pre 
llmlnary design developed. The storage unit consists of 
a cylindrical vessel, surrounded by thermal Insulation, 
containing the storage medium In which Is immersed a heat 
exchanger consisting of a number of helical colls of 
steel tubing, manifolded for parallal flow at the top and 
bottom, through which the Therminol -66 flows downward dur 
Ing charging and upward during heat extraction. 


The preliminary design summary Is as follows: 


Vessel height and diameter 

Insulation thickness 

System helirht and diameter 

System weight 

Amount of storage medium 

Unit cost 

Heat exchanger 

Storage output 
Specific cost 
Insulation loss 

Deficit heating 


2.35 m (92.5”) 

0.61 m (2il") 

3.6 m (l.i»2") 

27300 kg ( 60,200 lb) 
18000 kg ( 39,700 lb) 
$33,700 

.00635 m ( 1 /^") tub- 
ing; 280 coils, each 
27.7 m (91 ft) long 

3..:i X 10^ MJ 

$11.2/MJ 

3.25/5 of total stor- 
age output per day 

10.2% of storage out- 
put 
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INTRODUCTION 


In the generation of electric power from solar heat, by 
means of a heat engine, storage of heat is necessary 3n 
order to make the power generation system Independent of 
the solar insolation. The heat storage element of such 
a system is one of the more costlj' elements of the system 
and the cost goals for thermal storage have not yet been 
achieved. 

One approach to thermal storage for solar-electric sys- 
tems which offers a possibility for reducing the cost 
of storage is the use of thermal storage media which 
undergo phase changes between solid and liquid states, 
or between solid states, within the operating temperc-ture 
range, making the latent heats of the phase changes, as 
well as sensible heat, available for heat storage and 
retrieval. 


Objective 


The objective of the work reported herein has been to 
study the dynamic performance of a phase change thermal 
storage unit employing a passive heat exchanger, suitable 
for use in a solar total energy system, operating in the 
mid-temperature range, to develop a preliminary design 
for such a unit, and to estimate its cost. 


Concept for Phase-Change Thermal Storage 


The concept for the use of the phase change medium, which 
has been studied in this work, comprises a vessel contain 
ing the medium in which is immersed a single fixed heat 
exchanger, through which a single heat transfer fluid is 
passed for charging heat into storage when heat is avail- 
able from the solar collectors, and for retrieving heat 
from storage as required at other times. 

The particular configuration studied consists of a cylin- 
drical steel vessel, and a heat exchanger consisting of 
a nuL.jer of helical coils of steel tubing which are mani- 
folded together, inside the vessel, at the top, and at 
the bottom, forming parallel flow paths for the heat 
transfer fluid. When heat is being charged into storage, 
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the fluid flows downward through the heat exchanger, and 
upward when stored heat is being retrieved. 

A thermal gradient exists within the medium (with the 
higher temperature In the upper part of the vessel ) which 
moves vertically and changes shape during thermal cycling. 

The medium undergoes an Increase in volume as it melts, 
and to accommodate this a clearance space must be allowed 
at the top of the vessel. This space is open to the 
atmosphere through a "breather" tube which allows air to 
enter and leave during cycling, and which insures that 
the vessel operates unpressurlzed. The vessel is sur- 
rounded by thermal insulation, the outer surface of which 
is protected by a metallic shroud. 


Heat Storing Medium 


The heat storing medium used in the work repoi ted herein 
is tradenamed Thermkeep* and has the following composition 

Anhydrous NaOH (commercial grade) 91.8? by wt. 

Sodium nitrate 8.0? 

Manganese dioxide 0.2? 

The commercial grade of NaOH commonly contains 1-2? of 
sodium chloride, 1/2-1? of soaium carbonate, and smaller 
amounts of other salts and hydroxides. 

Salt baths of similar composition sold under several 
tradenames have been used Industrially for many years at 
tempei atares of 750 K and higher for metal cleaning and 
descaling. These baths are contained in large steel 
tanks. When the advantages of anhydrous NaOH as a phase 
change heat storage medium were recognized, this compo- 
sition was used because of this existing Industrial 
experience with the containment and stability of the 
material. 

As a heat storing medium, Thermkeep has the advantages of 

... high heat storing capacity with latent 
heats of fusion and a solid state phase 
change in the temperature range 566 K 
to 507 K, 


* Registered trademark Comstock & Wescott, Inc. 




. .. chemical stability and insignificant vapor 
pressure at temperatures up to 750 K, 

. . . s'cability in contact with air allowing the 
venting of containment vesse'’s to the at- 
mosphere, 

. . . compatibility with steel for containment 
and heat exchangers up to 750 K, 

... low cost and ready availability, and 

. . . large scale industrial experience as a 
metal cleaning salt. 


Background 


The use of Thermkeep as a high temperature phase change 
thermal energy storage (TES) medium has been under 
development by Comstock & Wescott, Inc. since 1958. The 
initial proposed use for TES units was for domestic space 
heating in conjunction witi heat pumps. Off-peak electric 
energy was to be used to chcrge the TES units which would 
assist the heat pumps during periods of high heating re- 
quirement. This development, sponsored by the Edison 
Electric Institute and the Philadelphia Electric Co. was 
carried through a successful field test from 1962-196^, 
but was then discontinued because the anticipated large 
scale use of heat pumps did not materialize. 

Beginning in 1964, in cooperation with the Hooker Chemicals 
and Plastics Corp., Comstock & Wescott undertook further 
development of the heat storage material and development 
of equipment for several anticipated applications of 
thermal storage. Extensive studies were carried out on 
chemical formulations and the compatibility of Thermkeep 
with containment materials, which has resulted in nearly 
complete elimination of corrosion of containers and heat 
exchangers fabricated from ordinary mild steel at tempera- 
tures up to 755 K (900 F) . 

A result of the cooperative effort between Comstock & 
Wescott and Hooker Chemicals and Plastics Corp. has been 
the development of the Thermbank Electric Water Heater, 
which was designed to provide hot water for all-electric 
buildings. Five Thermbanks were successfully carried 
through a one-year field test in a cooperative program 
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with the Hydro Electric Power Commission of Ontario. A 
Thermbank Electric Water Heater of this type has been 
supplied to the NASA Lewis Research Center under Contract 
No. NAS 3-202^10 for in-house testing and evaluation. 

The applications described above used electric energy to 
charge heat into storage, stored it at temperatures up 
to 755 K (900 F), and delivered it at approximately 3^0 K 
(150 P) for space and wate. heating. In these applica- 
tions the temperature of the heat delivered from stor- 
age can be allowed to drop substantially below the maxi- 
mum storage temperature as the thermal energy store is 
depleted without affecting the system efficiency. 

When heat from solar collectors is stored and retrieved 
from storage for power generation, a drop in temperature 
of the heat supplied to the heat engine reduces the effi- 
ciency of the power generating cycle. Therefore, the 
heat from storage must be delivered to the heat engine 
at a temperature close as possible to that at which it 
is charged into storage. This requirement for delivery 
of stored heat at a nearly uniform temperature imposes 
a new set of conditions on the design of the heat ex- 
changers, and experience with heat exchangers designed 
for the previously developed devices does not apply. 


Scope of Work 


The study has Included both analysis and the design, 
construction, and testing of a sub-size ■'■hermal storage 
unit for the acquisition of experimental data for cor- 
relation with the analysis. The analysis was based in 
part upon the results of prior work by Comstock & WescoLL 
under Purchase Order No. 87-5030 from Sandia Laboratories, 
Livermore, CA, to develop a mathematical model and write 
a computei^ program describing the dynamic behavior of a 
phase change thermal storage system using Thermkeep as 
the medium. In addition, the work has included the 
measurement of certain thermophysical properties of the 
medium. 

Specifically, the program has included the following 
tasks : 

1. Experimental determination of the effect 
of accidental contact between Therminol-66 
(heat transfer fluid) and Thermkeep at 
temperatures up to 613 K. 
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2, Determination of physical properties of Therm- 
keep: 

a. Thermal conductivity of the solid phase 
from room temperature to the melting 
point, and of the liquid phase from the 
melting point to 750 K. 

b. Density of the solid from room tempera- 
ture to the melting point, and of the 
liquid to 750 K. 

c. Heat capacity of the solid and liquid to 
750 K, including the heat of transition 
from one phase to the other. 

d. Viscosity of the liquid from the melting 
point to 750 K. 

3. The development of a Reference Design, based 
upon the previously developed Phase Change 
Computer Program (suitably modified for this 
purpose), for a phase change thermal energy 
storage (TES) unit suitable for use in a solar 
electrio power generation system, utilizing 
Therminol-66 as the heat transfer fluid for 
charging and retrieving stored heat . The TES 
unit has the following operating requirements: 

Storage capacity is 3.1 x 10^ kJ 
operating over the range of 516 K to 
58i< K. The heat transfer fluid for 
charging and extracting heat is Thermi- 
nol-66. The maximum charge rate is 
1.8 X 1Q6 kJ/hr with the Therminol-66 
temperature at 58* ** K ± 2 K, and the. 
maximum extraction rate is 1.0 x lo6 
kJ/hr with the Therminol-66 outlet 
temperature at 582 K ± 2 K and inlet 
temperature at 516 K i 11 K. 

*1. The design and construction of a subscale TES 
system, of one-tenth the scale of the Reference 

Design. The subscale system Includes the TES 
unit and a test bed which permits testing under 
operating condition., consistent with the opera- 
ting characteristics of the Reference Design. 
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5. The operation of the subscale TES system 
under a range of conditions and the collec- 
tion of experimental data. 

6. Comparison of the experimental data with 
results predicted by the analytical model 
via the computer program, and revision 
until satisfactory agreement was attained. 

7. A Preliminary Design of a full scale TES 
unit and estimation of the cost of this 
unit . 
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TECHNICAL DISCUSSION 


MATERIAL PROPERTY STUDIES 


A subcontract was Issued by Comstock & Wescott, Inc. to 
Dynatech R/D Company, Cambridge, Mass., to determine 
some relevant properties of the storage medium, Thermkeep. 
This was in conjunction with evaluating the feasibility 
of using Thermkeep in thermal storage devices for solar 
energy applications. Dynatech was selected because of 
the company’s experience in measuring properties of in- 
organic heat-of-fuslon heat storage materials (Ref. 1). 

A copy of Dynatech 's final report, "The Thermophysical 
Properties of a Thermal Energy Storage Material," dated 
January, 1978, is included as Appendix A. 

The work had two objectives: 

1. To determine the consequences of con- 
tact between the heat transfer fluid, 

Therminol-66 (T-66), and Thermkeep, 
and 

2. To measure specified thermophysical 
properties of Thermkeep. 

Regarding the compatibility of Thermlnol-66 and Therrn- 
keep, the Dynatech analysis showed that no detrimental 
reactions occur between the two materials. The majoj* 
effect occurred in the Therminol-66. A gradual increase 
in its viscosity was observed when the Therminol-Thermkeep 
mixture was in contact with air at elevated temperatures 
(593 K; 320 C). Dynatech recommended the use of a cover 
gas, such as dry He, over the Thermkeep to remove any 
water or T-66 which may accumulate in the Thermkeep due 
to a leak in the heat exchanger. This is not necessary 
from an operational standpoint and neglects the fact that 
air breathing is required to maintain the proper balance 
of NaNO^ which otherwise is slovuly reduced to NaNOg while 

the unit operates. 

The thermophysical property measurement phase required 
determination of the following properties of Thermkeep 
as a function of temperature: 
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1. The enthalpy (or specific heat). 

2. The thermal conductivity of the solid 
and liquid phases. 

3. The density in both phases. 

The viscosity of the liquid. 

Data are provided in the report on all but the enthalpy 
of the material. It was found that the Quantitative 
Thermal Analyzer (QTA) used for accurate calorimetric 
study of materials, owing to its design, would not provide 
accurate data in the temperature region of interest. 
Consequently, the data utilized for enthalpy vs. tempera- 
ture in the course of this work was based upon the best 
information which Comstock & Wescott had prior to this 
project. It is felt that at some point this ought to be 
Improved upon but correlation of analysis with experi- 
ment has shown that the present data are reasonably good. 


10 



REFERENCE DESIGN STUDY 


Computer Model 


A mathematical analysis and computer program were writ- 
ten to estimate the performance of a Thermkeep thermal 
energy storage device to meet the requirements of the 
Sandla Midtemperature Solar Test Facility and to design 
the one-tenth scale test model. It was expected to re- 
quire modification for correlation with actual data 
obtained from tests on the model. It was formulated as 
a design tool which would properly account for the thermo- 
dynamics of the process, approximately account for the 
heat transfer between the T-66 and the Thermkeep, and 
allow for easy handling of a variety of operating modes, 
from solar dally cycles to laboratory test modes. 

It is an analytical program rather than a design program 
in the sense that it requires a specific unit design as 
input, and outputs the calculated performance of that 
unit. This program is used for design by analysis of a 
large matrix of design options and separate parametric 
evaluation of the output by the engineer. 

The physical processes associated with the heat storage 
medium are quite complex and would be difficult, if at 
all possible, to model even given a large dedicated 
effort. The medium undergoes a solid-liquid phase change 
with a volume contraction associated with solidification. 

It is a non-eutectic multi-component material so that 
solidification, in part, probably occurs by precipitation 
of NaOH from the melt when the solubility is below the 
actual concentration. A solid-solid phase transformation 
occurs in the NaOH with a corresponding latent heat. 

Note: All symbols are defined in the Compute;- Analysis 

Symbol List presented as a fold-out in Appendix D. 

Another computer analysis was formulated earlier in a 
prior program to attempt to account for some of these 
effects. It is a considerably more costly program to use; 
it is less flexible, and it would not as easily lend it- 
self to the incorporation of empirical modifications. It 
is also based, in large part, upon a published phase diagram 
for the system Na 0 H-NaN 03 . It is not known, however, 
whether or not this is truly characteristic of Thermkeep 
which is composed of commercial grade chemicals contain- 
ing impurities as well as the basic two components above. 
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Consequently, it was decided that a simpler analysis, 
of the type to be described herein, would be more de- 
sirable. This earlier analysis is described In Appen- 
dix C. 


General Description 


The computer analysis Is formulated around the design 
concept of the proposed thermal storage unit. The design 
concept consists of a vertical cylinder containing the 
heat storage medium. A heat exchanger is immersed In 
the Thermkeep. The heat transfer fluid, T-66, flows 
through the heat exchanger tubing. (Any fluid which does 
not undergo vaporization can be modeled using its particu- 
lar thermophysical properties.) The heat exchanger con- 
sists of a number of parallel heat exchanger tubes which 
are identical geometrically and assumed to pass equal 
flows of T-66. The tubes are manifolded at the top and 
bottom of the unit, the hot end being at the top. Dur- 
ing charging, T-66 enters at the top and exits from the 
bottom, and vice versa during discharging. 

It is desired to know how the unit will respond to par- 
ticular stimuli, e.g., flow rate, temperature, and flow 
direction of T-66, both internally and as a "black box" 
characterized by what the outside world sees. The analy- 
sis has been constructed with two major sections. One 
contains the analysis of the heat transfer between the 
T-66 and medium within the unit. The other characterizes 
the stimulus and determines the status of the T-66 deliv- 
ered to the unit as a function of the operating time. 

This latter section is modified at will to allow modeling 
of various modes of operation. The mode described herein 
is the solar dally cycle which was used in developing 
the approximate Reference Design described in a subsequent 
section . 

Figure 1 describes the solar daily cycle which was taken 
to be typical. of the mode of operation of the Sandia Test 
Facility. The solar collectors generate heated T-66 which 
is supplied to the boiler of a toluene based Rankine power 
cycle. The boiler demands 250 kw for twelve hours a day. 
During the other twelve hours the system is shut down. 

The solar collectors produce a varying heat output over 
nine hours, peaking at 500 kw, starting at zero in the 
morning and dropping to zero in the late afternoon. The 
excess heat is stored to allow operation when the collec- 
tor output is deficient. 
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Figure 1. Typical solar daily cycle at Sandia Total Inergy Test Facility. 


During charging of the unit, the flow rate to the stor- 
age unit Is determined by the available excess solar 
heating, and the temperature change of the T-66 as It 
passes through the unit: . 


w 


f 


iex 

iT^rir 
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where w_ Is the T-66 flow rate to the storage unit, q 

I ex 

Is the excess thermal power from the collectors (l.e., 

above boiler demand), c^ Is the T 66 specific heat, and 

Tjj and are the T-66 Inlet and outlet temperatures, 

respectively, of the TES unit, w^ Is limited to some 

realistic value of pump capacity. 

During discharge of the unit, the flow rate Is held con- 
stant as If the fluid were heated over the full tempera- 
ture range. It will. In fact, exit somewhat cooler than 
the charging temperature. In this case, a boost heater 
Is assumed and the amount of this heating required is 
calculated. 

The unit begins at a fully charged condition and is per- 
turbed by the solar daily cycle. The calculation pi^o- 
ceeds until a repetitive cyclic operation Is achieved 
and the unit is characterized by Its steady static cyclic 
performance under this mode of operation. Cumulative 
heat balances over the cycle are made to define the per- 
formance. The important factors are the dally heat ab- 
sorption by the unit, the dally heat delivery, the daily 
heat loss to the environment, and the daily make-up heat- 
ing (boost heating) requirement. 

Based upon the best knowledge available at the time, a 
method was devised for estimating cost of materials and 
fabrication for a given design. These costs are computed 
at the beginning of each analysis for the particular de- 
sign configuration. The costs are printed out In an 
Itemized fashion for each part of the assembled unit and 
lead to a total unit cost. 

Further details of the computer model used to develop the 
Reference Design are given in Appendix B. 
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Reference Design Analysis 


The computer analysis which was described above and In 
Appendix B was used to produce a ’'reference design" for 
a Thermkeep thermal energy storage unit. This unit was 
specified as being subjected to the solar daily cycle as 
described on Figure 1. The unit was designed to provide 
a reasonable compromise between cost and efficiency as 
will be explained later. Aside from the heat flow speci- 
fications, the fluid utilized was Therminol-66. It was 
assumed that during charging the fluid would enter stor- 
age from the solar collector field at 584 K (311 C) and 
that during discharge It would enter storage from the 
boiler at 516 K (243 C). 


The design selected was to become the basis for the de- 
sign of a laboratory model of one-tenth scale in size, 
heat flows, fluid flows, etc., but with equivalent per- 
formance In terms of outlet temperatures, losses, and 
other factors. Eventually, after adjustments to the com- 
puter analysis had been made based upon empirical cor- 
relations, a new final reference design would be composed 
which would more accurately represent how a unit for 
Sandla Laboratories STETP would appear. 

Note: At the point In the project when the ref- 

erence design was composed, the available 
specific energy over the operating range 
was thought to be 450 kJ/kg. Subsequent 
investigation showed this to be closer to 
4?0 kJ/kg. However, the results to be 
discussed herein were based upon the 
larger value. They are shown to describe 
the method and trends but a roughly 1055 
addition to cost and size would be expected. 

The study was not rerun since the model was 
already In fabrication when the correction 
was discovered. The later analyses utilize 
the correct enthalpy vs. temperature relation- 
ship. rtny graphs and tables presenting re- 
sults u-iing the lower value are so noted. 

The procedure for design was as follows. Preliminary 
screening of some of the parameters had shown that the 
primary variables are the quantity of storage medium and 
the total heat exchanger surface area. Thus, data were 
generated for a specific tube size and a preliminary ref- 
erence design was produced. The secondary variables were 
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then tested to determine where Improvements could be made 
to arrive at a reference design. 

The unit Is presumed to be a shell and tube device with 
helical'' y colled parallel tubes with a specified tube 
size, number and length. The variable parameter was the 
total Inside surface area of the tubes. Heat exchangers 
were sized to produce about 703 g/cm^ (10 psl) peak pres- 
sure drop as a starting point. The reference tube size 
was 6.35 mm (1/4") o.d. steel tubing, with 4.57 mm ( 0 .l 80 ”) 
l.d. The ”atlo of vessel (sheJl) height to diameter, l.e-, 
the "aspect ratio" was fixed at 1 for screening. During 
cooling, flow enters the bottom end or element 1, anr 
vice versa during heating. 

For computational purposes, the solar collector and boiler 
demand characteristic as described wa^ Input. A suffi- 
cient number of points were Input and straight lines were 
assumed between them. This characteristic was adapted 
from the "winter cycle" of the Sandla Laboratories and 
provides peak collector output of 500 kV/, steady boiler 
demand of 250 kW and more excess heating available dur- 
ing charging than required during discharging. When 
steady state cyclic operation is achieved, the excess Is 
simply wasted. The total storage demand Is 3.45 x 10° kJ 
by this cycle, not 3.1 x 10° as specified in the contract, 
but a simple scaling could reduce this to the proper value. 

A 12-hour Idle period is specified which is consistent 
with Sandla Laboratories' mode of operation. This creates 
unusual situation where the capacity is low but the 
ind rate Is high since It occurs for a short period 
v-lme. It Is expected that a unit sized to operate 
, .roughout the idle period would f.'^^’form better over the 
complete cycle. 

The computer analysis was run over enough cycles to stabi- 
lize the unit performance. A criterion was found accept- 
able whereby the unit cycles until the storage cyclic 
InpiT minus the storage cyclic output is close to the nomi- 
nal stand-by or surface loss. This condition is approached 
by ftartlng with the unit fully charged to peak fluid tem- 
perature and then cycled until the stable full charge 
gradient is reached. 

Flow rate modulation during charg.ing was provided based 
upon fluid outlet temperature to attempt to maintain 
heat balance around the solar collectors. A maximum 
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total flow of 5 kg/sec was specified to collectors and 
boiler. In a thermocllne system, 2.95 kg/sec would be 
required at 500 kW output so that an increase of 70$ 
is allowed — a value which is adjustable. 


Figures 2 and 3 describe the results of analysis plotted 
in two different ways. Figure 2 being constant surface 
area lines with amount of Thermkeep as the coordinate 
axis. Figure 3 being the opposite. Overall efficiency is 
defined as 


where 



Q 

+ Qg + Q, 


= storage dally output 
Q = storage dally surface loss 
Q = fluid auxiliary heating 

a 


Q is computed as the cyclic amount of heat required to 

cl 

bring the fluid up to 584 K (311 C), the peak temperature. 

In fact, the specification allows as low as 580 K (307 ('} 

so t’nat the overall efficiency is somewhat pessimistic 

as will be discussed later. Q is based upon 45.7 cm (iR") 

s 

of rock wool Insulation. In application, Q and Q would 

S 3 

be weighted differently depending upon their effect on 
total system economics. Without knowledge of how to 
weight them accurately, the above definition of efficiency 
was used. It is not a pure thermodynamic efficiency but 
was felt to be adequate when comparing designs subjected 
to the same ducy cycle. 

Surface areas ranged from 37.16 sq m (400 sq ft) to 185.8 
sq m (2000 sq ft) in parallel tubes sized for 69 kPa 
(10 psi) peak pressure drop as discussed. Thermkeep quan- 
tities were varied from 6000 kg to 50000 kg. Figure 3 
shows the normal improvement expected with Increasing 
surface area with no significant increase in performance 
beyond 148.6 sq m (1600 sq ft). Figure 2 shows en un- 
expected decrease in performance with large amounts of 
storage medium. This is believed to be due to the forma- 
tion of permanent solid throughout the jnit with inade- 
quate surface to extract heat from it. The calculation 
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is believed to be pessimistic for solid heat transfer 
and real performance may be somewhat better. 

Cost factors were developed whereby materials costs were 
estimated at market values and fabrication costs were 
scaled from the Mod I Thermbank costs according to the 
normal "six-tenths" factor whereby fabrication costs 
increase as the 0.6 power of the size ratio. This unit 
is the accepted rule for process industry equipment such 
as heat exchangers, vessels, etc. 

Figures 4 and 5 show the costs of storage in dollars per 
thousand kilojoule based upon the resultant storage out- 
put which is taken to be the unit capacity. The important 
criterion, however, is not simply cost of storage but cost 
at an acceptable efficiency. Figure 6 is cross plotted 
from Figures 2 through 5 taking cost of storage for con- 
stant overall efficiency. From these curves and from 
some interpolation, a preliminary reference design was 
selected which uses 15000 kg of Thermkeep and 83.6 sq m 
(900 sq ft) of surface area. The cost was estimated at 
$27»880 with a total storage output of 3-225 x 10° kJ 
for $8.64/MJ. a total efficiency of 0.911 was predicted. 

All data were obtained with only 10 numerical elements. 

It was determined that the grid size, beyond 10 or so, 
had a small effect on the computed performance. Figure 7 
shows the stable temperature swings for a 10 element grid 
and Figure 8 shows these for a 25 element grid. Using a 
finer grid Improves the computed performance slightly to 
a total output of 3-237 x 10° kJ for $8.6l/MJ and an 
overall efficiency of .91^. Overlaying the two tempera- 
ture maps shows them to be similar except that the 25 
element grid has more definition (high frequency resolu- 
tion) as expected. T1 through t 6 are key time' as de- 
scribed on Figure 9, 

The major effect of the 12-hour stand-by period (T2 to 
T3) is a loss of high temperature heat through the top 
of the unit. This heat is mainly sensible heat of the 
liquid phase above the melting range of the NaOH-8? 

NaNO^ mixture (see Figure Bl). This aspect of perform- 
ance (and others) might be improved by decreasing the 
NaNO^ content and thereby raising the melting range of 

the mixture, so that the high temperature heat loss would 
come principally from latent rather than sensible heat. 
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Overall Efficiency 



Amount of Thermkeep, Mg 


Figure 2. Overall efficiency as a function of 

Thermkeep mass, showing effect of heat 
exchanger surface area. 
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Overall Efficiency 


Thermkeep 
Mass. Me 



Figure 3- Overall efficiency for constant Thermkeep mass 
at various heat exchanger surface areas. 




Figure 4. Cost of storage for constant heat exchanger 
surface area at various Thermkeep masses. 




Figure 5. Cost of storage for constant Thermkeep mass 
at various heat exchanger surface areas. 
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Cost of Storage, $/MJ 



Figure 6. Specific cost of storage for various TES 
designs at fixed overall efficiencies. 
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storage Temperature, deg K 



Element No. 


Figure 7. Reference design with 10 element grid 
(preliminary) . 
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Storage Temperature, deg K 



Figure 8. Reference design with 25 element grid 
(preliminary) . 
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Collector Output, kW; Boiler Demand, kW 
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Figure 9. Reference duty cycle. 


The fact that pure liquid with low specific heat exists 
above 578 K (305 C) is a contributing factor for the high 
temperature drops in this region during extraction of heat 
and during idle. In a unit design, it may be possible 
to use higher quality insulation only at the top and 
stainless steel outlet tubes to reduce loss. 

The preliminary reference design consists of I 5 OOO kg 
of Thermkeep; 252 parallel colled tubes, 0.635 cm o.d. x 

0. ^57 cm i.d. X 23.1 m long (0.25” o.d. x O.I 8 O'' l.d. x 
75*8 ft long); a vessel whose inside dimensions are 
2.23 m ( 7.31 ft) diameter by 2.23 m (7.31 ft) height. 

Total outside dimensions are 3*17 m (10.^ ft diameter 
by 3.19 m CIO. 5 ft) height. Total outlet capacity is 
3.237 X 10° kJ at 91 .^^5? overall efficiency. Of the 
8 . 56 ? inefficiency, about 2.5? is due to surface losses 
which could be Improved with high quality but expensive 
insulation. 

The next step was to evaluate the effects of minor param- 
eters. The preliminary reference design was used as a 
basis for all these parametric evaluations. Figure 1C 
shows the effect of increasing tube size. It can be 
shown from the cylindrical conduction equation that for 
equal surface areas and equal deposition of solid that 
small diameter tubes provide less conduction resistance 
so that the results were not unexpected. Clearly, the 
cost of storage Increases and the overall efficiency de- 
creases with increasing tube size. 0.635 cm (1/4") o.d. 
tubes are considered to be the smallest practical for fab- 
ricating the unit so that this size would be retained. A 
25 -element grid was used here but 10 elements were used 
for the rest of the studies. 

Figure 11 shows the effect of variation in aspect ratio 
(vessel height/vessel diameter). As aspect ratio in- 
creases for the same volume, surface losses Increase and 
axial conduction decreases. The results show decreasing 
performance with aspect ratio, suggesting that surface 
losses are more significant than axial conduction. In 
fact the decrease in efficiency can be attributed almost 
solely to increased insulation loss. From a pure surface 
loss standpoint, it can be shown a priori that an aspect 
ratio of 1 is an optimal configuration. The aspect ratio 
of 1 is therefore retained. 

Figure 12 shows the effect of increased fluid velocity, 

1. e., changing tube lengths and number to achieve the 
same surface area with higher velocities and consequently 
higher pressure drops and pumping powers. While there is 
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an expected performance Improvement with Increased veloc- 
ity, this must be weighed against increased pump power. 

If we assume that mechanical power in a total energy 
system is obtained at 10? efficiency, then in the extremes, 
viz., 76 kPa and 710 kPa (11 psl and 103 pel), the fol- 
lowing is true. It takes 57600 kilojoules-thermal to pro- 
duce the 5760 kilojoules-mechanical required to pump fluid 
for a day at 76 kPa (11 psl) and 551160 kllojoules-thermal 
to produce the 55116 kilojoules-raechanical at 710 kPa 
(103 psl), a net increase of ^93560 kllojoules-thermal. 

The additional 1? overall efficiency at 3.2 x 10° kllo- 
Joules-storage per day yields 32000 kilojoules-thermal 
which is not justified. Even from 76 kPa to 1^15 kPa 
(11 psl to 21 psl), an additional 52560 kilojoules-thermal 
is required to obtain 0.35? of 3.2 x 10° kilojoules or 
9600 kilojoules-thermal. Thus, the 70 kPa (11 psi) speci- 
fication was accepted. The peak pump power is now 0.33 
kW-ideal vs. 3.12 kW at 711 kPa (103 psi) and 0.64 kW at 
145 kPa (21 psi). Peak power may be reduced by reducing 
peak allowable flow which mav be worth considering. 

Figure 13 shows the effect of insulation thickness on 
the preliminary reference design. Performance improves 
with insulation thickness but with diminishing returns 
due largely to the cylindrical geometry. The increase 
in cost from 0.305 m (12") to 0.457 m (I 8 ") is $450 
which affords about 1?; the increase in cost from 0.457 m 
(lO") to 0.61 m (24”) is $520 which affords 0.5?; the 
increase in cost from 0 . 6 I m (24”) to O .76 m (30”) is 
$600 which affords 0.2?. The decision was made to go 
to 0.76 m (24”) of insulation but in reality a cost- 
effectiveness study for all parameters ought to be done 
which assesses the impact on the whole Total Energy System. 

Figure l4 shows the effect of peak allowable storage flow 
due to flow modulation during charging. Beyond a certain 
point, the excess flow doesn’t provide any gain since the 
unit can do no more than come to a full charge. Below 
that point inadequate flow modulation prevents insuffi- 
cient unit recharge. It would appear that something a 
bit smaller than the 5 kg/second total flow or 3.5 kg/ 
second storage flow would be tolerable but a rapid dete- 
rioration of performance eventually results. For now, 
the 5 kg/sec (3.5 kg/sec storage flow) will be retained. 

The final reference design is now essentially the prelimi- 
nary reference design with 0 . 6 I m (24”) of rock wool insu- 
lation. The estimated unit cost is $28,399 with a capacity 
of 3.24 X 10° kilojoules output, for $ 8,765 per thousand 
kilojoules. The total efficiency from a 25-element grid 
analysis is 0.919 and the insulation loss is 2.2? of th= 
total storage capacity. 
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Figure 10. Effect of tube size on cost and 
performance (25 element grid). 
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Figure 12. Effect of fluid velocity on cost 
and performance (10 element grid) 
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Efficiency 


Figure 13 . Effect of insulation thickness on cost 
and performance (10 element grid). 










Figure 15 is a map of the movement of the thermal gradient 
during stable cyclic operation, for a 25-element grid. 

The overall efficiency is based upon the insulation loss 
of 72,000 kilojoules and 213,500 auxiliary input to bring 
the fluid to 584 K (311 C). A modified efficiency can 
be utilized which involves bringing the fluid to 580 K 
(307 C), the specified minimum temperature, when it drops 
below this value. Then the modified auxiliary input be- 
comes 87,270 kilojoules and the modified overall effi- 
ciency becomes 0.953- Figure 16 is a graph of the fluid 
outlet temperature vs. cycle time during discharge. 

The configuration is now a vessel which is 2.23 m diam- 
eter by 2.23 m height (7 '4" diameter by 7 '4" height), 
containing 252 tubes, 0.635 cm o.d. x 0.4572 cm i.d. 

( 0 . 25 " o.d. X 0 . 180 " i.d.), and 15,000 kg ( 33.000 lb) of 
Thermkeep. The total overall dimensions including insu- 
lation are 3-474 m diameter by 3.50 ra height (' 11 * 5 " 
diameter x 11*6" height). The tubes would be constructed 
as individual identical helical spring-like coils, all 
parallel to each other along the vertical axis of the 
unit . 
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Element No. 

Figure 15. Final full-scale reference design — 
movement of thermal gradient. 
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Fip-ure l6. Final full-scale reference design — variation of 
Thermlnol-66 outlet temperature during discharge. 




SUBSCALE TES UNIT 


Following the completion of the reference design for the 
full scale unit, It was necessary to design a model thermal 
storage unit of about one-tenth scale. This would be a 
realistic engineering model whose test data would provide 
confidence that a full scale unit would perform as ex- 
pected, l.e., that a unit scaled up from the model would 
essentially perform like the model. If the original 
computer analysis were completely c. *rect, then the model, 
if directly scaled, would become full size design. 
However, since incorrect data were used in the reference 
design and model design, as previously discussed, and 
since some analytical modification would be expected, 
then a direct scaleup was not expected to truly provide 
a full size design. 

The subscale unit analytical design consisted of two parts. 
The first part was the exeunination of the theoretical 
criteria for scaling and the second part was the compu- 
ter analysis of model designs resulting in a choice whose 
output characteristics were in some direct relationship 
to those of the full scale unit. These tasks are described 
below. 


Scaling Procedure 

There are a number of methods for analyzing engineering 
processes and designing scale models. The best of these 
appears to be to write the differential equations describ- 
ing the process, non-dimensionalize the variables, extract 
coefficients which are dimensionless groups and set these 
groups equal to each other. 

Using the same assumptions in writing these equations as 
in writing the computer analysis, one can write the dif- 
ferential equations for infinitesimal elements in the 
height direction. An energy balance is made around the 
element of storage medium which includes axial conduction, 
surface losses, and heat transfer from (or to) the heat 
transfer fluid. An energy balance around the fluid then 
accounts for only transfer to (or from) the storage ele- 
ment . 

The differential energy balance around an element is writ- 
ten as 
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ds + (E ka) 


dx - UAj^(t^ - t)dx 


It 


ill 

3s* 


i!l 

3x* 


2xk. (t-t„) 
1 s 



while the energy balance around the fluid is written as 


3t- at. 

cpV^ Fr~ 3x~ UAj^(t - t^)dx. 


Note: All symbols are defined in the Scaling Procedure 

Symbol List presented as a fold-out in Appendix E. 


Note that the conduction down the tubes is taken along the 
tube locus since it is coiled. The derivatives along ds, 
a tube incremental length, may be transferred to the dx 
direction by 


ds 



dx 


so that 



38 



Moreover, making the substitutions 




2 

P 


V * 



0 



results in the energy balances being 
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The terms which now want to remain constant, i.e., the 
coefficients which are dimensionless groups are 


^t^t^ (Ska)P 
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and the 9(y,v) and 6^(p,v) will behave exactly the same in 

the model and in the full size unit. In fact, due to the 
temperature dependence of many of the properties, it is 
desirable to have t(p,v) and tj,(u»v) behave the same and 

then the properties at any point p,v will be identical. 
Were it not for this dependence, the experiment could be 
carried out over a different temperature range by select- 
ing a different tj^. 

While it is possible to manipulate the coefficients in a 
variety of ways and probably arrive at a reduced scale 
version which in some way represents the real process, 
consider the most important term in the process, the heat 
transfer between the T-66 and the Thermkeep. This is 
described by the overall heat transfer coefficient, U, 
which is composed of the fluld-to-wall resistance, the 
tube wall itself, the solid build up on the tube, and 
the solid to liquid phase resistance. 

This term is quite complex and it may be concluded that 
the best, if not the only, way to assure similitude would 
be to retain the cycle period, P, the tube configuration, 
the amount of Thermkeep per tube, and the flow per tube. 
Then the tubes in the model act essentially identically 
to the tubes in the full scale unit. This implies that 
the model should be a full height version of cross-sec- 
tional area reduced by the scaling factor as well as the 
number of tubes, total flow, and Implicitly the mass of 
storage medium. The duty cycle is, of course, the same 
with respect to time but with heat flows reduced by the 
scaling factor. 
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Consider what a scaling like this does to each of the 
dimensionless groups defined above. 


( 1 ) 


w 


This term has the same values of k, , P, c , and L, for 

t s t 

the model and the full scale unit. The total tube con- 
duction area Is reduced by the scaling factor since the 
total number of tubes Is reduced accordingly and the total 
mass of storage medium is reduced also by the scaling 
factor. Thus, this term Is unchanged. 


( 2 ) 


(E ka)P 


c M h 
s 


V 


The values of the thermal conductivities, k, P, c , and h 

* s V 

are identical for the model and the full size unit. If 
the thickness of the vessel is reduced by the square root 
of the scaling factor and since the circumference of the 
vessel drops by the square root of the scaling factor, 
then a terra due to the vessel wall goes down by this fac- 
tor. The medium axial conduction area also automatically 
drops by the scaling factor. In the denominator, M drops 
by the scaling factor so that the group remains unchanged 
from model to full scale. 


The values of U, P, and c are identical in both the model 

s 

and the full size unit, che rationale for equal U having 
been already explained. The tube surface area is scaled 
down because of the reduced number of tubes and the Therm- 
keep total mass is also scaled by the same factor so 
that the ratio of A/M is unchanged, and the group remains 
unchanged . 
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P and P are the same while the total flow m Is scaled 
and the total volume of T-66 In the tubes Is scaled so 
that the group remains unchanged. 


(5) 


UAP 

pcV 


The term U, P, p, and c are Identical for both the model 
and the full scale unit while the tube surface area and 
the total volume of T-66 In the tubes are scaled Identi- 
cally. Thus the group remains unchanged. 


( 6 ) 


2-nk. h P 
1 V 



The values of k. » h , P, and c are identical for the model 
1 V* * s 

and the full scale unit. Retaining the equality of the 
groups requires setting 



Noting that the storage medium masses are in the ratio 
of the scaling factor, R, this requires that 



model 


>\S. 






or 



Referring to the reference design discussions 



= 3-^7 
2.23 


1.55 


which is to say that 



(1.56) 


10 


85. 


The model shroud would be 60 m (197') in diameter! This 
is caused by the cylindrical conduction geometry and the 
fact that the vessel surface decreases only as the square 
root of the scaling factor. 

The solution to this problem was to reduce the surface 
losses in a different way. Rather than Increase the 
insulation thickness, we can use a reasonable insulation 
thickness and reduce the temperature difference between 
the model vessel and the model shroud. This Involves the 
heating of the shroud to an adequately high temperature. 
Consequently, for safety reasons and to conserve the 
amount of shroud heating required to maintain a desired 
temperature, another layer of insulation and an outer 
unheated shroud are used. 



The use of the dlfferentlui equations to scale the full 
size unit neglects the heat transfer at the top and bot- 
tom ends of the unit. However, It should be clear that 
since the cross-sectional area of the unit has been re- 
duced by the scaling factor, and since the number of 
tubes has also been reduced accordingly, the end losses 
approximately go down by the scaling factor as well. 


Subscale Unit Design 


As was discussed In the previous section, the scaling 
guidelines are straightforward and the model was designed 
as a full height cylindrical section of the full scale 
unit. The cross-sectional area is one-tenth of that of 

-I 

the full scale so that the diameter is or .316 

VTo 

of that of the full scale unit. The number of tubes was 
reduced from 252 to 25 but the length was identical. 

The major analytical design task with respect to the model 
was that of selecting a temperature to which to bring the 
heated Inner shroud in order to sufficiently reduce the 
surface losses from the vessel. The result of this was 
that the model would contain ^5. 7 cm (18”) of mineral 
wool insulation between the vessel and the heated shroud 
and the shroud temperature would be maintained at ^98 K 
(225 C) . The unit would then have unotner 15.2 cm (6") 
of insulation contained by an outer shroud. 

The criterion was used that the surface losses oe the 
same fraction of the total capacity in both the model and 
the full size unit. While this criterion was met, there 
is a difference between the model and the full size unit. 
Since the effective ambient temperature has been increased 
and since there exists a gradient from the bottom to the 
top of the vessel, the distribution of temperature dif- 
ference between the vessel and ambient is distorted. 

In the full scale unit, with an ambient temperature of, 
say, 296 K (23 C), a top end temperature of 58^ K (311 C) 
and a bottom end temperature of 516 K (2^3 C), the ratio 
of AT*s comparing the top to the bottom is (584-296)/ 
(516-296) = 1.31. In the model, with an effective ambient 
temperature of 498 K (225 C) and the same top and bottom 
end temperature, the ratio of AT's is ( 584-498)/( 516-498) 
4 . 78 . Consequently, while the total surface loss is In 
the right proportion, the top end of the unit is penalised 



somewhat In that It loses relatively more heat than it 
does In the full size unit. Since overall the surface 
losses are only a couple of percent of the total capacity 
per day, the redistribution should be unimportant. 

P5gure 17 shows the storage response to solar cycle dy- 
namics for the 1/10 scale model. This may be compared 
with Figure 15 (page 69), the response for the full scale 
unit. It Is substantially identical. Figure l8 shows 
the fluid outlet temperature at the hot end (during dis- 
charge) vs. cycle time for the model. This may be com- 
pared with Figure l6 (page 71) for the full-scale unit. 
Again, the two are substant'* ally the same. Note, however, 
that these results were obtained with the incorrect en- 
thalpy as was the reference design. They are shown to 
demonstrate similarity between the two designs. 

The model is basically an insulated vertical cylinder. 

The storage vessel is 0.70 m (27.5") diamater by 2.23 ni 
(87.8") tall. Since it has 0.152 m (6") of insulation 
between the heated shroud and the outside shroud, the 
overall size is 1.93 m (6'^") diameter by 3-^5 m (11 '4") 
tall. This is somewhat large for the Comstock & Wescott 
building so it was installed in an area where a three 
foot deep pit had been previously constructed for a simi- 
lar project. 

Provision was made to remove the heat exchanger (should 
a T-66 leak develop). The cold side inlet comes into 
the unit from the top (or hot side) via an insulated 
duct so that heat is not transferred to the T-66 passing 
through it (due to contact with hot end Thermkeep). 

Heating elements were installed to keep the Thermkeep 
liquid or bring it back to the liquid state to facili- 
tate replacement of the heat exchanger after repair. 

The vessel contains 1500 kg of Thermkeep and 25 tubes, 
0.635 cm (1/4") o.d. by 23.1 m ('75’9'*) long, through 
which the T-66 passes. Measurements of temperature were 
made on the vessel surface and via horizontal wells ex- 
tending into the vessel. These gave a measure of the 
Thermkeep thermal gradient. Their meaning and position- 
ing will be discussed '.iider data analysis. 



Construction of Subscale TES Unit 


The storage unit Is comprised of: 

... a steel vessel 

... a he?v exchanger which does not pene- 
trate the bottom and Is hung from the 
lid of the vessel 

... a number of thermocouples set Into 
wells which penetrate the vessel 

... a number of thermocouples bonded to 
the surface of the vessel 

. . . strip heaters attached to the outside 
surface of the vessel 

... ^5*7 cm (l8”) of mineral wool Insula- 

tion 

... a light gage sheet metal Inner shroud 
to which are fastened flexible Insu- 
lated heating wires 

... 15*2 cm (6") of mineral wool Insula- 

tion 

... an outer shroud 

. . . support legs 

. . . support struts 

The vessel was purchased from an outside shop. It stands 
2.i|i| m (96") tall, has an o.d. of .718 m (28"); the side 
wall is .63 cm (1/4") thick, and the base plate is 1.27 cm 
(1/2") thick. Twenty holes are drilled into the vessel 
for insertion of the thermocouples, figure 19 is a photo- 
graph of the vessel lying on its side with the thermo- 
couple wells pa'-tially inserted. 

Figure 20 shows how the thermocouple wells are formed. 

Type J (iron-constantan) thermocouple wire was stripped 
near thp end and Inserted into the steel well. The end of 
the well was pinched to clamp the ends of the wires and 
the flattened end was welded to seal it. The wells were 
inserted from Inside the vessel to the proper depth and 
welded from the outside to the vessel wall. 
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Figure l8. 1/10 model design — Thermkeep Thermal Energy Storage Unit. 
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This design was selected to Insure aaxlmum contact between 
the thermocouples and the surface to be measured. If 
thermocouples were merelj inserted, there was a possibility 
of poor themal contact. On the other hand, removal Is 

difficult and would require replacement by simply insert- 
ing thermocouples If a failure were to occur, but, due 
to the double shroud construction and the large quantity 
of insulation, the wells are not easily accessible anyway . 

There are 20 thermocouple wells, arranged in two rows of 
10 each, 15^® apart from each other and at different 
depths plus 52 on thj surface of the vessel at other posi- 
tions . The heat exchanger was designed to allow the tube 
coils to clear the rows of thermocouples should removal 
of the heat exchanger be necessary. As It turned out, 
at one point , the heat exchanger had to be removed due 
to a leak which was discovered before adding the Therm- 
keep and a couple of wells did not clear the coils and 
had to be broken to allow removal. These were replaced 
and the heat exchanger was additionally constrained to 
prevent this from occurring again. 

The heat exchanger consists of twenty-five 23.1 m (76 ' ) 
long steel tubes, .635 cm (1/M”) outside diameter by 
.7 mm ( .03") wall thickness. They were coiled on a 10 cm 
(M”) i.d. with a 3 cm (1.2”) pitch, for total coil height 
of 2.2 m (87”). They were arranged Inside the cylindrical 
vessel to be as evenly distributed as possible allowing 
for thermocouple well and lower manifold inlet tube clear- 
ance . The coils were provided with wires between turns 
to prevent sagging due to solid weight so that the pitch 
would be maintained. Wire lashing was provided between 
tubes to hold the lateral spacing between them since they 
are springlike otherwise and due to flexing would have 
uncontrolled separation. 

Upper and lower manifolds were fabricated from two steel 
tubes each of 1,7 cm (.6?”) outside diameter and 2,3 mm 
( .09”) wall thickness rolled Into circles of 58,3 cm { 23" ) 
and 35.8 cm (IM”), then cross-connected by a straight 
length of 2.7 cm (1 .05”) o.d. steel pipe . The lower mani- 
fold was split In two spots to allow for clearing the 
thermocouple wells for removal . Tubes were welded both 
to the circular and the straight parts of the manifolds. 

The two manifolds were spaced In the height direction by 
a steel bar welded to the top manifold and threaded to a 
fixture on the bottom manifold. 



Figure 21 Is a photograph of the two manifolds connected 
by the bar. Figure 22 is a photograph of the nearly 
completed heat exchanger being welded. Figure 23 Is a 
photograph of the heat exchanger looking from the bottom 
end which shows the split lower manifold and the clear- 
ance spaces. 

The heat exchanger was Inserted into the vessel from the 
open (top) end, three struts were welded across the top 
of the vessel and the heat exchanger was hung by hooks 
from these struts . The thermocouple wells were inserted 
to their approximate depths and welded to the vessel wall. 
Figure 24 is a photograph of the assembled heat exchanger 
inside the vessel looking Into the top. 

The bottom ends of the two shrouds were assembled onto 
the supporting feet for the vessel in the pit. The ves- 
sel was then set on top of its feet. Marinite Insulating 
blocks were used to support the vessel and prevent a 
direct steel conduction path from the bottom of the ves- 
sel to the outside. Figure 25 is a photograph of the 
vessel standing on the above-described assembly showing 
the heating cable already fastened to the Inner shroud. 

The thermocouples were then bonded to the vessel surface 
and nine strip heaters of .75 kilowatts each were also 
attached to the vessel surface. Insulation was packed 
into the bottom of the shrouds and the sidewalls of the 
heated shroud were attached save for a section allowing 
space to walk Inside to continue the insulatin’* task. 
Figure 26 shows the completed vessel with part of the 
sidewall of the heated shroud in place . The inner insu- 
lation and the completed heated shroud was assembled . 
Figure 27 shows this completed assembly sitting In the 
pit . 

The connections for the transfer of T-66 from the test 
rig to the storage unit were then made and the outer 
shroud and its Insulation were then completed. Figure 28 
shows a portion of the completed storage unit. Once the 
test rig was completed and checked out and it was possible 
to flow hot T-66 through the heat exchanger, the Thermkeep 
was added to the vessel. Thermkeep comes In the form of 
solid flakes and filling was done in steps with partial 
loads melted down to fill in the void spaces before add- 
ing more. Finally, the top of the . unit was assembled 
and insulated. 
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SUBSCALE TEST SYSTEM 


A facility has been built which was used to test the 
performance of the model storage unit. This facility 
is required to deliver hot T-66 to the unit during charg- 
ing with provision made to reheat the outlet T-66 as 
would be done by the solar collectors in the full-scale 
installation. On the other hand, it is required to de- 
liver cool T-66 to the unit during discharging with pro- 
vision made to recool the outlet T-66 as would be done 
by the toluene boiler in the full scale installation. 

The test modes, including the solar daily cycle simula- 
tion, were to be characterized by constant inlet temp- 
erature of the T-66 both during charging and discharg- 
ing, but, mainly in solar cycle testing, by variable 
flow rate. The operation of the collector field in the 
Sandia Laboratories' system maintains a fixed collector 
outlet temperature of 584 K (311 C) by variation of T-66 
flow rate or dumping of excess heat if flow control can- 
not keep the temperature down due either to intense heat- 
ing or lack of available storage capacity. 

Consequently, T-66 is supplied to the model storage unit 
from a high capacity, well mixed reservoir controlled to 
the appropriate delivery temperature. Figure 29 is a 
schematic diagram of the test apparatus. It does not 
show incidental details, such as safety re"'-ief valves, 
over-temperature safety shutoffs and .low liquid level 
safety shutoffs, pressure gauges, arid temperature sensors 
installed in the actual system. 

The reservoir is a 66 gallon horizontal cylindrical steel 
vessel. It is contained in a rectangular sheet steel box 
insulated by mineral wool. Set point temperature is main- 
tained by three 10-kilowatt, 480 V tubular electric immer- 
sion heaters. This allows supply of the rated 25 kW of 
heating to the T-66 plus a margin to compensate for losses. 

Control for the heaters to maintain set point is provided 
by a proportional plus reset plus rate type of driver 
and a solid state zero-crossover-firing power controller 
rated at 4l kVA. The rate function provides fast response 
while the reset feature provides long term accuracy on 
set point. A high temoeratui'e water cooled centrifugal 
circulatin'" pump is used to mix the T-66 in the reservoir 
to prevent thermal stratification and the consequent poor 
temperature control which was experienced before installa- 
tion of the pump. 
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When the mode of operation of the storage unit is that 
of charging, the T-66 exits from the bottom manifold 
cooler than it entered. Then the heaters are all that 
are required to maintain delivery temperature. On the 
other hand, during the discharging mode, the T-66 exits 
from the top manifold hotter than delivery temperature. 

Here the cooling whici. the toluene boiler would cause 
must be provided. 

This cooling is done .ii a heat exchanger immersed in 
water in a 55-gallon drum. The heat exchanger is formed 
by two helical colls of 2.5^ cm (1") copper tubing mani- 
folded in parallel. Facility water flows into the drum 
at the bottom and exits hotter from a drain on the side 
of the drum near the top. The heat exchanger and water 
flow rate were chosen to inhibit boiling of the water 
(assuming a 584 K (311 C) T-66 inlet temperature), so 
that energy balances could be made from measured water 
temperature and flow rates. 

The performance of the heat exchanger was calculated so 
that, without having to adjust the flow of water, the 
T-66 would be at least cooled below set-point delivery 
temperature. Moreover, it should never overcool so much 
that the 30 kW of heating capacity would be insufficient 
to return the T-66 temperature to set point. Therefore, 
no cooling control is required, the cooler is completely 
passive, and the heating control, which is quite depend- 
able, is used both for charging and discharging tests to 
hold a constant delivery temperature. 

The expansion tank, which compensates for the change in 
density of the T-66 with temperature changes, was fash- 
ioned after the guidelines described by Monsanto Company 
brochure, "Therminol Heat Transfer Fluids, a Design, 
Operating, and Maintenance Guide for Low-Cost, Low-Pres- 
sure Heat Transfer Systems.” The capacity of the tank is 
such that the level of the T-66 stays within the tank 
regardless of the temperature. It has a side-irounted 
sight glass to observe the T-66 level. A slight pressure 
of dry nitrogen gas is maintained above the T-66 to pre- 
vent oxidation if air were present. A low level safety 
switch which shuts off the reservoir heaters was installed. 

The T-66 is pumped to the storage unit by a high tempera- 
ture rotary pump with graphite asbestos packing. Maximum 
flow is 7.5 gpm. The pump is driven by a shunt wound d.c. 
motor controlled by a motor controller which can be oper- 
ated either manually or by command from an analog control 
system. The latter option was required to allow scheduling 
of the T-66 flow rate for solar eye e simulation. 
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Figure 29. Test cli'cult schematic 







A high temperature turbine flowmeter is used in conjunc- 
tion with a flow rate monitor (D/A converter) to generate 
a flow signal. The flow signal is used as a measure of 
the T-66 flow rate and also as a feedback signal for an 
analog flow control circuit commanding the motor control- 
ler for the delivery pump. 

Two groups of valves can be seen on the schematic diagram 
(Figure 29). One group (V1-V4) is for regulation of 
function of the expansion tank. During startup, VI and 
V2 are open, V3 and V4 are closed. This allows T-66 to 
pass through the expansion tank for deaeration. After 
a short period of time, V3 is opened, VI and V2 are 
closed, and V4 which is smaller and on a smaller diameter 
pipe is open. This valve and pipe are sized smaller so 
that the expansion is allowed but thermal syphoning of 
the T-66 is inhibited. 

The second group of four valves (V5-V8) controls the flow 
direction to the thermal storage unit . During a charging 
test (heating the storage medium), T-66 flows from the 
reservoir through valve V? (V8 is closed), into the top 
manifold of the unit, through the heat exchanger, out of 
the bottom manifold, through valve V6 (V5 is closed) and 
back to ‘■he reservoir for reheating. During a discharge 
test (cooling of the storage medium and heating of the 
T-66), T-66 flows from the reservoir through valve V5 
(V6 is closed), into the botton manifold, through the heat 
exchanger, out of the top manifold, through valve V8 (V7 
is closed) , through a heat exchanger which cools the T-66 
with facility water to below the reservoir control tempera- 
ture and back to the reservoir for reheating. 

Figure 30 is a photograph of all of the aforementioned 
components of the test facility except for the expansion 
tank which is elevated beyond the upper edge of the photo- 
graph. The drum, containing the heat exchanger for cool- 
ing the T-66, is in the left foreground with the drain to 
the left side. The rectangular box to the rear contains 
the T-66 heated reservoir. Connections to the upper and 
lower manifolds of the storage unit are in the right 
foreground with the edge of the pit containing the unit 
also visible. This photograph was taken before the pipes 
were Insulated with Fiberglas. 

Pressure measurements using standard Bourdon-tube gauges 
were made at various points in the circuit: upstream of 

the high temperature rotary pump; between the pump and 
filter; between the filter and the flowmeter; upper heat 
exchanger manifold; and lower heat exchanger manifold. 
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In addition to the temperature measurements inside the 
storage unit, certain temperatures were measured in the 
external circuit: water supply temperature, water drain 

temperature; T-66 temperatures into and out of the cooler 
heat exchanger; upper and lower manifolds of the storage 
heat exchanger; and two points inside the heated reser- 
voir plus the control point which is near the delivery 
tube. All temperature sensors are type J (iron-constantan) 
thermocouples . 

Figure 31 is a block diagram schematic of the circuit 
logic for the analog flow control circuit. The feed- 
back loop compares the measured flow rate from the out- 
put of the flowmeter with a desired flow rate and causes 
the motor controller to adjust pump speed until the error 
is effectively zero. For steady flow rate testing, a 
fixed value of desired flow is input to the comparator. 

The rest of the circuitry is used to program the f}ow 
using the cam program controller. 

The block function q(t) represents the cam program con- 
troller. To use it, one takes a circular disk, which 
rotates at a predetermined rate in the controller, and 
forms the desired function by cutting the disk accord- 
ingly. A cam follower displaces an amount determined by 
the disk pattern and generates a 0 to 8 VDC analog signal. 
Normally this is coupled to an appropriate temperature 
controller which then could, for example, regulate a 
furnace temperature to achieve a certain controlled pro- 
cess. 

In this case the controller output is taken to be a rate 
of heat flow, either the excess heat from the solar collec- 
tors which one is attempting to put into storage, or the 
heat deficiency of the solar collectors which one is 
attempting to extract from storage. As was discussed in 
the section on analysis, tne flow rate to the storage 
unit is determined differently for charging than for dis- 
charging. 

The storage unit is, in essence, a large thermal regenera- 
tor. It is characteristic of such a device that, as heat- 
ing or cooling proceeds, the outlet temperature of the 
fluid changes. During discharge, the delivery tempera- 
ture to the boiler will drop and during charging, the 
delivery temperature to the solar collectors will rise. 

In the former case, the temperature deficiency is accepted 
and it is assumed that, possibly, e -r.eans of boosting the 
temperature, e.g., a fuel-fired burner would be provided. 
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Figure 31. Main T-66 pump control circuit schematic. 
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In the latter case, it is assumed that the flow rate 
through the collectors would be Increased to attempt to 
absorb all the heat collected in the face of a diminish- 
ing temperature rise. The additional circuitry shown 
accomplishes this additional flow rate adjustment. 

For discharging, a cam is cut which outputs the differ- 
ence between the fixed flow rate demanded by the boiler 
and the flow rate which is expected from the solar col- 
lectors as described by the characteristic curve. The 
multiplier (M2) shown on the block diagram would be by- 
passed and the flow demanded would be in direct propor- 
tion to the output of the program controller. The shape 
of the flow demand is such that flow begins at zero in 
the late afternoon when the collector output becomes 
equal to the boiler demand. Thereafter, the storage flow 
increases until the collector output ceases and all the 
boiler flow passes through storage. The flow rate is 
constant until the beginning of the nighttime idle period 
at which point it becomes zero. In the early morning, 
after the idle, storage flow begins at full boiler flow, 
stays constant until collector output begins and decreases 
steadily to zero flow which occurs when the collector out- 
put rises to the value of the boiler demand. 

During charging, one would like to absorb all the output 
of the solar collectors. If this occurred, one could 
write the equation 
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where T^^ is the fixed delivery temperature from the collec- 
tors (584 Kj 311 C), Tj^ is the variable outlet temperature 
from storage, c^ is the average T-66 specific heat over the 
range T to T„, w« is the required flow and q is the 
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heat available to go into storage. The flow rate function 
becomes 


w^ 




Of. (T„ - Tf_) 
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The control circuit takes the value of T^, as measured 

at the lower manifold, and uses It first to form the dif- 
ference, - Tj^, In the summer shown, and second to form 

the average specific heat In the function block c^ (T^) 

where this can be described as a function of only with 

fixed. The outputs of both of these operations are 

Input to the multiplier to form the product - Tj^) 

which Is then Inverted to form l/cj,(Tj^ - T^^). This value 

together with q(t) Is Input Into the multiplier to form 
the demanded flow according to the aforementioned equa- 
tion. Lastly, it is expected that the pump would be 
limited to some realistic flow rate, at least to limit 
the pressure loss through the system so that the func- 
tion w(w^) is essentially a limiter whereby w = w^ up to 

w at which point It no longer increases. Moreover, 
max ^ ® ’ 

the flowmeter cannot function below a minimum of about 
0.75 gpm so that the solar characteristic cams were modi- 
fied to demand no less than this. 

Debugging the facility involved the usual problems asso- 
ciated with such a system. Erratic behaviour of the 
turbine flowmeter was traced to two problems. First, 
the wrong sensor head was provided. Rather, one rated 
for only 355 K (l80 P) was received and had to be replaced 
by the factory. Second, problems were traced to the D/A 
converter and certain factory recommended adjustments 
had to be done. 

It was expected that the test flow rate would be adequate 
to agitate the T-66 in the reservoir. However, the T-66 
formed a thermocllne (thermal stratification) and a cen- 
trifugal circulating pump had to be ordered. T-66 leaks 
developed in the plumbing, a problem which is known to 
be quite prevalent with this fluid, and these had to be 
repaired. Other minor problems were also rectified. 

Overall, it is felt that a very useful test facility now 
exists. It can be used to model any test duty cycle which 
requires fixed temperatures and variable flow rates. The 
next section discusses the performcnce testing which was 
done on the model storage unit using this facility. 
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SUBSCALE TEST RESULTS AND EVALUATION 


^ .Preliminary Discussion 


The steps which led to the final form of the analysis 
will be briefly discussed. Then all the data and the 
computer predictions will be presented and compared, the 
computer predictions being those from the final form of 
the analysis. 

A number of changes were made to go from the design 
analysis used In the reference design tc the test analy- 
sis. Some relate to simulation of test conditions and 
will not affect the analysis to be used for design of a 
full scale unit. These modifications involve the ability 
to input the proper starting conditions of the unit (i.e., 
temperature distribution) , to input a variable T-66 in- 
let temperature for runs where this temperature had not 
yet stabilized by the start of the test, to input the 
proper constant flow rate, and to override the solar 
cycle schedule normally in the program. Three decks of 
computer cards were produced from the original program. 

One simulates a single-pass heatup, one simulates a sin- 
gle-pass drawdown, and one simulates the square-wave 
cycling. 

Other changes relate to updated knowledge of material 
properties. Since the reference design was produced, new 
Information was received from the Monsanto Company on the 
properties of T-66 and new information was provided by 
Dynatech R/D Company on the properties of Thermkeep. 

These changes were Incorporated into the program. 

New tables of the enthalpy and solid fraction of Therm- 
keep as a function of temperature were input. These were 
composed on a rational basis but are not based upon labo- 
ratory measurements. As mentioned, the h(T) data ex- 
pected from the subcontractor were never properly pro- 
duced and C&W used the best information previously avail- 
able . 

The enthalpy vs . temperature curve was adapted from a 
temperature vs. time cooling curve done in the C&W labo- 
ratory (see Figure 32). The use of this has produced 
reasonably good results, as will be seen, but it is be- 
lieved that some of the lack of agreement between test 
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and analysis can be attributed to the lack of accuracy 
In the enthalpy-temperature correlation. The solid 
fraction, (|i(T), was adapted from a phase diagram for the 
binary system NaOH-NaNO^ (Ref. 3). This also is not con- 
sidered to be highly accurate, but It is not considered as 
significant a factor as the Inaccuracy in the h(T) curve 
In contributing areas of lack of agreement. 

The only Important change made In the calculation proce- 
dure was the inclusion of a correction to the fluid side 
heat transfer coefficient caused by the winding of the 
tubes in a springlike helix (Ref. li ) . The original com- 
putations were based upon straight tube correlations. 

The new correlation says that the critical Reynolds num- 
ber for transition from laminar to turbulent flow is 
Increased according to 
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where d. is the tube diameter and d. is the helix diameter. 
z n 

Moreover, the heat transfer coefficient in laminar flow is 
increased as follows. Define 


and 





Where Pr is the fluid Prandtl number, then 


Nu 
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Figure 32. Thermkeep enthalpy from C&W cooldown data. 





where Nu Is the Nusselt number for curved tubes and Nu„ 
e s 

Is the Nusselt number for straight tubes. Using this 
correction for the low flow rate tests brought the agree 
ment from poor to good and this verifies its usefulness. 
Other than this change, the method of the computer analy 
sis stands essentially as it was for the reference de- 
sign analysis. 


Description of Test Procedures 


The model TES unit was Installed into the test facility 
and then subjected to a variety of tests to determine 
its ability to accept heat from and deliver heat to the 
heat transfer fluid. These tests fell into four basic 
categories, each with its own objective. 


1. One-Pass Charge/Discharge Testing 


In this type of test, which was the basic sour-ce of 
data for tuning the computer analysis, the unit is 
bT*ought to either a typical charged or a typical dis 
charged condition and then either discharged or 
c..arged at a constant flow rate. The change of the 
internal condition of the unit was observed as a 
function of time as well as the change in delivery 
temperature of the T-66. A number of different flow 
rates were checked ranging f*rom O.CZ kg/sec (1 gpm) 
to 0.29 kg/sec ( 6.3 gpm). 

2. Idle Tests 


In this type of test, the unit was brought to a hot 
or charged condition, flow was stopr^d, ard the 
change in internal condition was recordel as a func- 
tion of time. The purpose of this test Is to elimi- 
nate the effect of heat transfer with the fluid and 
observe only the effects of axial conduction and 
losses to the environment. 



3. Square-Wave Cyclic Testing 


In this type of test, the unit is brought to a 
typical charged condition and then alternately 
discharged and charged at constant T-66 flow rates 
with minimal changeover times. The purpose of 
this test is to verify that the unit approaches 
a stable condition under cyclic operation and that, 
without specifying the starting condition for each 
segment of the test, except for the initial start- 
ing condition, the analysis could predict the 
performance after a number of repetitive cycles. 


4. Solar Dally Cycle Testing 


These tests were run to describe how the unit 
would function under solar cycle conditions where 
flows would vary according to a prescribed sched- 
ule and a control mode. This continual variation 
would create difficulty in using this type of 
test for tuning the analysis but would be a proof 
test of the unit under realistic conditions. 


The test conditions were selected to be 1/10 scale of the 
typical full scale conditions. T-66 temperature during 
a charge was held at around 584 K (311 C) while during a 
discharge, it was held at around 516 K (243 C); however, 
these changed somewhat from test to test. 

With a 25 kW heat load and a 68 K available overall fiT, 
a T-66 flow of about 0.14 kg/sec (3 gpm) is required. 

The test flows ranged from 0.05 kg/sec (1 gpm) to 0.29 kg/ 
sec (6.3 gpm) but most tests were run in the 0.09 to 0.l6 
kg/sec (2 to 3.5 gpm) range. A range of flows is desired 
since in an operational system, the flows would vary as 
the available collector output changes. 

The primary objective of testing was to determine if the 
computer analysis to be used for the design of large 
scale TES units is accurate in its description of how the 
unit would respond to a given stimulus. If reasonable 
accuracy can be obtained, then the analysis can he used 
for determination of efficiencies, delivery temperatures, 
pressure drops, and any other factors considered important 
in designing a full-scale storage unit . 



Comparison of Computed and Experimental System Response 


One-Pass Charge/Discharge Test Results 


For a steady flow discharge test, the unit Is brought to 
a reasonable charged condition and a flow of T-66 (at 
low cycle temperature) Is Introduced Into the bottom mani- 
fold. The hot-end outlet temperature of the T-66 as a 
function of time Is the critical characteristic In that 
temperature deficiency here Impacts directly on the tur- 
bine performance and must be compensated. Figures 33 
through ^0 plot the T-66 outlet temperature during draw- 
down tests for flow rates increasing from 0.05 kg/sec 
(1 gpm) to 0.29 kg/sec (6.3 gpm) . 

Except for test No. 007, which started at quite a low 
state of charge, the computer predictions agree reason- 
ably well with the measurements. The only observation 
that seems at least somewhat consistent is that between 
577 K and 566 K, the unit delivers more heat than pre- 
dicted and vice versa below 566 K. This suggests that 
the h(T) curve used should be reshaped to have le s of 
a plateau at 566 K with slightly more energy available 
above 566 K and a more rounded shape below 566 K. 

Figures 4l through 49 are comparisons of measured thermal 
distribution with computed distribution. The computer 
technique assigns one value of temperature to any depth 
and Ignores lateral gradients. As may be seen by some of 
the figures, especially in the lower temperature regions, 
the measured values do not necessarily fall on a single 
line. There are two axial x’ows of thermocouples, each 
row at a different radial position in the unit. Thus, 
the existence of radial (lateral) gradients would cause 
disagreement between each set. Radial gradients tend 
to occur in lower temperature regions because more of the 
Thermkeep is solid. Circulation in the liquid phase pre- 
vents this from occurring. 

Note that the C'^mputed axial temperature profiles do not 
extend to the top of the vessel. The initial height of 
the Th'rmkeep depends on the initial temperature of :;he 
Thermkeep, and actually varies in some cases over a 
range of roughly 20cm during operation of the system. 

The computed points are associated with the estimated 
total height of Thermkeep defined by the initial measured 
temperature profile. The estimated distance of the Therm- 
Keep surface from the top of the vessel is indicated on 
the figures. Thus, m.easured points above this estimated 
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Thermkeep surfaces are interpreted as being in the air 
space above the surface and are not computed by the 
analytical model. 

Again, except for test No. 007, the data and analysis 
agree quite well. The analysis predicts a longer dwell 
at 566 K than actually seems to occur, suggesting that 
the plateau on the h(T) curve may be too strongly de- 
fined. 

For a constant flow charging test, the unit is brought 

to a low state of charge and a flow of T-66 (at high 

cycle temperature) is Introduced into the top manifold. 

The cold-end outlet temperature of the T-66 as a func- 
tion of time is the critical chai’acteristlc in that 

rising outlet temperature requires additional flow to 
absorb the collector output without exceeding a speci- 
fied upper temperature limit. Figures 50 through 5^ 
plot the T-66 outlet temperature during charging tests 
for flow rates varying from 0.05 to 0.3 kg/sec (1 to 
3.5 gpm). 

Overall, the agreement is reasonably good. The compu- 
ter prediction of the outlet temperature is higher than 
measured outlet temperature. This is possibly due to 
over-estimating the resistance of the solid Thermkeep 
layer. The rate of rise of the temperature is consist- 
ently good which could be indicative of a fair agreement 
in h(T) in the low temperature regions. 

Figures 55 through 59 are comparisons of measured thermal 
distributions with computed distributions. Again, the 
computer seems to predict a stronger than measured pla- 
teau at 566 K which also may be due to incorrect shaping 
of h(T) in the vicinity of this teraperaLure. Otherwise, 
the agreement is reasonably good. The fact that both the 
cold-end outlet and cold-end distribution are off in the 
same way, i.e., that the computer predicts higher tempera- 
tures than the data, could indicate an error in h(T) 
rather than in the resistance of the solid layer. 
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Figure 35- T-66 outlet temperature discharge 

test 003 (? gpm) . 
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Figure 36 . T-66 outlet temperature discharge 

test 007 (2.5 gpm) . 
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Figure 37. T-66 outlet temperature discharge 

test 009-1 (2.8 gpm) . 



Figure 38. T-66 outlet temperature discharge 

test 001 (2.83 gpm) . 
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Figure 45 . Thermkeep temperature profile, dis- 
charge test 009-1 ( 2.79 gpm), at 2 hrs. 
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Figure 46. Thermkeep temperature profile, dis- 
charge test 001 ( 2.83 gpra) , at 2.33 hrs 
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Figure 47. Thermkeep temperature profile, dis- 
charge test 001 (2.83 gpm) at 5.33 hrs. 



Figure M 8 . Thermkeep temperature profile, dis- 
charge test 005 (3.5 gptn) , at 2 hrs. 
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Figure 57. Thermkeep temperature profile, 

charge test 6 (2.5 gpm) at 2 hrs. 



Figure 58 . Thermkeep temperature profile, 

charge test 002 (2.83 gpm) at 3 hrs. 
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Figure 59. Thermkeep temperatu 
test OQiJ ( 3.5 gpm) , 






Idle Tests 


For Test No. 010, the unit was charged and allowed to 
stand for 115 hours, during which time the thermal 
distribution was recorded. Figures 60 through 64 describe 
the computed and measured thermal distributions at the 
start of the test and after 25, 45, 65, and 115 hours. 

Two characteristics are evident. First, the real tem- 
perature drop Is greater than the computed. This may be 
due to the fact that the losses are actually higher than 
computed or due to the possible error In the h(T) curve 
which has already been suggested. In general, stand-by 
loss should be predictable although there are aspects of 
the model design which are not taken Into account In the 
analysis and which could Increase the stand-by losses. 

First of all, the vessel has an Internal height which Is 
155? greater than that occupied by the storage medium and 
causes more surface losses. Second, there are shroud 
support struts and support legs which cause heat loss but 
are not described by the program. Third, there is air 
circulation which can have a cooling effect on the ex- 
posed part of the heat exchanger in the open space at the 
top of the vessel and on the top surface of the Therm- 
keep. Moreover, values for insulation conductivity are 
not precise. There Is also a capacitive effect associa- 
ted with the insulation. When the vessel is heated, the 
Insulation, which has a low thermal dlffusivlty, lags 
behind In temperature rise. During a standby, the Insu- 
lation absorbs heat from the vessel In an attempt to come 
to equilibrium with It. 

Overall, there seems to be no single reason for the more 
rapid temperature drop In the data than In the computer 
predictions. Rather, a combination of effects all of 
which move the results in the same direction, seem to be 
Involved. 
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Square Wave Cyclic Testing 


Two tests were run with repetitive cycling and minimal 
changeover time, thereby minimizing the effect of erro- 
neous prediction of stand-by heat loss, as described 
above. During a flow test, the heat exchange between 
the fluid and the storage medium far exceeds the effect 
of stand-by losses and errors in stand-by loss are in- 
significant in terms of correlation. 

The two tests were run at flows of 0.1 kg/sec (2 gpm) 
for hours discharging, 4 hours charging, with one-hour 
periods for changeover. Test 017 ran two cycles, and 
test 0l8 ran three cycles. Pour cycles were desired each 
time but operational problems prevented completion of the 
full four cycles. 

Figures 65 and 66 describe the T-66 outlet temperatures 
during the first and second cycles of the two-cycle test 
(No. 017). The computed results are based on the ini- 
tial measured temperature profile at the beginning of 
the first cycles, l.e., only at the start of the test 
sequence . 

Figures 67 through 70 are the computed and measured 
thermal distribution at the end of each of the four 
segments. Agreement appears acceptable. 

Figures 71 through 73 describe the T-66 outlet tempera- 
tures during each of the three cycles of test No. 0l8. 
Figures 74 through 79 describe the measured and computed 
thermal distributions at the end of each of the six seg- 
ments. Again, agreement is acceptable. 
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Figure 66. T-66 outlet temperature, cyclic 

test 017 (2 gpm) , cycle 2. 



Figure 67 . Thermkeep temperature profile, cyclic 
test 017 (2 gpm), at end of first 
discharge phase. 
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Figure 72. 


T-66 outlet temperature, cyclic 
test 018 (2 gpra), cycle 2. 
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Figure 73. T-66 outlet temperature, cyclic 

test 0l8 (2 gpm), cycle 3. 
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Therrakeep temperature profile, cyclic 
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Thermkeep temperature profile, cyclic 
test 018 (2 gpm), at end of second 
charge phase. 







Figure 78. Thermkeop temperature profile, cyclic 
test 018 (2 gpm), at end of third 
discharge phase. 
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Figure 79. Thermkeep temperature profile, cyclic 
test 018 (2 gpm), at end of third 
charge phase. 







Solar Cycle Simulation Tests 


Two solar cycle tests were run with data being measured 
every 3C minutes. The first (test 020) was run as a 
system shakedown with the charging and discharging sched- 
ules compressed to 6/10 of normal, so that the twelve 
hours of boiler demand were decreased to 7.2 hours. This 
allowed the test to be run in 'he course of two normal 
workdays, and as there were no system malfunctions, the 
data are being presented. The second test (test 023) 
was run with the insolation and boiler demand just as 
described in Figure 1 (page 13). 

The order of events in the computer program was modified 
so that a starting profile was input and computations 
began at the predawn start of boiler demand (time T3), 
and a cycle runs from then through the end of the over- 
night idle period. Next, the program used an input 
table of actual flow rates rather than calculating flow 
on the basis of the solar characteristic and the inlet 
and outlet temperatures. This modification removed any 
discrepancy between the actual flow as determined by 
the test bed electronics and the computed flow. Lastly, 
a change was made to allow for the fact that the test 
bed has a minimum flow rate of .75 gpm (about .04 kg/sec). 
Below this rate, the flowmeter is not accurate and the 
load on the pump causes it to pulsate. Therefore, the 
flow rate will stay at the minimum rate in response to 
any signal to pump at .75 gpm or lower. The computer 
analysis allowed an idle period between charging and 
discharging phases, because during the tests the flow 
was cut off early at the end of a phase or started late 
at the beginning of a phase. The times of cutoff or 
startup were calculated so that the total desired amount 
of heat was still transferred. With these modifications, 
the computer analysis duplicates the actual events. 

The morning discharge period in which flow starts before 
dawn and later decreases as the sun rises is referred to 
in the data as the Phase 1 discharge to distinguish it 
from the afternoon, or Phase 2 discharge during which 
the sun sets and the unit continues to supply heat to 
the boiler unassisted. 

During each phase the T-66 inlet temperature was main- 
tained constant within limits — near 516 K during dis- 
charge and 534 K during charging. The actual inlet 
temperatures were input into the computer to eliminate 
variations due to external effects on the control system. 
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Figure 80 shows the computer correlation to T-66 outlet 
temperature for Test 020 and Figures 81-84 compare the 
computed and actual Thermkeep temperature profiles for 
the first dally cycle. Figure 85 shows the varying rate 
of T-66 flow during the charge cycle. Note that the 
highest rate attained was 0.31 kg/sec, whereas the maxi- 
mum possible rate Is about 0.38 kg/sec. This Is because 
the flow rate Is electronically adjusted to account for 
the Insolation characteristic and the return tempera- 
ture, just as the computer analysis did previously. In 
a relatively low state of charge, the return tempera- 
ture will be low and maintain a large At which requires 
less flow. 

Figures 86-90 show the computer correlations for the sec- 
ond cycle of Test 020. Agreement has improved for both 
charging and discharging phases and is very close at the 
end of the test (Figure 89). Figure 90 shows that the 
maximum flow rate was momentarily reached, and although 
the return temperature was high, the flow rate dropped 
off as the collector output dropped. 

Solar Cycle Test 023 (Figures 91-101) was a full-scale 
simulation of two solar daily cycles. This test closely 
approximates the actual performance conditions of a full- 
scale storage unit in the field. 

Figure 91 shows the correlation between computed and 
actual fluid delivery temperatures for both cycles. 
Agreement is very good, with two exceptions: during both 

charge phases, when the T-66 temperature approaches 560 K, 
computed and actual data begin to diverge. Similarly, 
the Phase 1 discharge of the second cycle (from 24 to 
26.5 hours) Is slightly divergent in the 560 K range. 

The reason for this deviation has not been determined. 
Some factors which may be contributing are a misshapen 
enthalpy curve, some hysteresis which raises the latent 
heat plateau on heating and/or depresses the latent heat 
on cooling, possible segregation of the Thermkeep compo- 
nents which changes the specific heat and melting point 
of the bottom most elements, and expansion and circula- 
tion of hot liquid Thermkeep which shifts energies and 
volumes of the elements in ways unaccounted for by the 
computer. 

Although the predicted return temperature was lower than 
actual during charging, which might suggest that less 
energy was absorbed by the unit, the profiles at the end 
of the charge phases (Figures 94 and 99) show that 
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temperatures are slightly higher than predicted Inside 
the unit. In particular, Figure 99 shows an abrupt 
dip In the computed temperature at about 150 cm from 
the top of the vessel. This again Is a reflection of 
the strong plateau in the enthalpy curve, which experi- 
mental data suggests are not as strong as we assumed 
Below the level of the dip, the discrepancy is roughly 
15 K, or the same discrepancy which appears in the re- 
turn fluid temperature during charging (see Figure 91) 
at about 32 hours elapsed time. 

The final discharge phase showed good agreement between 
computed and actual data on the T-66 delivery tempera- 
ture, and the test ended with the temperature profile 
also agreeing closely (Figure 100). 

Figure 101 shows the variation in flow rate during the 
final charge phase. Although peak Insolation occurred 
at about 30 hours el8 ed time, the flow rate continued 
to rise and held near ..icxximura (.35 kg/sec) for about an 
hour, while the collector output continued to drop. 
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T-66 outlet temperature for 6/10 scale 
solar cycle simulation test. 
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Figure 81. Thermkeep starting temperature 
cycle 1. 






3ol»r Cycle Te»t OJO 

Themkeep 




Elapsed Time *1.6 hours 

»i ■— Surface 





o 

^0 — <« 0 0 



o 

Computed 





(KS9 SC020 8-15-70) 




\ 

/.o 0 




Measured"^ 

•C 0 




0 ® o 










N. ft 





X. ft 

U. wi 





JL r,L 


OIJ'TANCE PROM TOP Op STORAGE VESSEL, cA 


Figure 82. Thermkeep temperature at end of 
Phase 1 discharge, cycle 1. 
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Figure 83 . Thermk<»ep temperature at end 
of charge, cycle 1. 
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Figure 8^1. Thermkeep temperature at end of 
phase 2 discharge, cycle 1. 
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Figure 92. Thermkeep temperature. Test 023 
starting profile, cycle 1. 



Figure 93. Thermkeep temperature at end of 
phase 1 discharge, cycle 1. 
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Figure 96. T-66 flow rate variation during 

first charge cycle. 
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Figure 100. Thermkeep temperature at end of 
phase 2 discharge, cycle 2. 
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Figure 101. T-66 flow rate variation during 

second charge cycle. 
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PRELIMINARY DESIGN 


With the computer model now reliably predicting experi- 
mental results, it was used to produce an updated de- 
sign for a unit meeting the performance specifications 
as listed under Task 3, Page 7. 

The primary variables affecting system performance are 
the quantity of storage medium and the total heat ex- 
changer surface area. These two parameters were used 
to produce a first optimization. A figure of merit 
(P.O.M.) which is an efficiency based upon the deficit 
requirements of the unit (the amount of heat required 
to maintain the outlet temperature at a minimum of 580 K 
(307 C) during discharge) was used to describe the per- 
formance of each design. 


where 


P.O.M. 





* storage dally output 

Q = storage daily surface loss 

o 

Qjj * fluid deficit heating 


Surface areas from 55.7 m^ (600 ft*) to 167 m* (l800 ft*) 
were run in units ranging from 10,000 to 50,000 kg of 
Thermkeep. All were run with the optimum aspect ratio 
of 1 (vessel height to vessel diameter) and standard 
0.635 cm (1/^”) tubing of a given length and number to 
produce about 69 kPa (10 psi) of pressure drop in the 
heat exchanger with a flow of 5 kg/sec of Thermlnol-66 . 

Plgures 102 through 105 show the heat exchanger surface 
area and amounts of Thermkeep plotted against the figure 
of merit and the specific cost of storage in dollars per 
MJ. With some interpolation, the amount of Thermkeep 
was then cross plotted against specific cost at constant 
figures of merit (Fig. 106). 
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At this point, a first optimization is selected with one 
additional arbitrary limit — reasonable size. A vessel 
with an aspect ratio of 1 tends to exceed 12 feet in 
diameter at about 20,000 kg, above which the unit would 
probably have to be fabricated and shipped in more than 
one piec.e. The highest F.O.M. at lowest specific cost 
for a reasonably sized unit falls at an F.O.M. of .88, 
as specific cost of $11.20/MJ and uses 18,000 kg of Therm- 
keep. The heat exchanger has 111.5 (1200 ft^) of 

surface area. 

This design was obtained using 10 numerical elements for 
computer analysis; a finer grid of 25 elements Improves 
the computed performance- slightly and brings the cost 
down to $11.17/MJ. 

Next the design was tested to evaluate the effect of five 
minor parameters: tube size, fluid velocity, aspect ratio, 

insulation thickness, and peak flow through storage. 

Figure 107 shows the effect of increasing tube size both 
against efficiency and cost. 0.635 cm (1/4") is believed 
to be the smallest practical size to work with; and as tubes 
get larger, the cost increases while efficiency decreases. 
The 0.635 cm (1/4") tubes are therefore retained. 

Figure 108 shows the effect of Increasing fluid velocity 
by Increasing the length of the tubes and decreasing 
their number, retaining the same surface area. Perform- 
ance does increase with velocity as expected; however, if 
we assume a 1055 thermal-to-mechanical efficiency for driv- 
ing the pump, any Increase over 11 psl cannot be Justified. 
To go to 145 kPa (21 psi), for example, the pump would 
use an additional 54,000 kJ output. Therefore, the 76 kPa 
(11 psi) speicification was accepted. 

The effect of the vessel aspect ratio (height/diameter) 
on the specific cost and on the figure of merit is shown 
on Figure 109. Although the figure of merit is highest at 
the ratio 1.0, and the specific cost has a minimum at the 
ratio 2.0, the aspect ratio has so little effect on these 
parameters that it can be chosen on the basis cf other 
factors. The aspect ratio 1.0 is used in the Preliminary 
Design. 

Figure 110 shows the effect on the Figure of Merit and on 
specific cost of increasing Insulation thickness. Perform- 
ance does Increase but with diminishing returns. The 
Increase in total cost of the TES unit by Increasing the 
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insulation thickness from 0.457 m (18") to 0.61 m (24”) 
is $470 which affords about 1.5S5 increased performance. 

The increase from 0.457 m (18”) to 0.6l m (24") costs 
$560 and affords about 0 . 7 % increased performance. From 
0.61 m (24") to 0.76 m (30") costs $630 and affords less 
than 0.5?. Meanwhile, the loss from 0.457 m ( 18 ") of 
insulation is 3 •958 of the daily storage, and from 0 . 6 I m 
(24"), it is 3.3^. 0.76 m (30") makes the loss 2 . 9 %. 

Thus, the 0 . 6 I m (24") thickness was selected as the 
point at which returns begin to diminish most sharply. 

The effect of varying the peak allowable storage flow is 
shown in Figure 111. 5 kg/sec seems to be the point b ' ow 

which efficiency begins to deteriorate rapidly and cost 
begins to escalate markedly. Therefore, 5 kg/sec is re- 
tained in the design. 

The quantity of Thermkeep and the size of the heat ex- 
changer, and consequently the size and cost of the TES 
unit are strong functions of Q^j, the deficit heating. 

Figure 112 shows the relation between the specific cost 
of the TES unit and the deficit heating as percent of the 
dally output from storage when the unit is performing 
according to the solar cycle. The points on the curve 
correspond to a range of quantities of Thermkeep from 
12,000 kg to 27,000 kg. Throughout, the aspect ratio of 
the vessel is 1.0, and the total inside surface of the 
heat exchanger tubes is in the same proportion to the 
mass of Thermkeep as in the Preliminary Design, namely, 

0.00619 m*/kg. 


Preliminary Design 


Vessel height and diameter 
Insulation thickness 
System height and diameter 
System weight 

Amount of storage material 
Unit cost 

Heat exchanger 


Figure of Merit based upon 
a 25 -element grid 

Storage output 
Specific cost 
Insulation loss 
Deficit heating 


2.35 m (92.5") 

0.609 m (24") 

3.6 m (142") 

27290 kg (60,200 lb) 
18000 kg ( 39,700 lb) 

$33,700 

.00635 m (1/4") tubing; 
280 coils, each 27.73 
(91 ft) long 

.882 

3.01 X 10^ kJ 

$11.2/MJ 

3.2558 of storage output 
10.255 of storage output 
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Figure 102. Specific cost of storage for 

various heat exchanger surface 
areas at constant Thermkeep 
masses. 
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Figure of Merit 



Amount of Thermkeep, Mg 


Figure 103. Figures of merit for various Therm- 
keep masses at constant heat ex- 
changer surface areas. 
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Figure of Merit 







Figure 105. Specific cost of storage for 
various Therrakeep masses at 
constant heat exchanger surfaces. 
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Amount of Thermkeep, Mg 


Figure 106. Specific cost of storage for 
various Thermkeep menses at 
constant figures of merit. 
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Tube Outside Diameter, cm 


Figure of Merit 
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Figure 107. Specific cost of storage 
and figures of merit vs. 
tube diameters. 
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Figure of Merit Specific Cost, $/HJ 



Peak Pressure Drop, kPa 

Figure 108. Specific cost of storage, peak pump 
power, and figures of merit at vari- 
ous heat exchanger pressure drops. 
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Figure of Merit Specific Cost, $/MJ 



Vessel Aspect Ratio 


Figure 109. Specific cost of storage 
and figures of merit vs. 
aspect ratio. 
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Figure of Merit 



10.8 11.0 11.2 11.4 11.6 

Specific Cost, $/MJ 


Figure 110. Effect of various insula 
tion thicknesses on cost 
and performances. 



Specific Cost, $/MJ 



Pump Flow, kg/sec 


Figure 111. Specific cost of storage and 
figures of merit at various 
peak pump flows. 
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Deficit Heating, % of Dally Output 


_L 


Figure 112. Specific cost versus deficit heating. 
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Cost Estimates 


The estimated cost of the TES unit is calculated by the 
computer program by methods described in the preceding 
section titled Reference Design Study. Input data con- 
sists of the cost of material for each component, and 
a fabrication cost for each component which differs from 
part to part depending on the type of fabrication re- 
quired. The input data for cost calculations are given 
in Table I. The detailed result? »f the cost calcula- 
tion for the Reference Design a”c shown in Table II. 

The materials costs were those obtained at the time of 
the construction of the subscale unit, i.e., 1977. The 
fabrication costs are based upon the actual manufactur- 
ing costs of two Thermbank electric water heaters of the 
type described in the Background section of the Intro- 
duction. 

The estimated cost of the Preliminary Design TES unit 
is therefore representative of the cost to be expected 
if a very small number of units were to be fabricated. 
The fabrication component of the cost can be expected 
to decrease substantially if units are manufactured in 
large enough numbers so that the cost advantages of 
quantity production can be realized. 
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TABLE I 

COST CALCULATION INPUT DATA 


Shroud material cost, $/kg 0.509 

Shroud side fabrication cost, $/sq m 10.1 

Shroud side fabrication cost, $/kg 1.04 

Shroud top fabrication cost, $/sq ra 77.4 

Shroud top fabrication cost, $/kg 7.92 

Shroud bottom fabrication cost, $/sq m 99.8 

Shroud bottom fabrication cost, $/kg 7.67 

Shroud top thickness, m 0.00152 

Shroud side thickness, m 0.00122 

Shroud bottom thickness, m 0.00340 

Vessel wall material cost, $/kg 0.472 

Vessel end material cost, $/kg 0.481 

Vessel fabrication cost, $/cu m 6l5. 

Vessel fabrication cost, $/kg 2*45. 

Tube material cost coefficient 0.0280 

Tube material cost, $/kg 4.4l 

Tube fabrication cost, $/kg 11.9 

Tube fabrication cost, $/sq m 964. 

Storage material cost, $/kg 0.441 

Insulation cost, $/kg 0.287 


TABLE II 

COST CALCULATION OUTPUT DATA 


Tube material cost, $ 4100 
Total tube cost, $ 19200 
Storage material cost, $ 7940 
Insulation cost, $ ll80 
Shroud side material cost, $ 199 
Shroud side total cost, $ 377 
Shroud bottom material cost, $ 138 
Shroud bottom total cost, $ 641 
Shroud top material cost, $ 62 
Shroud top total cost, $ 346 
Shroud total material cost, $ 400 
Shroud total cost, $ 1360 
Vessel material cost, $ I 65 O 
Vessel total cost, $ 3940 
Total materials cost for TES, $ 15300 
Total fabrication cost for TES, $ 18400 
TES total cost, $ 33700 
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c 0 N C L UDINQ REMARKS 


The use of a phase change thermal storage media for the 
storage of heat from solar collectors for use in electric 
power generation may have economic advantages over other 
types of heat storage media. Previous experience with 
phase change thermal storage has been limited to electri- 
cally heated applications in which the stored heat is 
used for low temperature water and space heating. Devel- 
opments of this kind have only been carried to point of 
small scale field testing of water heaters and space 
heaters . 

The work reported herein has been undertaken to provide 
an assessment of the characteristics of a phase change 
thermal storage unit meeting the requirements of a mid- 
temperature range solar powered electricity generating 
system. 

An analytical model and computer program have been de- 
veloped for predicting the performance and estimating 
the cost of a phase-change TES unit using a passive heat 
exchanger. These have been applied to a phase change 
TES medium having the composition 


NaOH (Commercial grade) 91.8?? (wt) 

NaNO^ 8.0 

Mn02 0.2 

The heat transfer fluid for charging and extracting heat 
is Therminol-66 . 

By relatively simple modifications, the model and compu- 
ter program can be applied to any phase change medium 
for which temperature-enthalpy data are available. 

The computer-predicted performance has been correlated 
with experimental data obtained from tests on a TES unit 
and test bed which have a thermal storage capacity (in- 
put and output per cycle) of 0.3 x 10° kJ at a maximum 

Input rate of 0.l8 x 10 kJ/hr, and a maximum extraction 

rate of 0.10 x 10 kJ/hr, and contains 1800 kg of the 
TES medium. 

Satisfactory agreement has been obtained between the 
computer predictions and the experimental data. 
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The program is now believed to be a reasonably reliable 
tool - for predicting the effect on the operation of phase- 
change TES units of this type, of variations In the 
operating parameters, and variation In size and storage 
capacity. 

The program has been used to predict the characteristics 
of a TES unit which has the following performance speci- 
fications : 

g 

Heat storage capacity of 3.1 x 10 kJ, oper- 
ating over the range of 516 K to 584 K. The 
heat transfer fluid for charging and discharg- 
ing heat Is Theimlnol-66. The maximum charge 
rate Is 1.8 x lO® kJ/hr with the Thermlnol-66 
temperature at 584 K ± 2 K. The maximum dis- 
charge rate Is 1.0 x 10“ kJ/hr with the Therral- 
nol-66 outlet temperature at 582 K ± 2 K and 
Inlet temperature at 516 K ± 11 K. 

The analysis showed that supplemental heat of about lOJS of 
the heat extracted from storage, added to the Thermlnol-66 
during discharge to help maintain the outlet temperature 
within the specified range, had a large beneficial effect on 
the size and cost of the storage system. Therefore, such 
supplemental ("deficit") heating was employed In the design 
of the system. 


The characteristics of a TES unit employing the NaOH-NaNO, 
heat storing medium predicted by the analysis and compu- “ 
ter program are as follows: 


Preliminary Design 


Vessel height and diameter 
Insulation thickness 
System height and diameter 
System weight 

Amount of storage material 
Unit cost 

Heat exchanger 


Figure of merit based upon 
a 25-element grid 

Storage output 
Specific cost 
Insulation loss 
Deficit heating 


2.35 m (92.5") 

0.609 m (24") 

3.6 m (142") 

27290 kg (60,200 lb) 
18000 kg (39,700 lb) 
$33,700 

.00635 m (1/4") tubing; 
280 colls, each 27.73 m 
(91 ft) long 

.882 

3.01 X 10^ kJ 
$11.2/MJ 

3.25/t of storage output 
10.2% of storage output 
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Report on 

THE THERMOPHYSICAL PROPERTIES OF 
A THERMAL STORAGE MATERIAL 


Port Comstock and Wescott, Inc. 

765 Concord Avenue 
Cambridge, Massachusetts 02139 


Section 1 


Introduction 

The work described In this report was carried out at Dynatech R/D 
Company for Comstock and Wescott, Inc. under Contract No. 23385. 

The purpose of this Investigation was twofold. The first task was 
to experimentally determine the effect of accidental contact between Thennlnol-66, 
a heat transfer fluid, and Thermkeep, a potential thermal energy storage material 
over the approximate temperature range of 20 to 480C (70 to 900F). Once a com- 
patible material was chosen, measurements of various properties using standard 
equipment and techniques available or de'*eloped especially for energy storage 
materials were utilized to generate baseline data on Thermkeep. The properties 
analyzed were differential thermal analysis, thermal conductivity in the solid 
and molten phases, density in the solid and molten phases, specific heat and heats 
of transformation, and viscosity. The temperature range of interest was 20 to 480C 
(70 to 900F). 
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Task 1 


Section 2 

PROPOSED PROGRAM FOR THE 
COMPATIBILITY OF THERMINOL-66 AND THERMKEEP 

Dynatech proposed the following study fbr determining the compatibility 
of Theminol~66 (Monsanto Chemical Co.) with Thermkeep. Design and fabricate a 
special mild steel closed container vented to the outside* fabricated to contain 
the salt and Therminol-66 . The ThermlnolobS sample will be analyzed using a gas 
chromatograph prior to being heated in contact with the Thermkeep. The container 
will be heated at a fixed rate to one of three or four successive increasing 
temperatures up to the order of 300C. During the heating and at each temperature 
condition any vapors given off will be collected and analyzed. After the heat- 
ing the salt sample will again be analyzed. Infra red spectra of the Thermlnol-66 
will be obtained before and after heating with the Thermkeep. IR spectra of the 
condensed vapors will also be obtained. Following this the container will be 
heated once more to the next higher temperature and the procedure repeated. 

Thermlnol-66 is stated to be a hydrogenated terphenyl with a boiling 
range of 340‘C to 396*C. Thermkeep is stated by Comstock and Wescott to be 
largely anhydrous caustic soda (technical grade) containing variable amounts 
of sodium carbonate, 8% by weight of sodium nitrate, and 0.22 of Manganese dioxide. 
It is said to melt from 282"C to 292*C. 
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Section 3 

EXPERIMENTAL PROCEDURE 

A container was made from a long nipple cut from standard 

2-*l/2" MPT iron pipe, threaded into a 2-1/2" NPT iron pipe cap with a 3" x 3" 

X 1/8" mild steel plate welded to the top of the nipple. A 3/8" - 3/4* reducing 
bushing was welded to the center of this plate and adapted from 3/8" NPT to 
3/8" OD stainlesA steel tubing with Swagelock fittings. The resulting vessel 
is 4-3/4" high as measured from the too of the steel plate to the bottom of 
the 2-1/2" IPS pipe cap. The stainless steel tube passes into a 500 ml 
ErlenmeTer flask submerged in a trichloroethylene-dry ice mixture. The flask 
is vented through Tygon tubing into a jar of water so any rapid evolution of 
gases can be Bxmitored. The apparatus was later modified by incorporating 
fittings for a second 3/8" OD stainless steel tube so the vessel can be 
flushed with inert gas during the run. 

The vessel was placed on a 4-1/2" diameter resistance heater which 
was mounted on a heat sink. Tmnperature instrumentation was spot-welded on the 
container to monitor the temperature through the sample / container assembly. 

The entire assembly was surrounded by a cylindrical furnace with the interspace 
filled with expanded vermiculite to minimize tmnperature gradients. 

Figures 1 and 2 are respectively pictures showing the vessel sur- 
rounded by the Heaters and the modified vessel connected to two lengths of 
3/8" OD stainless steel tubing in the oven with the upper heating element 
removed. The results of the compatibility tests are summarized in Table I. 

The results indicate the following: 

The chemical reaction between the Thermincl-66 and Tliermkeep at 
temperatures up to 320*C (close to the boiling point of the Therm- 
inol-66) are mild and produce no visible evolution of gases. Such 
reactions as do occur (darkening and thickening of the Therminol) 
should not have a detrimental effect of the Thermal Energy Storage 
(TES) system. 
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When a dry inert gaa (helium) is circulated above the Therakeep* 
it sweeps out any Therminol-^h present by the time the temperature 
of the Thermke>3p has reached 326"C. 
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Section 4 

FURTHER ANALYSIS OF THE THEBMINOL-66 AND TUBtMINOL-66 TESIDUES 
4«l Gas Chromatography 

After 9-1/4 hours a very sharp peak for Thermlnol-66 was obtained 
in the first run using our Cow Mac 69-750 Gas Chromatograph with a 6* long 
1/8" stainless steel column packed with 3Z OV-101 silicone oil on chromosorb H 
obtaincsl from Applied Science Laboratories State College, PA. The column was 
heat«i to 290"C. He were unable to obtain a peak on subsequent runs and con- 
cluded that either the column had been inactivated by long exposure at 290”C 
or that the peak observed was due to "noise” of sane kind. 

A second column was obtained only 2* long packed with Tenax 6-C 
on 60/80 Chromosorb H. This column is claimed to be capable of prolonged 
heating up to 350^C. However, we ob*-ained no Therminol peak during 14 hours 
operating at 300°C with maximum nitrogen flow. 300“C is the maximum temperature 
available with our Gow Mac 39-750 apparatus. We concluded that our gas chro- 
motagraph is net capable of analyzing Therminol-66 for breakdown because of 
its low vapor pressuz'> at 300°C. 

4.2 Viscosity Changes 

A qualitative assay of the viscosity of Therminol-66 at.d Thermlnol- 
66 from runs 21745 and 21746 summarized in Table I was made by putting a drop 
of each in three flat dishes, simultaneously tipping the dishes and observing 
the comparative rate of flow of the drops. From this, we estimate the viscosities 
to be as follows: 




CPS 

Therminol-66 as received 


150* 

Therminol-66 heated to 320®C 

in contact with Thermkeep 

225 

Therminol-66 heated to 320“C 
with Thermkeep 

for 16 hours in contact 

300 


* Monsanto data 
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4.3 IR Spectra 


Figures 3« 4, 5, and 6 are respectively IR spectra of: 

Figure 

3 Heptane 

4 lOZ Thermiiu>l-66 in Heptane 

5 lOZ TherDinol-66 taken from Experiment 

21745 (Table I) in Heptane 

6 lOZ Thenainol-66 taken from Experiment 

21746 (Table I) in Heptane 

All non heptane peaks observed in Figures 4, 5, and 6 are typical 
of aromatic hydrocarbons. There are no peaks indicating any take up of oxygen. 

4.4 Conclusions from Analytical Data 

The changes that take place in Therminol-66 when placed in contact 
with Thermkeep at 320’’C in contact with air appear to be slow darkening and a 
slow increase in viscosity. Any oxidation of the Therminol is too low to show 
up in the IR spectra. 
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Section S 


CONCLDSIONS 

We conclude from our study that accidental contact between Therm- 
inol**66 and Thermkeep in the TK will not result in a reaction which will be 
detrimental. 
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Section 6 


RECOMMENDATIONS 

Provide a means for sweeping the top of the Themkeep with a dry 
inert gas. This will sweep out any water present in the Thennkeep and also 
any Therminol-66 which may accidentally leak into the Thormkeep in the TES 
unit. The end of the vent should be lower than the exit point and cool enough 
so that any Therminol>66 will condense and collect in a vessel positioned 
below the end of the vent. 

Our experiment demonstrating the above was made with dry helium. 
Dry air or nitrogen should work equally well but if they are to be used tests 
should be made similar to the one using helium to confirm this assumption. 
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8IMMRY OF THERMINOL 06 > THEmKEEP COHPATIBILITY TESTS (1) 


• 

Run 

No. 

Weight of 
Therakeep 

K 

Weight of 
Thenilnol 
« 

Hok. 

Temp. 

•c 

Time at 
Hax. TMp. 
Hours 

Inert 

Cas 

Blanket 

Volatiles 

In 

Trap 

Comments 

27142 (2) 

100 

10 

273 

>1 

No 

Hooe 

A sample of Themlnol ohtainec 
from Che vessel using a hypo- 
dermic needle after the run 
appeared similar to the 
original material except for 
a slight Increase in viscosity. 

21743 (2) 

u 

u 

238 

>1 

Mo 

None 

No sample could be obtained. 
Thetmkeep too solid to permit 
sampling fluid. 

21744 (2) 

•• 

M 

310 

>1 

Mo 

None 

The Thermkeep had apparently { 

melted and solidified. 

21743 (2) 

M 

M 

320 

>1 

Mo 

None 

After cooling the apparatus 
was disassembled. The Therm- 
keep had solidified with a 
seml-solld brown mush of 
Thetmlnol and Thermkeep on the 
Surface. There was a small 
amount of viscous liquid con- i 
densate In the lowei end of | 

the stainless steel tube. | 

Samples were taken for further ' 
analysts and the vessel cleaned 1 
out by Immersing in hoc water. i 

21746 

.« 

20 

320 

16 

No 

None 

The vessel was cooled ard { 

examined as above. The results > 
were similar except th.«t the ' 

mush was darker brown in color. | 
Samples were taken for further 
analysis. | 

27130 

to 

«g 

320 

>1 

Tes (3) 

Ves (4) 

! 

After cooling the vessel was I 

opened. The Thermkeep had 
solidified to a solid blue 1 

mass with no evidence of any 1 

Thermlnol remaining. j 


NOTES: 


(1) Se« CMC for description of apparetus used. 

(2) The vessel was not opened after these runs. 

(3) Modified apparatus used with a flow of hellusi passing over contents and through condenser and 
Into water so flow rate could be monitored. 


(4) The condensate was the color and viscosity of Thermlnol 66. There was a few drops of 
water suspended in the condensate. 
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Figure 1 


Core of Heating Apparatus SRuwliillllie Rl-Thenukeep 


Corapatlbility Test Fixture In Place 



Figure 2 

Therffllnol M**Therffikeep Compatibility Test Fixture Modified to Permit 
An Inert Gas Sweep. Lower PartiRf: Heating Apparatus Shown Also. 
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Figure 4 10% Therminol 66 in n-Heptane 
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Figure 5 10% Thennlnol 66 in n-Heptane (Therminol 66 from run 21745 heated to 320“C) 
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Ta^k 2 
Section 7 

PROPOSED PROGRAM 

Dynatech proposed the following techniques for the measurement of 
the required therraophyslcal properties: 

a. The measurement of differential thermal analysis utilizing a Netzsch 
DTA with a Netzsch Temperature Programmer was undertaken on three 
replicate samples of the material with each sample being heated and 
cooled at an agreed heating and cooling rate. 

b. The measurement of thermal conductivity in the solid and molten 
phases by the Comparative Flat Slab Technique as modified for the 
requirements of the material using a Dynatech Model TCFCM-N8 Thermal 
Conductivity Instrument. 

c. The measurement of density and coefficient of thermal expansion In 
the solid phase by dllatometry utilizing a Netzsch Electronic Auto- 
matic Recording Dllatometer and the measurement of density in the 
molten phase by standard buoyancy techniques. 

d. The measurement of specific heat and heats of transformation In the 
solid and molten phases using a Dynatech Model QTA-N7 Quantitative 
Adiabatic Calorimeter. 

e. The measurement of viscosity utilizing a Brookfield Viscometer (with 
an extended splndlf*). 
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Section 8 


EXPERIMENTAL PROCEDURES 

8.1 Differential Thermal Analysis 

A amall amount of the crystalline material was placed in a specially 
shaped nickel DDK measuring cup centered around a protected thermocouple. A 
similar amount of powdered Kaolin as a reference material was similarly placed 
In a comparison measuring cup. The outputs of the two thermocouples were 
connected differentially. A separate temperature measurement protected thermo* 
couple was placed centrally between the two cups. The cups were placed at the 
center of a resistance heated temperature enclosure which could be heated or 
cooled at constant temperature rates. After allowing the system to come to 
equilibrium it was heated automatically at a constant rate of 5C/min up to 230C 
(450F). The heating rate was then decreased to 0.1 C/min up to 350C (660F). 

The heating rate was then increased to SC/min to 480C (900F). The system was 
then allowed to cool to 350C at SC/min and to 230C at O.lC/mln. 

The melting point and temperatures of transformation were determined 
from the intersection of lines drawn through the portions of the temperature 
time curve immediately below and above the peak. 

The system had been pre-calibrated with standard reference materials 

KCIO, , S10_, and K.CrO, obtained from National Bureau of Standards. The meas- 
4 2 2 4 

ured melting points were within + 0.5C of the reference melting points. 

8.2 Thermal Conductivity 

8.2.2 Solid Phase 


The basic method is shown schematically in Figure 7. For the 
measurements Pyroceram 9606, an inert ceramic glass, was chosen as the most 
suitable available reference material to cover all of the materials over Che 
total temperature range. The thermal conductivity of Pyrocc.'oro 9606 is in 
Che range 4 to 3 W/mK over the temperature range of interest and the salta 
were expected to be in a similar thermal conductivity range. 















For the purpose of test, fine gauge chromel/alumel thermocouples 
in a high nickel alloy protective sheath having overall dimensions of 0.25 mm 
were fitted tightly into the fine grooves cut across the surfaces of the sample. 
The instrumented sample was sandwiched between two similar instrumented samples 
of the reference material and of similar dimensions. 

The composite sample was placed between two similar sized heater 
units and mounted on a fluid cooled heat sink under a uniform load applied at 
the top of the composite stack. A metal and ceramic heater guard tube some 
75 mm diameter and 100 mm long which could be heated at various positions 
alongs its length and cooled at the bottom end was placed around the composite 
test stack such that the sample stack was positioned centrally with the guard 
tube. A further large metal fluid cooled shroud was placed around this stack 
and the interspaces around the sample and between the metal tubes filled with 
a thermal insulating powder which had been pre-dried under vacuum at elevated 
temperatures. A lid was fitted over the metal shroud and a large glass bell 
Jar was placed around the whole assembly. The system was then evacuated and 
backfilled with pure dry argon. 

A steady temperature equilibrium was established in the system by 
means of adjustment of the power to the heaters in the sample stack and the 
rate of flow of cooling water through the heat sink. The temperature gradient 
along the length of the guard tube was matched approximately to that on the 
composite test stack by automatic control of the heater along the length of 
the guard tube. 

At equilibrium the temperatures in different sections of the stack 
were obtained from the various thermocouples In different locations In the sys- 
tem. The thermal conductivity was derived from a knowledge of the heat. flow as 
determined from the mean value calculated in the top and bottom reference 
materials, the temperature difference across the sample and the known dimen- 
sions as follows: 
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2 (a 

sample 

^*Xt’ 1 - ^ fkAl 

-I bottom reference 




- 


«ffaere A « thermal conductivity 

Ax ■ thickness 
at » tenperature difference 

k *= thermal conductivities of the reference at respective mean 
temperatures 

Measurements were made at regular temperature intervals up to as 
close to the melting point as could be obtained without any part of the test 
sample being at a temperature above the measured melting point. Following 
these measurements a repeat determination was made at a lower temperature to 
check whether the sample had changed during the heating or for possible contam- 
ination of the thermocouples. 

8.2.2 Molten Phase 


The method chosen was based on Che comparative technique modified 
for Che molten state. It is shown schematically in Figure 8. 

Essentially the test sample is contained w. thin a thin walled 
cavity in a larger piece of a reference material. The bottom solid section 
of the material acts as the lower portion of the comparative stack while an 
upper section of the same reference material fits tightly into the sample 
cavity to seal the test sample within the cavity. The upper reference material 
contained a vent hole to allow for expansion of the test sample on melting 
and subsequent increase in temperature. Sheathed protected fine gauge thermo- 
couples were placed in the wall of the cavity, in the surfaces of the refer- 
ences touching the molten material, in holes along the length of the reference 
materials and in the sample itself. In all other respects the experimental 
techniques were the same as for the solid phase. 

In deriving the thermal conductivity of the molten sample the same 
principle was used as that for the solid material. However allowances had to 
be made for the heat flow through the walls of the sample container. This was 
found to be In the order of 20 to 25Z of the total heat flow measured in terms 
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Figure 8 

Molten Thermal Conductivity Test Stack Schematic Assembly 
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of that in the upper and lower reference materials. 

The test cell was fabricated out of Inconel 600. The Inconel 600 
has a much hlguer conductivity (10 to 20 than the salt, but this factor 

was compensated for by using a thin wall section and therefore decreasing the 
area ratio betwen the containment vessel and the specimen. The specimen was 
cast directly in the Inconel 600 container by melting the specimen externally 
and pouring It into the cavity. 

8.3 Linear Expansion and Density 

8.3.1 Solid Phase 


For the measurement of linear expansion the initial length of a 
smaple was measured accurately with a micrometer. The sample together with 
a tmnperature measurement themoucouple was placed in a calibrated fused 
quartz push rod measuring system of a Retzsch Electronic Automatic Recording 
Eilatometer. The system was placed at the center of a resistance heater envir- 
onmental chamber and allowed to equilibrate. Power was then supplied to the 
heater in a regular manner such that the sample temperature Increased at a 
constant slow rate 6f 1 C / min up to as close to the melting point of 
the material as could be obtained before sottenlng of the sample was noted. 
During the whole length of the experiment the continuous length and tempera- 
ture changes of the sample were recorded. 

The linear thermal expansion of the saroole was obtained ae fallows: 


L 

«C » " - 

Lo T 

where « & the coefficient of linear expansion 

Lo B the Initial length at 20C 

L = the change in length for a particular temperature Interval 

obtained from a curve of al versos temperature, allowance being 
made for the length change of the quartz system calibrated 
previously 

A-?S 
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The volume coefficient was obtained as the sum of the three respective linear 
coefficients and the density calculated from knowledge of the volume coefficient 
and the mass of the specimen. 

8.3«2 Molten Phase 

The basic buoyancy method was used to determine this property. A 
solid cylindrical piece of nickel some 6 mm diameter and 18 mm long weighing 
approximately 40g at 20C was used as the suspended bob for the measurements. 

The expansion of a separate rod of the same nickel was measured in the Netzsch 
dllatometer in order to determine the change in volume of the metal bob at any 
particular tonperature vrithln the desired temperature range. 

A sample of the material was melted into a large steel container 
which could be maintained at any desired temperature level. The nickel bob 
was suspended from tae arm of a balance into the molten salt. Shields and 
insulation were placed around the wire and between the furnace and balance to 
minimize radiation, convection and updraft effects. The mass of the bob was 
evaluated once the system had attained equilibrium. 

Measurements were made in this manner at successive Increasing 
temperatures some 200C (360F) above the melting point. 

The density of the salt was determined as follows; 

0 ^Ni^“20 ^ V 

® ”20 

where Ps“ density of salt 

= density of nickel bob at test teinp-»rature (T) 

*"20 s mass of nickel bob at 20C 

ntp “ mass of nickel bob at test temperature (T) 


A-26 


20 



8. 4 Viscosity 


600 grams of Thermkeep were melted In the vessel described in Task 
1 after modifying the vessel by opening up the center port and welding a 4“ 
length of 1 Inch IPS Iron pipe In place. A slow stream of a dry nitrogen 
gas was used to blanket the salt while heating. The gas was turned off 
during the actual viscosity measurements. 

The viscosity was measured using a LUT Brookfield Viscometer with 
a # 1 spindle at 60 RrH. Two spindle extensions were required to have enough 
length to permit the spindle to be immersed in the molten Thermkeep to the 
notch on the shaft and still have the instrument clear the top of the 1” pipe 

The Brookfield Viscometer was calibrated using glycerine the vis- 
cosity of which was measured using a calibrated Cannonubbeholde Viscometer 
#400 A 35 at 35*C. 
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Section 9 


TEST RESULTS 

9.1 Differential Thermal Analysis 


Table 2 

The Differential Thermal Analysis of 3 "Thermkeep" Materials In Air 


Test 1 

Test 

2 

Test 

3 

Heating 

Cooling 

Heating 

Cooling 

Heating 

Cooling 

O.lCmin”^ 

005 C mln"^ 

0.1 Cmin"^ 


0.1 Cmln"^ 

[X05Cmin“^ 

Endotherm 

Exotherm 

Endotherm 

Exotherm 

Endotherm 

Exotherm 

260 C 

277 C 

276 C 

280 C 

277 C 

280 C 

500 F 

531 F 

529 F 

536 F 

531 F 

536 F 

Endc *-.herm 

Exotherm 

Endotherm 

Exotherm 

Endotherm 

Exotherm 

306 C 


290 C 

222 C 

287 C 

228 C 

583 F 


554 F 

432 F 

549 F 

442 F 
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9.2 Thermal Conductivity 

9.2.1 Solid Phase 


Table 3 

The Thermal Conductivity of ''Thermkeep" in the Solid Phase 


Temperature 

Thermal Conductivity » \ 

C 

F 


-1 -2 -1 
Btu-in hr ft F 

50 

122 

1.30 

9.0 

150 

302 

1.10 

7.6 

250 

482 

1.08 

7.5 

150* 

302 

1.15 

8.0 


*on cooling 
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9.2.2 Molten Phase 


Table 4 

The Thermal Conductivity of "Thermkeep" 
in the Molten Phase 


Temperature 

Thermal Conductivity , X 

T 

C 

F • 

Wm“\"^ 

Btu-in hr ^ft ^ 

350 

662 

0.90 

6.2 

400 

752 

0.92 

6.4 

450 

842 

0.935 

6.5 
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9.3 Linear Thermal Expansion 


Table 5 

The Linear Thermal Expansion of "Thermkeep'' In Three Mutually 
Perpendicular Dir ‘ctlons 


Specimen 

4 

Thermal Expansion » (AL/L^) x 10 @ 

20 

50 

100 

150 

200 

250 C 

1 1 CO Laminations 

0 

6.0 

19,2 

37.0 

60.0 

90.0 

J_ CO Laminations 

0 

8.5 

25.2 

45.5 

69.5 

98.8 

1 1 CO Laminations 

0 

m 

19.0 

36.6 

59.2 

89.0 
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9.3.1 The Density of **Thermkeep** ic* the Solid and Molten Phases 


Table 6 


1 Temper 

ature 

Oens 
(Test Sp< 

ity 

scimen) 


ity 

itlcal) 

C 

F 

kg m“^ 

lbs 

-3 

kg m 

.3 

lbs ft ^ 

20 

68 

1883 

117.5 

2140 

133.5 

SO 

122 

1879 

117.2 

2136 

133.3 

100 

212 

1871 

116.8 

2127 

132.7 

150 

302 

1861 

116.1 

2115 

132.0 

200 

392 

1848 

115.3 

2100 

m 

250 

482 

1832 

114.3 

2082 

129.9 

300 

572 

1808 

112.8 

1808 

112.8 

350 

662 

1782 

111.2 

1782 

111.2 

400 

752 

1756 

109.6 

1756 ■ 

109.6 

450 

842 

1731 

108.0 

1731 

108.0 

500 

932 

1705 

106.9 

1705 

106.4 


^solld 

molten 


^melting point 


1.08 
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9. A Viscosity 


Table 8 


The Viscosity of ''Thermkeep" 


Temper! 

iture 

Viscosity, 

Centlpolses 

C 

F 

328 

622 

7.9 + 0.3 

398 

748 

6.7 + 0.3 

480 

896 

5.2 + 0.1 
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APPENDIX B 


COMPUTER MODEL USED FOR REFERENCE DESIGN STUDY 
AND PRELIMINARY DESIGN 


Main Program 


The main program determines the fluid condition which Is 
transferred to subroutine ELANAL, an element by element 
analysis of the heat transfer within the unit. The nature 
of the analysis Is such that It assumes that the temperature 
within the unit stays constant over a very small time. The 
size of the time step Is determined also In the main program. 

The time step Is computed as follows. Wf Is the computed 
fluid flow rate at the time of the calculation. Then: 





where 


z ■ 0.06 


^TK 



so that the time step Is determined by the ratio of the 
thermal capacity of the Thermkeep In an element to the 
thermal capacity of the T-66 passing through it. The 
ratio M^j^/Ng Is simply the Thermkeep mass per element. 

c^j^ varies throughout the unit as a function of temperature 
so that a c^j^ is selected to be as small as will be exper- 
ienced (liquid phase). Since the object is to compute a 
sufficiently small time step to satisfy the assumption of 
negligible temperature change in the Thermkeep, the time 
step will then be satisfactorily small for regions where 
c^j^ is larger and larger time steps could be tolerated. 

The factor of 0.06 was selected so that in the worst case, 
the error In heat transfer by this assumption will be no 
more than 2 %. Smaller time steps would reduce this error. 

A table of solar collector output is read into this sec- 
tion of the program. It is based upon key times within 
the twenty-four hour day, plus other tabulated points to 
fully describe the shape of the output. The key times 
are determined in such a way that the beginni ng of the 
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cycle occurs at the beginning of discharge of the unit, T6, 
i.e., the time when the solar collector output has 
dropped to exactly the boiler demand (late afternoon). 

(See Pig. 1.) 

From there to Tl, the collectors and storage combine to 
produce the desired flow to the boiler, 

sired flow rate to the boiler, is computed by 


'’b 

“min - Of(T„ - ’ 

where is the heating rate required by the boiler. 

The flow which can be provided by the collectors is: 

coll C(.(T„ - • 

where q^ is the heat output from the collectors, and the 
net flow from the storage unit during this period is: 


w 




W 4 

min 


w 1 -t 
coll 


At Tl, the collector output has dropped to zero, and all 
the flow required passes through the store ^e unit. There- 
fore, from Tl to T2: 


w 


f 


''min 


P. T2, the boiler demand goes to zero and a nighttir.e 
idle period (12 hours) occurs. At Tj, the boiler demand 
begins again (early morning). From T? to T3, >-he only 
calculations perform- h on the unit are those which are 
not flow related, i.e., aoss to the environment and axial 
heat conduction. 


:-2 



Prom T3 to T4, the collectors produce no output and all 
the boiler flow originates in the storage unit. For this 
period, as from T1 to T2, From T4 to T5, the 

collector output is less than the boiler demand and the 
flow, as between TO and Tl, comes from both storage and 
the collectors, the amount from storage again being 
determined by: 


w«=w. -w 
f min coll 


At T5, the collectors begin to put out excess heating 
over that demanded by the boiler. The excess heating Is 
computed by: 


q 


ex 



q 


B 


The flow to the storage unit during this period is calcu- 
lated by: 


«f = 


"ex 


and is limited to some predetermined maximum value. 

A key, called K7, is recorded to be transferred to ELANAL 
which tells it what the flow condition is so that the cal- 
culation process is directed to the correct part of ELANAL. 
If K7 = 2, then the flow enters che top or hot end of 
the unit and charging is occurring. If K7 = 3, then the 
flow enters the bottom or cold end of the unit and dis- 
charging is occurring. If K7 = 1, the flow is zero and 
no flow related calculations occur in ELANAL. This applies 
to the overnight idle and the points in the cycle where 
the collector output exactly equals the boiler demand. 

After determination of w^, /it and K7, the computation 

enters ELANAL and computes the internal heat transfer dur- 
ing the time seep. Upon completing the computations in 
ELANAL, the calculatic t returns to the main program where 
it adds incremental ' ' flows to the cumulative heat 
balances for the b^ heating, heat absorbed, etc. Upon 
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conviction of a dally cycle, the total cumulative heat 
balances are output, reset to zero, and the dally cycle 
begins anew with the state of the storage unit being as 
it was at the end of the previous cycle. 

Other information is produced such as peak flow rate, 
peak ideal pump power, peak p’^essure drop, and total 
elapsed operating time. The computation tests the num- 
ber of cycles against a preset number of daily cycles 
and either stops or begins analysis of a new design after 
the desired number of cycles has been completed. 

The heat balances are computed for each time step and 
summed to produce cumulative heat balances over the 
daily cycle. Q , the daily total solar collector output, 

is dependent upon how the time increments are determined 

since q is calculated by linear interpolation of a table 
0 

of vs. time and is assumed constant over a time step. 
Thus, 


cycle 

i 


^abs* amount of heat absorbed by the TE3 unit for the 

cycle, is computed by summing the enthalpy changes of the 
T-66 during each time increment of the charging phase. 


charge 

cycle 

«abs ' E “M '“l 

i 

Q, , the amount of heat which is wasted, i.e., which is 
xoss 

available frcm ohe solar collectors but not absorbable by 
the unit, due to reaching maximum pump flow, is the dif- 
ference between these. 


Q 


loss 


charge 

cycle 

E - “r,i 'r'-H 

i 
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Q . , , the amount of heat delivered by the storage unit 
del 

over a cycle Is also computed by summing the enthalpy 
changes of the T-66 during each time Increment during 
the discharge phase: 


Q 


del 


discharge 
cycle 

E « 


t,l °f'‘’'e,l " ^^1 


S)urn* compensatory or boost heating to be applied to 

bring the T-66, which exits at T^, up to for delivery 

to the boiler is the sura of the incremental deficiencies 
in each time step during discharge: 


Q 


burn 


dlscnarge 

cycle 



1 


CfCTfl 



Lt 


1 


Q , the total heat lest to the environment, is equal 

to the difference between Qg^jjg end once steady state 

is established and no net change of state is experienced 
within the unit over a cycle : 


Q = Q w - Qj 1 at steady state 
env abs del 


Subroutine ELANA L 


This subroutine is called by the main program each time 
step and performs the heat transfer calculations between 
the T-66 and the Thermkeep. However, one section of it 
is utilized only once for each design, the fl^-st time into 
ELANAL, and bypassed thereafter. This section perform.^ 
initialization of certain groups of values which are un- 
changed throughout the analysis, reads most of the inpvt 
to the program (geometry, cost factors, fluid properties. 



etc.), and conputea the estimated unit cost. It indirectly 
sets up a table of Therrakeep properties by calling sub- 
routine TKPROP which is used at this point to read in 
values and set up tables of element total enthalpies vs. 
temperature and solid fraction of Thermkeep vs. tetnpera- 
uure. These tables are referred to when properties are 
sought within the rest of ELANAL. 


Initialization Section 


The following calculations are performed to determine the 
geometry of the unit from input values. 

Displaced tube volume, 

\ubes * -T- Wo' 

where is the tube length, is the number of tubes, and 
d^ is the tube o.d. 

Thermkeep volume. 


V^K “ 


Total internal volume, tubes plus Thermkeep, 


V 

Total 


\r + Y 

‘tubes TK' 


An aspect ratio, the vessel Internal height divided by 
the internal diameter is specified so that, from knowledge 
of the total volume, the height and diameter of the ves- 
sel may be computed. 


shell 




h 
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If the height of the vessel happens, by poor selection, 
to be less than the total tube length, L^, then It Is 

set exactly equal to L-, and the diameter Is determined 
by the volume, * 


d ® 2 

shel!! 


The Inner diameter of the shroud containing the insulation 
is determined from the insulation thickness, and the 

vessel wall thickness, 


^shroud ” *^shell ^ ^ ^ins * ^ ^shell. 


The total shroud inner height is determined from the insu- 
lation thickness and the vessel and plate thicknesses, t^^^. 


^shroud ^shell ^ ^ ^ins ^ ^ ^end. 

The weigats of the vessel, shroud, heat exchanger, and 
insulation are determined from the density of steel (heat 
exchanger, vessel, and shrouds), density of the 

insulation, and the shroud thicknesses, ^ .bottom* 

^shroud, top* ^shroud, side* 

^vessel “ ’^^st ^end * '^shell ^shell ^shell 


'"^heat exchanger ~ ^steel ^i J ’ 
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w 


Ins 


IT 

ir p 


Ins 1 


^^shroud “ ( 


d 


shell 


+ 2t 




‘shell 


+ 2d 


2 

shroud 



> 


U 


shroud 


s irp 


k 


steel 1 shroud shroud shroud, side 


shroud 


+ 1 

V shroud, top shroud, bottom j 


The cost of the unit (without going Into each equation used) 
Is determined In two ways. Cost of materials Is based upon 
the material weight and a cost per unit weight. Costs of 
fabrication are determined by a well accepted approximation 
known as the "six-tenths factor” rule. This says that, 
based upon some parameter, weight, volume, surface area, 
etc., the cost of a particular part will vary as the 0.6 
power of the ratio of values of this parameter where the 
cost for one value is known: 


C = C 

unknown k lown 


0.6 


unknow n 

> 

known 


where C, is the known cost for a unit of value, P, 

known * known 


for the parameter used, and P 


is the value of that 


unknown 

parameter for the unit to be estimated. For example, a 
heat exchanger 3f twice the surface area of another nay 


increase only as (2)*^’^ or 1.5 cost. In most cases, 
the calculation uses both welr,ht and surface area as the 
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scaling parameters and averages the two values which were 
calculated. 


In this section certain groups of numbers which contribute 
to heat transfer calculations and do not change each time 
Incremei t are computed to save repetitious calculations. 

The total capacitive effect of the storage unit Is based 
upon the Thermkeep, the tubes, and the vessel. The insu- 
lation and the T-66 are ignored. Although the Insula- 
tion may have a large capacitive effect. Its thermal . 
diffusivity is so low that during daily cycling. It should 
be unimportant . 


The unit is divided into a number of numerical elements 
which are horizontal slices through the total unit. Each 
element, except for the first and last, have an equal 
weight distribution due to Thermkeep, tubes, and vessel 
sides. The weight of Thermkeep in each element, is 


the total weight divided by the number of numerical ele- 
ments. The weight of steel in all but the end element, 
AM , is the total tube weight plus the total vessel 

weight less the ends divided by the number of numerical 
elements. The weight of steel in the end elements, 

is ^M . plus the weight of an end plate for each, 
s u m 


9 


During the course of the analysis, the calculation works 
with the total enthalpies or energies of each element 
where 


E^(T) 

m 




AM . c ^ ,T 
st,m stee^ 


E^(T) 


AMTKhTj^(T) 


^st ,e^steel*^ 


These two tables of element enthalpies are established 
during the first pass through subroutine TKPROP when the 
enthalpy table for Thermkeep is read. In the initializa- 
tion section of ELANAL, the starting temperature is known 
and the energies are computed. Thereafter, new energies 
are determined from element heat balances and TKPROP de- 
termines the updated temperature of the elements. 
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Heat Transfer Section 


The unit la modeled as a vertical cylinder with the nu- 
merical elements being horizontal slices each assumed to 
be at uniform temperature. In this way, a step-wise 
temperature gradient Is calculated which approaches a 
smooth curve as the number of elements Increases. The 
computed temperature of each element Is taken to be the 
midpoint temperature for purposes of axial conduction 
calculations and plotting of the output. 

A heat balance Is made around each element consisting of 
three components; 

. . . the axial conduction from one element to 
the next consisting of conduction through 
the Thermkeep, through the tubes, and 
through the vessel wall; 

. . . the conduction to the environment which is 
primarily conduction loss through the 
thermal insulation; and 

. . . the heat transfer between the T-66 and the 
Thermkeep . 

The last component Is, of course, set equal to zero when 
no T-66 is flowing. The rates of heat flow by these 
modes are multiplied by the time step and summed, a net 
energy change per element is computed, and the new tem- 
perature Is found from subroutine TKPROP. 

The key, K7, selects which part of ELANAL is used for 
these computations, depending upon the direction of flow. 
The computations within each part are Identical but the 
direction of movement through the numerical elements is 
reversed causing a reversal of the indexing, and a change 
in the definition of which T-66 condition is specified 
and which Is sought. Consequently, the computational 
method will be described only once with differences between 
the parts being specified where necessary. 

The numerical elements are numbered from 1 to N„, with 1 

E 

being at the cold or bottom end and Ng being at the hot 

or top end. Fluid temperatures are numbered from 1 to 
Ng + 1. One is the lower manifold temperature which is 
the inlet during discharge when it is fixed at Tg and 
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Is the outlet during charging and has been called T^. 

1 4 Ng Is the upper manifold temperature which Is the 

Inlet during charging when It Is fixed at Tjj and In the 

outlet during discharge and has been called T^. 2 through 

Ng are the T -66 temperatures between the elements. Thus, 

for each element, T -66 Inlet and outlet temperatrres are 
defined, the outlet temperature from one el'^ment oelng 
the Inlet temperature to the next element In the flow 
direction. 

T -66 thermophysical properties, viz., specific heat, vis- 
cosity, denslt , and thermal conductivity, are assumed 
constant and are taken to be the valve at the nominal 
arithmetic mean temperature. The properties do vary with 
temperature but these variations are small over the tem- 
perature range of operation of the unit and accounting 
for this variation would add unnecessary complexity to 
the analysis. 

The calculations on the T -66 side are as follows. The 
fluid Reynolds number is defined as 


"re - “f 57^ 

where the flow area is 

^fl " f ^1* ^T* 


The T -66 side heat transfer coefficient is calculated as 


h 


f 


0.022 N 


0.8 


RE 


N 


O.h 


PR 



(Ref. 2) 
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Where N 


the Prandtl number is 


PR* 



M<.c 




In the laminar flow region, N_„ < 2000, the heat transfer 
coefficient is described by 


h 


f 


3.66 



« 


The object is to compute an overall heat transfer resist- 
ance between the T-66 and the Thermkeep. One component 
of this is the resistance between the T-66 and the inside 
wall of the tube which is 


hf 


where the total inside surface area, A^^, of the tubes in 
a numerical element is 


A 


i 


N 

^^iS N 


T 

E 


The fluid side pressure drop through the heat exchanger 
is : 


. 0.102f(N^) 



* In the final form of the analysis, this has been modi- 
fied to account for the effect of helically coiled 
tubes. See pages 72 and 7^ of this report. 
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where Is the friction factor and Is a user sup- 

plied function of the Reynolds number, and the coeffi- 
cient, 0.102, represents conversion factors to arrive 
at AP In units of kg/m*. fhe Ideal pumping power Is 


w- AP 

p 8 ■■ — 

P 102 pj. 


In kilowatts. The peak values of Ap and P are output 

P 

at the end of each cycle. The total cyclic pump work Is 


cycle 

% “ 3 ^ E ^P,l^^l 
1 

In kilowatt hrars. 

All calculations are performed In the mks system of units 
which Is essentially the SI designation. Mass Is ex- 
pressed In kilograms, length In meters, temperature In 
Kelvin degrees, time In seconds, energy In kilojoules, 
power In kilowatts, etc. 

The other components of heat transfer resistance between 
the T-66 and the Thermkeep are the ”lateral” resistance 
through the tube and the resistance between the tube out- 
side surface and the Thermkeep. The tube has an "axial" 
component also which Is associated with axial heat conduc- 
tion. The "lateral" tube resistance Is computed by the 
usual cylindrical conduction equation; 
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The only resistance calculated on the Thermkeep side 
Is the resistance of the solid layer built up on the 
tube surfaces. The amount of solid In a numerical ele- 
ment Is taken to be a function of the temperature of 
the Thermkeep and Is Input to the program together with 
the temperature variation of the enthalpy In tabular 
form. From knowledge of the solid fraction, the 
total solid layer size may be estimated. The volume of 
the solid layer around each tube In an element Is 


♦ AM-k 

^solld " * 


from which the diameter of the solid layer Is determined. 


^solld 



The heat transfer resistance through this solid layer Is 
also computed by the cylindrical conduction equation. 


In 


"solid 


2wk 


solid 


solid 




The three resistances are summed to determine the total 
resistance between the T-66 and the Thermkeep: 


^total “ "'f ^T,£ ^solld 
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The number of transfer units, NTU, of the elemental heat 
exchanger Is 


NTU «* 


1 

'^f ®f ^totaJ 


where the T-66 side capacity is taken to be the minimum 
capacity. This is in agreement with the assumption that 
the temperature change of the storage medium is negligible 
over a time step. With this assumption and the assumption 
of uniform Thermkeep temperature within an element, the 
effectiveness of the heat transfer process, e. Is calcu- 
lated by 


e 


I 


^-NTU 

e 


For a fluid with constant specific heat, the effectiveness 
represents the degree of approach of the fluid temperature 
to the Thermkeep temperature. Where Tj^^ is the temperature 

of the fluid into the element, is the temperature out 

of the element, and Tg is the temperature of the storage 

material in the element. 


T 

out 



+ e(Tg 



). 


The rate of heat transfer between the T-66 and the storage 
material in the element is 




WfCf 


(T 


out 


- 


Excluding the end elements, the •'lateral” insulation resist- 
ance to heat transfer is determined also by the cylindrical 
conduction equation. 
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c 


3 



and the rate of heat loss from an element at Tg through 

the insulation in the lateral direction to ambient tempera- 
ture, which is an approximation to the shroud surface 

tenq)erature is 


*^ins,f 


amb 


^InSjf 



The axial conduction term for all except the end elements, 
i.e., between the Ng storage temperatures of which there 

are (Ng - 1) conduction calculations, consists of conduc- 
tion through the heat exchanger tubes, the shell, and the 
Thermkeep. The resistances for each mode are computed 
separately and summed as parallel resistance paths. The 
tube axial conduction resistance is 


Rq., 


a 


f ^steel ^*^0 


- d,MN, 


while the shell axial conduction term is 


*^shell 


^shell 



^ ^steel ^shell shell 
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Both of these teinns are the same for every conduction 
path and for all times and are not changed after being 
set initially. 

The conduction through the Thennkeep is rather complex 
due to the geometry of the buildup of solid on the tubes. 
The conductivity in the solid is slightly higher than 
that in the liquid phase but not enough to warrant consid- 
eration of the above-mentioned geometry. Rather, an aver- 
age resistance within an element Is calculated assuming a 
net thermal conductivity which is the average of that in 
the solid and liquid weighted by the fraction of each. 


*^TK * ‘'solid + (1 - ♦> ‘'liquid* 


and these conductivities are available as temperature 
dependent functions within the program. 

The axial conduction resistance within a uniform element 
is calculated from the average cross-sectional area of 
the Thennkeep Itself and the conduction length as: 


shell 


TK,a 


?• Ic 

f “tk 




shell 


Nmd 
T o 




and is different for each element. Since the element 
temperature is considered to occur at the axial midpoint 
of each element for purposes of conduction, the actual 
conduction resistance through the Thermkeep between ele- 
ments is the series sum of the half values of resistance 
of adjacent elements. Between element J and element J + 1, 
e . g. , 




+ 


‘^TK.a 


(J) 


2 
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The total axial conduction resistance between tliese ele- 
ments would now be 




1 

1 — + 1 

*^shell ^,A 


+ 


1 


and the rate of heat flow In the axial direction is calcu- 
lated by 




Tg(J-fi) - Tg.(j) 
^TOT,a^*^ 


The total rate of heat flow associated with these elements 
is therefore 




For each of the end elements, one axial conduction term 
must be modified. Both of these elements are assumed to 
conduct in one direction to ambient through a series of 
two resistances. One is the half resistance of the ele- 
ment and the other is the net resistance between the end 
plate and the shroud which consists of contributions from 
the insulating material and the pipes and supports which 
penetrate the Insulation. The other axial conduction 
term for the end elements is the normal conduction to the 
adjacent element. 

The result of all of this is a net ^Q.pU) for each of 

the Ng elements. The total heat flow associated with an 

element over a time step is the net energy change of the 
whole element and is 
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These energy changes are applied to the total energies 
of each element to arrive at a new energy for each. 

These values of energy are used, by calling subroutine 
TKPROP, to determine the updated temperature array, 

Tg(j) and solid fraction array, for the storage 

unit. This updated state then becomes the steu?tlng state 
for the next time step. The state of the storage unit 
Is printed at key times as a table of Thermkeep tempera- 
ture, solid fraction, and T-66 temperature for all the 
elements. 

After solving for the new updated state of the storage 
unit and the exit temperature of the fluid, the calcu- 
lation returns to the main program for determination of 
the new fluid condition to be delivered to the unit. 


Secondary Subroutines and Functions 


In addition to the two major sections of the computer 
analysis there are three functions and one subroutine 
which are auxiliary in nature. These are used in sub- 
routine ELANAL for the purposes already described. 

Subroutine TKPROP is a routine used to determine the up- 
dated temperature and fraction of solid in a numerical 
element after the energy change over a time step has 
been applied. During the initialization section of ELANAL, 
TKPROP is called for the purpose of reading in a table 
of Thermkeep enthalpy and solid fraction as a tabular 
function of temperature, setting up tables of total ele- 
ment energy for internal and end elements as a function 
of temperature, and determining the starting value of 
energy for each element from the starting temperature 
by linear interpolation between the table values . There- 
after, the subroutine is used in the reverse mode to 
determine, as mentioned above, the updated temperature 
and solid fraction in each element as a function of the 
updated element energy arrived at in the elemental heat 
balances in ELANAL. This also is done by linear inter- 
polation within the tabulated values. 

Function TKCS is used to determine the thermal conduc- 
tivity of solid phase Thermkeep from the element tempera- 
ture. Function THCL is used to determine the thermal 
conductivity of liquid phase Thermkeep from the element 
temperature. 
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At the time of the preliminary reference design analysis 
the temperature dependence was not known and these 
values were assumed constant. 

Function FF Is used to determine the friction factor In 
the heat exchange tubing as a function of Reynolds 
number. It Is a curve Ht approximating the friction 
factors on the well-known Moody chart for an estimated 
relative roughness of 0.00012. 
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APPENDIX C 


ORIGINAL COMPUTER ANALYSIS 


Introduction 


Prior to the start of the current project, a detailed 
computer analysis was formulated. This analysis is simi- 
lar to the one generally used in the project. It does, 
however, attempt to account for additional effects not 
considered in the analysis described in the main text. 

It was not used for a number of reasons: 

... It is considerably more costly to use. 

... It is less flexible. 

... It would not as easily lend itself to 

the Incorporation of empirical modifica- 
tion. 

... It is based upon a published phase dia- 
gram for the binary system NaOH-NaNO^ 

which may not be characteristic of Therm- 
keep . 

Being made from commercial grade chemicals, Thermkeep con- 
tains impurities which probably alter its phase diagram. 

This analysis does compute for each element during each 
time step an overall heat transfer coefficient account- 
ing for the temperature dependence of all properties. 

It uses the phase diagram to include the effect of solu- 
bility of the two components in each other. It accounts 
for redistribution of mass as the density changes due to 
the solid-liquid phase change. The major purpose of this 
appendix is to describe the parts of this analysis which 
differ from the simplified version used in the report. 

This analysts consists of a larger number of functions 
and subroutines for the purpose of describing the tem- 
perature dependence of fluid properties and Thermkeep 
properties, and the characteristics of the phase diagram, 
and of handling the different conditions of Thermkeep 
distribution which occur in the numerical elements. How- 
ever, the basic program structure, method and order of 
calculations, etc., are the same in both analyses. 
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Description 


The TES unit Is assumed to be a vertical, cylindrical, 
Insulated vessel containing Therrakeep Into which Is 
Immersed a heat exchanger consisting of colled parallel 
tubes manifolded at the top and bottom. The analytical 
approach Is to divide the unit Into a number of numeri- 
cal elements which represent horizontal slices taken out 
of the cylinder. In the simplified analysis, density 
change Is not considered. Consequently, elements con- 
tain both fixed mass and fixed volume. Here, the volume 
of each element Is considered fixed but mass Is allowed 
to redistribute Itself. 

Figure Cl describes the way the Important part of the 
NaOH-NaNO^ phase diagram was approximated. Thermkeep Is 

roughly 8 % sodium nitrate. It Is assumed that the unit 
never operates at or above the melting point of pure 
sodium hydroxide. This can be changed but was not nec- 
essary for this project since NaOH melts at about 592 K 
(319 C) and the specified charging temperature Is 584 K 
(311 C). Thus, the orly process for liquefying NaOH 
Is that of dissolving Into the melt. 

The first euLectlc point for this system occurs at about 
33 % NaNO^. Consequently, the concentration of NaNO^ will 

never exceed this. Once enough NaOH precipitates out of 
solution, the eutectic behaves like a pure substance and 
freezes isothermally with no change in concentration. 

The numerical elements within the unit are assumed to be 
at uniform, constant temperature during a ^ime Increment. 
This was also done In the simplified analysis. The con- 
dition of the liquid Thermkeep Is determined by the equi- 
librium concentratlori at the element temperature and any 
excess NaOH precipitates out as a solid. Determination 
of the new equilibrium condition, however, requix -s an 
Iterative procedure. 

The cool-down process can be described In the following way. 
Assume one starts with every element in the liquid phase 
and occupying the full volume of the eloment. In fact, the 
size of the numerical elements is determined by the liquid 
phase volume. The net heat transfer by convection, axial 
conduction, and loss to the environment is computed gener- 
ally as it is in the simplified program with an exception 
to be described later. The program generally works with 
element energies also as in the simplified program. 
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An energy loss Is computed yielding a net energy at the 
end of a time step. The total mass of material In an 
element Is not adjusted until the end of a time step. 

The program keeps track of the amount of mass of each 
component (NaOH, NaNO^) In an element. The Llquldus line 

of the phase diagram defines the solubility of NaOH In 
the melt of NaOH-NaNO^. However, the amount of solid 

NaOH Is dependent upon the way mass has been transferred 
and the solid, of course, stays In Its element since It 
Is assumed to be frozen to a tube wall. The liquid phase 
is distributed to fill in any void space caused by solidi- 
fication with consequent higher density. 

Ignoring for the moment the way in which liquid Is re- 
distributed, the analysis takes the updated energy con- 
tent, knows the amount of each component as well as the 
steel associated with the vessel wall and the heat ex- 
changer and computes a new temperature. This temperature 
is one where the total of the energy associated with an 
equilibrium solution of melt plus the contributions from 
solid Thermkeep and steel agrees with the known element 
energy. Finding the proper temperature is an iterative 
process . 

The Thermkeep changes state in the following way. Pure 
liquid exists down to the Llquldus line. No volume change 
occurs. Between the Llquldus and the a-B transition, 
3-NaOH precipitates out as a solid assumed to adhere to 
the cool tube walls. Again, putting redistribution of 
liquid aside for the moment, the cooling continues to pre- 
cipitate solid B-NaOn until the a- B transition tempera- 
ture is reached. At this point, the melt is about 10% 
NaNO^ from which one can deduce the amount of solid B-NaOH 

in the element. 

The energy extraction process at the ouB temperature 
(566 K; 293 C) is then attributed only to the latent heat 
associated with the transformation of the NaOH not in 
solution in the melt. Once all the solid B-NaOH has trans 
formed to a-NaOH, the cooldown continues. The solid NaOH 
which precipitates out below 566 K (293 C) appears in the 
a phase. Thus from 566 K (293 C) down to 529 K (c56 C) 
(tt.e eutectic temperature), the analysis combines energy 
contributions from a-NaOH and an equilibrium melt to 
determine the updated temperature. 



At 529 K (256 C), the Thermkeep begins as a melt of eutec- 
tic concentration, 33% NaNO^, and a solid of a-NaOH. The 

eutectic melt then freezes as if it were of a pure material, 
l.e., isothermally , and no energy is removed from the 
a-NaOH or the steel. When all the latent heat associated 
with freezing the eutectic has been removed, the mater- 
ials then cool sensibly with a layer of eutectic composi- 
tion solid over the layer of a-NaOH frozen to the heat 
exchanger tubes. 

During heating, the reverse process is assumed to occur: 
sensible heating of the eutectic plus the a-NaOH solid 
phases; melting of the eutectic at 529 K (256 C); heat- 
ing of the solid a-NaOH and the melt with NaOH dissolv- 
ing into the melt; transformation of the solid a-NaOH on 
the tube walls to solid 3-NaOH at 566 K (293 C); further 
heating of the melt and the solid 8-NaOH with the NaOH 
dissolving into the melt; and, finally, sensible heating 
of the liquid phase Thermkeep up to 584 K (311 C). 

The volume change during cooling is treated by computing 
resultant void spaces in the constant volume elements 
and filling downward so that all the void space is trans- 
ferred to above the liquid level. At the end of a time 
step a new fraction of liquid and solid results from the 
new energy content and consequent new equilibrium tem- 
perature. An array of volumes of liquid is composed for 
the elements as well as an array of the total volumes in 
each element not occupied by solid. The liquid is then 
used to fill the dead space in each element taking liquid 
from sequential upper elements as needed to completely 
fill lower elements until all the liquid is accounted for. 
The revers-e procedure is used during heating where excess 
volume is generated by melting and must be transferred up 
to successive volume elements. 

The relative amounts of each chemical component in ele- 
ments of different temperatures are different and the 
actual amount of each being transferred is calculated to 
arrive at a new concentration of each in the melt at the 
end of a time step together with a new total energy con- 
tained in the element. In reality, this might imply a 
new equilibrium concentration after mixing but it was 
assumed that, since the time steps are quite short and 
since the transfer of material should be proportionally 
quite small, the resultant temperature of each numeri- 
cal element is unchanged and the subsequent time step 
uses the temperature distribution as computed at the end 
of heat transfer for the previous time step. Thus, the 
available AT for heat transfer is infinitesimally 
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Inaccurate t but the energy and material contents are 
accurately accounted for« so that In all respects the 
laws of thermodynamics of the process are obeyed. 

In the write-up on the simplified program. It was ex- 
plained that a criterion was used for computing the time 
step whereby the capacitances of the fluid and the stor- 
age medium were compared. The resultant time step would 
be small enough to Justify the assumption of negligible 
temperature change In the storage medium and the use of 
the simplified relationship. 


E - 1 
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to characterize the heat transfer process. Where no loss 
of mass is experienced, the use of a minimum value of 
c^jj was considered a sufficient criterion. 

In this analysis. It is possible to lose enough material 
In a particular element to cause this assumption to be 
inaccurate. When this is the case, the time step is still 
determined as before but for the elements at the upper end 
of the unit where insufficient mass may exist, calcula- 
tions are performed which account for the changing element 
temperature during the heat transfer process. The calcula^ 
tion is separated into sections of the phase diagram to 
include the knowledge of temperature arrests, etc. The 
general formulation for variable temperature is as follows 

The temperature of adjacent elements ought to rise or 
fall similarly to the one in question. Moreover, the 
axial transfer is assumed to be small. Therefore, the 
axial heat flows are based upon starting temperatures and 
are assumed to be constant during a time step. A heat 
balance around the element in question says: 


dT 

®"c^ dt 


E 


q^ - q^ - Wj.Cf.e (T^. - - 


'^E ~ '^amb 
^Ins 


where and q^^ are the axial heat flows computed in ad 

vance, the next term is transfer to the fluid, and the 
last term is Insulation loss. The calculation of e may 
be unchanged as long as the time variation of T^., the 
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element temperature, is considered. The solution to this 
first order differential equation is slmpl;y 


T 


E 



a 

/ 


At 


e 


X 


where 




_1 

^ins 


^2 " ‘^l R 


amb 


Ins 


+ Wj.Cj.e 


Em 

T B _.C. 
a 


and ) represents the total heat capacity oi‘ the stor- 

w 

age materials in the element. 

The computer analysis uses two subroutines for this pro- 
cess, one for heating and one for cooling, in order to 
properly move along the phase diagram. Comparing the 

g 

value of - with T- the initial element temperature, 
a D 9 o 

determines which of these to use. 

If use of this relationship over the whole time step. At, 
passes a transition temperature, such as a-6 transition, 
then a reduced At is computed to get to the transition 
temperature and the element is assumed isothermal at the 
transition temperature. The computations proceed this 
way until the whole time step. At, is traversed. Some 
simplifications are made in the way the capacitance term 
is handled but in the end the thermodynamic balances are 
always preserved and no serious errors arise. 

As one will now realize, the computations are considerably 
more involved than those In the simplified program which 
is why it is so costly to run and so hard to modify em- 
pirically. The following sections will try to describe 
each of the subroutines and functions used. At this 
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point, the program works reasonably well; but since It 
was never used In detail for design or test analysis. It 
Is possible that there are areas requiring modifica- 
tion which would be discovered through extensive use. 


Main Program 


The main program Is, In essence. Identical to the main 
program of the simplified analysis Already described. 

It does refer to a subroutine, COLL, to determine the 
solar collector characteristic but this Is an unimportant 
difference. 

It Is, as before. In this section where the time step, 
the fluid flow rate, and direction are computed. The 
cyclic heat balances are performed. All the data are 
printed. 


HTCALC 


This subroutine performs all of the heat balances except 
for the case where the element temperature is assumed to 
change during the time step. In that case, basic heat 
transfer Information Is still computed here but the calcu- 
lation is transferred to either VRTHS or VRTCS to deter- 
mine the energy change In the element. In certain cases 
where sufficient material exists to Justify it, the pro- 
gram uses a sequence of up to three calculations assum- 
ing constant temperature but with the time step divided 
and an updated temperature computed between each reduced 
time step. 

The calculation me chod is essentially identical to that 
used in ELANAL In the simplified program. The differ- 
ence is that the program keeps track of the masses and 
volumes of solid and liquid In each element and the masses 
and volumes of the chemical components, so that accurate 
accounting is made of the process of transfer of mass 
from one element to another. 

At the end of the energy balances, regardless of which 
method was used to accomplish them, the analysis has 
arrays of variables describing the material content of 
the elements. These are listed below with both algebraic 
and computer notation; 
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i * AL(I) ■ mass of liquid In element i(I) 

Ag ■ AS(I) « mass of solid in element i(I) 

^nf,l * ANL(I) * mass of NaNO^ (liquid) in element i(I) 

A^g « ANS(I) » mass of NaNO^ (solid) in element i(I) 

• EL(I) a energy of liquid in element i(I) 


These values together with knowledge of the liquid and 
solid densities (one value only is used for each, and the 
total fixed element volume can be used to compute the 
material transfer. The computations are as follows: 

The volume of solid in element i is 




s 


the volume of liquid is 




t 


The volume available for liquid is the difference between 
the element volume, V„, and the solid volume. 
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From this Information Is created a K-long array of avail- 
able liquid volumes and associated enthalpy and NaNO, 
concentration 


■’j ' 


where j * 2 is the first of the i array of elements with 
available void space, and the total available liquid vol 
ume is 


K 

^TOT,f * YL ^f,J* 

The total available void space (not solid space) is 

^TOT,V “ C ^v,i 


and the volume remaining at the top of the unit which is 
empty becomes 


The liquid is then redistributed in the following way. 
The values of void and liquid are compared sequentially. 
If the void is larger than the liquid, all the liquid is 
assigned to that void and the next liquid volume is 
examined, etc. until the value of liquid is larger than 
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the remaining void. Then the void is filled, the volume 
of liquid from which it was taken is reduced by that 
amount, and the next void is examined and compared with 
this ’-ext available liquid volume. 

The arrays AL(I), EL (I), and AKL(I) are then reevaluated 
as the sums of the contributions from each of the araountc 
of liquid deposited into the particular element. In this 
way, all the mass, energy, and chemical components are 
accounted for. The calculation, on the assumption that 
the transfer of material is small, does not readjust the 
element temperature to allow for the new mixed tempera- 
ture since it would involve new considerations of solu- 
tion equilibrium, etc., and new fractions of liquid and 
solid, requiring a repetition of the above procedure. 

This did not seem necessary. 


FULLEL 


This subroutine is used to determine the updated states 
of the elements when the normal calculation of heat bal- 
ances can be used in HTCAtC. It takes the newly computed 
energy of the element together with the total mass of 
material and the knowledge of the amount of NaNO^ in the 

element and computes the new equilibrium temperature and 
relative amounts of solid and liquid phase. 

The computation determines where the total energy of the 
element lies with reference to key identifiable points 
in the phase diagram. These points are: the Liquidus 

line temperature at the concentration of NaNO^ which is 

input, the ot-B transition temperature, with all B and all 
a points of the NaOH component; the eutectic temperature 
both all liquid and all solid eutectic component. Above the 
Liquidus line, between all a and all B“NaOH at the transi- 
tion temperature, between the all liquid and all solid 
eutectic points, and below the eutectic temperature, the 
calculations are straightforward since they involve 
either pure sensible or pure latent heat effects. 

Between the Liquidus temperature and the a-B transition 
temperature or between a -3 transition temperature and 
the eutectic temperature, an iterative approach is used. 

This is a standard method known as "the method of regula 
falsi" and is described in mathematics textbooks (Ref. 5 ). 
The calculation begins by assuming a straight line between 
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reference points and Iterates with successive straight 
line approximations until the error is small. 


VRTCS, \T^THS 


These two subroutines compute the element energy change 
when the amount of material in the capacitance of the ele- 
ment is small enough to require considerations of the 
variation in temperature during the time step. It assumes 
that only solid phase remains and consequently considers 
only pure sensible or pure latent heat transfer. 

In the case where some liquid resides in the element, it 
will underestimate the heat transfer but it is assumed 
that the error is not significant. Energy and material 
balances are maintained at all times. It should be noted 
that at any time there would be only one element contain- 
ing solid plus some liquid and it is only here, and only 
when the liquid is quite deficient in terms of filling 
the available space, when the error would occur. 

VRTCS is used during cooling and VRTHS is used during 
heating. The method of calculation has already been 
described. As discussed, the variable temperature analy- 
sis is used when sensible heating occurs. When a cross- 
over into a latent heat region occurs, the time step is 
divided and the calculation for variable temperature is 
used when needed. Two subroutines were written for 
control purposes to allow selection of moving either up 
or down in temperature in the phase diagram. 

These routines compute the energy change in the material 
directly, and secondarily compute the energy change of 
the fluid to determine the nominal fluid temperature out 
of an element. 

Most calculations are done in the reverse order where the 
fluid energy change is computed first and then assigned 
to the element. 


COLL 


This subroutine reads in the solar collector characteris- 
tic and is used to determine collector output in the main 
program. The simplified computer program does this 
directly in the main program. 
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RLCALC 


This subroutine was to be used to describe an additional 
heat transfer resistance between the T-66 and the Therm- 
keep which Is associated with the liquid phase trans- 
ferring heat to the solid on the tubes. At present It 
Is assumed to be zero and the subroutine Is not used. 


PP 


This function Is used to calculate the flowing T-66 fric- 
tion factor as a function of Reynolds number. It is 
used for pressure drop calculations. 


PLTEMP 


The temperature variation of the specific heat of the 
T-66 is considered in this analysis. During HTCALC, the 
specific heat at inlet temperature to an element is used 
to compute the heat transfer. Subsequently, the computed 
enthalpy change of the T-66 is applied to the inlet 
enthalpy and a corrected outlet temperature is computed. 
This function computes the temperature of T-66 as a func- 
tion of that enthalpy. 


CPF, VISE, ENTHF 


These functions compute the specific heat, viscosity, and 
enthalpy of the T-66 as a function of temperature. 


CPA, HA, HB 


These functions compute the specific heat of a-NaOH, the 
enthalpy of a-NaOH, and the enthalpy of B-NaOH as func- 
tions of temperature. The specific heat of B-NaOH is 
relatively constant so that the enthalpy is a linear 
function of temperature. The a-NaOH appears to have a 
temperature dependence but within the temperature range 
of Interest here it is not significant, so that for now 
this has been assumed constant and the enthalpy is a 
linear function of temperature. 
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THCA. THCB, THCE 


These functions compute the thermal conductivities of 
o-NaOH, 3-NaOH, and eutectic solid phase as functions of 
temperature primarily foi purposes of solid resistance 
to heat transfer between the T-66 and the liquid phase 
Thermkeep. At the present time these are assumed to be 
equal and constant. 


CLIQ, TLIQ, HLIQ. CPT 

The first two functions are the inverse of each other. 
CLIQ computes the concentration of NaNO^ in the melt as 

a function of temperature along the Liquldus line. TLIQ 
computes the temperature of the melt as a function of 
concentration of NaNO^ along the Liquldus line. The line 

Is assumed to be strai" , at present, as shown on Figure 
Bl, since the accurate xocus was not required at this 
point. 

Energies derived from the phase diagram are estimated 
based upon the appearance of the phase program and all 
compositions are assumed to have the same total enthalpy 
change distributed according to the phase diagram. HLIQ 
computes the enthalpy of the melt as a function of tem- 
perature for the Liquldus line. No solution effects are 
accounted for. Eventually, good h(T) data should be 
obtained for a range of concentrations so that a good 
map of h(T,C) could be produced. 

CPT computes the specific heat of the liquid phase as a 
function of concentration of NaNO^. This implies that 

no temperature dependence is expected. At this point the 
function does not provide even the concentration depend- 
ence since data are not available. 
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Conclusion 


At the completion of tho work on this computer analysis, 
Comstock & Wescott submitted an interim report, "Computer 
Program Development," dated December 7, 1976, to Sandla 
Laboratories (their P.O. No. 87-5030). 

This report discusses the general program philosophy and 
the way by which data are fed into the program. Sample 
output was presented. Decks of computer cards were de- 
livered for both the FORTRAN and the binary programs 
along with listings of the routines in both FORTRAN and 
machine language. Following this discussion are list- 
ings of the FORTRAN code for all the routines. 
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f%ftrtnn^r 


FORTRAN CODE 


ORIGINAL PAGE IS 
OF POOR QUALITY 


HEAT-OF-FUSXON THERMAL STORAGE ANALYSIS 

THIS PR06RAM ANALVEES THE DYNAMIC RESPONSE OF A TMO-COMPONcNT* 
HEAT-OP-PUSION* COUNTERPLOli NEAT STORAGE DEVICE 
INPUT HEAT 1N« IS DETERMINED OY A SOLAR COLLECTOR CHARACTERISTIC 
OUTPUT IS TO A CONSTANT DEMAND BOILER POP A RANKINE CYCLE TURBINE 


too PORMATC3PIO«8»F10«3«P10«E*P4«O*/PI0»3»Pl0«7*lJ»P10*3*3PI0.6l 

101 P0RMAT€4P10«4«13»P10«4| 

102 PORMAT(P10«B*P10«7«2P10«3I 

103 FORMAT! 7P10*3/3F10«6*P10«7| 

201 PURMATI1H1»SX«« HEAT^OF-PUSION THERMAL STORAGE ANALYSIS* »///) 

202 FORMATCIBX** 0ATA«»//| 

NO. PARALLEL TUBES 
TUBE X.O.. N. 

TUBE O.O.. M. 

TUBE LENGTH. M. 

TUBE NATL DENSITY. KG/CU.M. 

TUBE THERMAL K. KA/N.-K 
AMT STORAGE NATL. KG. 

VESSEL DIAMETER. M. 


203 FORMAT! SX. 

204 FORMAT! SX. 

208 FORMAT!SX. 
20o FORMAT!aX. 
207 FORMAT! SX. 
20« FORMAT! SX. 

209 FORMAT! SX. 

210 FORMAT! SX. 


211 FORM AT! //.SX.* NO. ELEMENTS 


212 FORMAT! SX. 

213 FORMAT! SX. 

214 FORMAT! SX. 

215 FORMAT! SX. 

216 FORMATI5X. 

217 FORMAT! SX. 
216 FORMAT! SX. 

219 FORMAT! SX. 

220 FORMAT! IHl 

221 FORMAT! lOX 

222 FORMAT! lOX 

223 FORMAT! 1 OX 

224 FORMAT! 1 OX 

225 FORMAT! lOX.* 

226 FORMAT! lOX.* 

227 FORMAT! 10X>* 
226 FORMAT I /// 
229 FORMAT!* !11 


MAX. TINE INCR.. SEC. 
BOILER INLET TEMP.. K 
BCILER OUTLET TEMP.. K 
MAX. PUMP PLOW. KG/SEC 
SOLAR COLLECTim OUTPUT NO. 
BOILER DEMAND. K« 

TOTAL MEI6HT. KGlAPPROXl 
VESSEL HEIGHT. M. 

9X.* ELAPSED TIME. NR 
YCLE TIME. SEC 
LOST COLLECTOR HEAT. KJ 
COLLECTOR OUTPUT. KJ 
STORAGE OUTPUT. KJ 
STORAGE INPUT. KJ 
AUXILIARY INPUT. KJ 
BOILER REQUIREMENT. KJ 


1S!1> 


ANL!II 


TE!I1 


»//l 


TF!1| 


*.F10.0> 
*.P10.4) 
*.F10.41 
•.F10.21 
*.F10.31 
*.F10.51 
*.FI0.2l 
*.F10.3I 
*.F10.0) 
*.F10.2I 
•.P10.2I 
*.F10.2I 
*.FI0.3I 
*.6X.14) 
•.F10.31 
*.FI0.1I 
*.F10.1| 
•.P10.21 
•.FIO.IA 
* .E10.4) 

• .£10.41 
• .£10.41 
• .El 0.41 
*. £10.41 
• .£10.41 
SYSTEM MAP*. //I 
VS! 11 


VL!I> 


AN 


230 FaRMAT!lX.I3.1X.P9.2.3X.P9.2.3X.P9.5.3X.P9.S.3X.F9.5.3X»F9.SI 

231 FORMAT! IX.I3.13X.F9. 21 

232 FORMATISX.* ESTIM. MATL COST. DOLLARS *.P10.21 

233 FORMAT!SX.* OVERALL HEIGHT. M. *.P10.11 

234 FORMATISX.* OVERALL DIAMETER. M. *. FI 0.21 

238 FORMATISX.* AMOUNT OF INSULATION. KG *.F10.2l 

236 FORMAT! SX.* AMBIENT TEMPERATURE. K *. FI 0.11 

237 FORMAT!////.* VS. VL.ANS.ANO ANL SHOULD BE MULTIPLIED BY THE NUMBEk 
I OF TUBES WHICH IS ••FlO.ll 

236 FORMAT! SX.* SHROUD THICKNESS. M. 

239 FORMAT! lOX.* PEAK PUMP POWER. KW 


* .FI 0.61 
• .clO.41 


240 FORMAT! lOX.* TOTAL PUMP WORK. KW-HR 

241 FORMAT!! OX.* PEAK PRESS DROP. KGF/SM 

242 FORMAT! iOX.* PEAK STORAGE PLOW. KG/S 
24S FORMAT! SX.* VESSEL WALL TH.» M. 


*•£10.41 
*•£10.41 
* .F10.4I 

*.F10.6l 


mCEDING PAGE BLANK NOT RLMO 



2«6 rafIMATCSX** VESSEL BASE TH«* M* ••FI0.6I 

2«7 FOBMATCSXt* INSULATION THICKNESS* M. •»F10«4| 

24S FOAMATiSX** TM* CQNO* INSUL** KM/M-K •*F10*6) 

290 FONMAT(SX*« LlOU* PHASE OENS* KG/CU*M* •*FI0*2| 

2S1 F0RMATI8X** SOLID PHASE OENS* K«/CU*M* •*F10*2I 

292 FORMAT! 8X*« INSULATION OENS** KG/CU*M* •*FI0*2) 

293 FORNATCSXa* H*T* FLUID OENS* K6/CU*M* •*FI0*2) 

COMMON TECIO 1I*TF( 101 1 *VSC 101 1* VLC t01l*ANSCl 01 1 *ANL(101 1 
CCMMON/OATA/RHOL*RHOS* TTR * TCU*HM0 *MEUT*HABT *CPEU *CPS * C V * TKL 
C0NM0N/VC0N/R«AX*R1NS*R1NEN* TAMB*HIT 

C COMPUTE THE UNIT GEOMETRIC CONFIGURATION 

REAOCSalOOl 0O*0I*0V*XLT*ASTOR*XN*RHO«*TK9*N*0TMAX*TVV*T9V*TSH 
HcA0(9*102l TINS*TNCIN*RHOF*RH01 

REA0CS*103) RMOL*RHOS*TTR*TEU*HMO*HEUT*HABT*CPEU*CP9*C«*TKL 
C EASIC UNIT GEOMETRY ANO FIXED HEAT FLOW RESISTANCES 

XM«M 

VSTORsASTOR/RHOL 

VTUEsO .78844 XN4XLT400442 

HITaC VSTOR4 VTuei/ C 0* 789440V442I 

XLTEbXLT/XM 

A1ETS3.1416940I4XLTE 

VcLsVSTUR/l XM4XNI 

AlFsO * 78S44C 00442-0 1442)«XLTE4RHO« 

RVsALCGCOQ/OI I/I6.28324TKM4XLTEI 
RVWSH1T4XN/€3*I4I8940V4TMV4TK«M>XM| 

RTU»*XLTE/f 0*78944(00442-01442|4TKV| 

RMAXsl./Cl*/RV44I./RTU8> 

R INSs T 1NS4XM4XN/ C 3. 141 594 THC 1 N»0 V*H1 T I 
R£N0«T 1 NS4XN/ CO. 78S44THC IN40 V«42 1 
RPENsT INS/C 0 .78844TKII4(004«2-014«£ I ) 
RlNENsl./Cl./RlNS41.yREN041./RPEN| 

AFLTsO. 7884401442 
AVHS3 .141 894RH0M40 V4TH V4HI T 
AVEL«AV4/ CXN4XMI 
A«sA44AVEL 

AVSsO .78S44RH044T8V40V442 
wesAVB/XN 
HSHSHIT42.4TINS 
CSHSCV42 .4TINS 

■ INS30.78S44RH014UOSH442-OV442I4HIT42.4T1NS40SH4421 
HSHsQ . 78844RH0H4TSH40SH4C 4.4HSH40SHI 
tlTF«AFLT4XLT4XN4RH0F 

H E I Ts A44XN4XM42 • 4 AVS4A ST0R4 4 INS4 4SH4 MTF 

COS T« 1 . 1 4( AM4XN4XM42 .4AV84IISH |40 .444ASTOR4 0 . 22 4« 1 NS 

REAOCS.lOll TU*TL*«MAX*CL1QR*N0*TAN8 

HUsENTHFCTUJ 

HLsENTHFfTLi 

OHMAX«HU-HL 

C COEFFICIENT FOR DETERMINATION OF TIME STEP 

2s0 .14VEL4CLI0R4RH0L/CPFC TUI 
Jls99 

CALL COLLCCT*OC*T1*T2*T3*T4.08.J1I 

«MlNaQ8/0HMAX 

MR1TEC6.201I 

MRITEC6.202I 

WR1TEC6.209I ASTOR 

HRITEC6*203I XN 

«R1TEC6*204I 01 
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ttiRlTeC6«20SI oo 
«<UTEC6*2a«l XLT 
MftlTE(6»207) RHO« 
llRITE<6«2S2) RHOI 
«RlTe(o«2S3l RHOF 
«RiTe(t>»£SOI RHOL 
MHlT«C6t2&i) RHCS 
WRITE16*249I THV 
WR1TEC6*246| TBV 
■RITE(6»20«) TKH 
«R|TE<4*238> TSH 
WR1TEC6*246) THClN 
WR1TE<6«2I«I HIT 
»R1TE<6*210I OV 
MRlTECo»247| TINS 
«RlTEl6»2i3) HSH 
MR1TEC6»234I OSH 
WR1T£(6*2J91 MINS 
MRiTE(6»218J HEiT 
MR1TEC6*232I COST 
WRlTECotEin XM 
MR1TEC6.212I OTNAX 
«R1T£<6*213I TU 
WRlTEf6*214| TL 
•R1 TEC 6« 2361 TAM6 
MR1TEC6*215I MMAX 
WRlTEC6*2t6> J1 
•K1TEC6«£17| OB 
NCsO 

ET«0. 

CT»T4 
TFC n»Ti. 

0B011.sQe«T4 

OO TO 3 

C ENTER SUBROUTINE FOR INTERNAL HEAT TRANSFER ANALYSIS 

10 CALL hTCALCC WT»6*0TlME*K3«M*01*Rfe«AI£T»AE»D0«XLTE*YEL» JSvIliBl 

C COMPUTATION OF PUMP REQUIREMENT FOR PUMP EFF»1 .0 

TLMs(TFCM4l)«TF41l)/AL0GCTFCM4l}/TFC II) 

REAS64C1/V1SFCTLM) 

P0«0 • 1 024FFC REA J4XLT4044 2/C 2 •«Ul*RHOF I 
lFCPO*GT«POMAX) POMAXsPO 
PUWERshAPO/C 102«4RHOFJ 
IFf P0HiER.6T«PMAX} PMAX«POHER 
HPUMPS2 . 723E-64P040T lME/RHOF-«>llPUMP 
iFC M.OT.MPEAKi «PEAK*« 

60 TO CS0«91I»K3 

C COMPUTATIONS ASSOCIATED HITH CHARGING PHASE 

SO QLOSS«QLOSS4C 0C«Q8-M4 C HU- ENTHFC TF 1 1 1 > >) 40TI ME 
QSOL«OSOL-»OC«CT IME 

QCH6>0CH6«M4CNU-ENTHFCTFC 1 )||40T1ME 

GO TO S2 

C COMPUTATIONS ASSOCIATED MXTH DISCHARGING PHASE 
b 1 QBURNsaeURN4ll«C HU-ENTHFC TF f M4l ) |) 4DT1 ME 
Q01 SsQOIS-*-««CENTHF<TFCM4I43-HL)«DTXME 
QSOLsQSOL ♦QC 40T IME 
52 IF(CT«GE.T4J GO TO 2 
30 CALL C0LL<CT*QC*Tl»T2«T9»T4*0ai 


oRiQii^Ai. f mge » 
0|t pedR QUALTCt 
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lFCCT«Ge*m 60 TO 6 

C IMA066UATB COLLECTOR OUTPUT AT ENO OP OAV» EVENIN6 

MCOLL»OC/OHMAX 
M>«MlN-liCaLL 
MT«II/XN 
6«»T/AFLT 
KJ «2 

1PCCT«£Q*0.I 60 TO 2 
0TlME«2/llT 

IPCOT1ME»LE«OTNAX| 60 TO 11 
OTlM£«OTMAX 

11 CTaCT40TlMF. 

IP<CT«LE*T1I 60 TO 12 
DT IMEoOT INE-CCT»T1 } 

CT*T1 

12 ETsET«OT 1NE/3E00* 

TP( 1 l«TL 

60 TO 10 

6 1PCCT«6E.T2I 60 TO 7 

C NO NEAT1N6 AVAILABLE PROM COLLECTORS 
ttaHMIN 
■T*W/XN 
6««T/AFLT 
K3*2 

CTlMEsZ/HT 

1P(0T1NE«LE«0TMAX1 GO TO 13 
CT1ME«0TMAX 

13 CT«CT«OTlME 
1F(CT.LE«T2) 60 TO 14 
CT1ME»0T1ME«CCT>T2I 
CT«T2 

14 ET«ET4OTlME/3600. 

TPI 1I«TL 

60 TO 10 

7 1PICT.GE.T3I 60 TO S 

C INADEQUATE COLLECTOR OUTPUT AT BEGINNING OF OAV* MORNING 

WC0LL«0C/0I1MAX 
W«MIM1K-MC0LL 
MT««/XN 
GbMT/APLT 

K 3 s 2 

0T1ME»Z/ET 

1FC0TJME.LE.0TNAXI GO TO IS 
OTIME«OTMAX 
16 CT«CT4CT1ME 

1FCCT«LE«T3| 60 TO 16 
OT1MESOT1ME-CCT-T3) 

CT»T3 

16 ETSET4CT1ME/360C. 

TF< 1I*TL 
GO TO 10 

C EXCESS SOLAR NEATING AVAILABLE FOR STORAGE 

6 OSTsQC'OB 

W*OST/C NU-'ENTNP ITPC 1 1 > I 
XFfCT.E0«T3) GO TO 2 
WC0LL«W46N1N 

IFC VCOLL.GT*«MAX> GO TO 20 
HlCsCV4ENTHPCTFf 1II4MM1N6HLI/K0LL 


c-?n 



c 


c 


c 


di 

dz 

20 


26 

27 

2 


60 

3 


T1C«FI.TEM*|H10 

MTsli/XN 

6«HT/APi.T 

K3«l 

0T1ME«2/6I’ 

iP(CTlM£«l.£«OTMAX| 60 TO 21 

0T1ME*DTNAX 

CT«CT60T IMS 

IPCCT«I.2»T4) 60 TO 22 

OT 1 ME«0T IMS- C CT> T4 1 

CT»T4 

ET«eT40TlME/3eOO* 

TPCM41|«TU 
60 TO 10 

INSUPPlCieNT PUMP FLOW CAPABILITY 
WSWMAX-WMIN 
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WTsW/XN 

6«WT/AFLT 

K3«l 

OTIMcsZ/WT 

iPCOTlME.LE*OTMAXI 60 TO 26 

0T1M£«DTMAX 

CTBCT40T1NE 

1PCCT«LE«T4I 60 TO 27 

OT lME»0TlME-iCT*T4l 

CT»T4 

£T«£T40T IME/3C00« 

TPCM41I«TU 
60 TO 10 
OUTPUT DATA 
liRlT6(6*220l 
WRlTEf 6.221) 

WR1TEC6*222) 

WR1T£(6*223I 
WRlT£C6t224) 

WHlTEC6*225i 
ttRlT£C6*226J 
WR1T£I6*227) 

•RITEC 6*2421 
«R1TEC6*2411 
WR1T£C6*2391 
«RITEC6*2401 
«RlT£C6*22dl 
WftlT£(6*229) 

00 60 l«l*M 
■RlT£f6*2301 
K7«M41 

WRITE! 6 *2311 K7*TFCK71 
WRlT£t6*2371 XN 
OTIME«OTMAX 
cT« ET «OT MAX/3600 • 
lPiCT*LT*T41 60 TO 61 
0E6IN NEW DAILY CYCLE 


ET 

CT 

OLOSS 

OSOL 

«>1S 

QCH6 

oeuRN 

OBOIL 

WPEAK 

POMAX 

PMAX 

WPUMP 


X*TECll*TPf 1)*VS€11*VLC11*AK6(11*ANLC1) 


CT«0* 

OLOSS>0 • 

0S0L«0B4CT1ME 

00IS«0. 

0CM6»0 • 
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aSURM^O* 

P0MAX«0« 

PMAX«0* 

•PUMP«0* 

WPEAK«0» 

NOiNC«i 

61 CT«CT40TMAX 

iPCNC*EO*NO> 60 TO 99 
60 TO 30 

99 STOP 
END 
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c 

c 

c 


3UaROUTlNli hTCACCCWTa6*OTlM£*KJ*M»OI AlcT*«M»eO»XLTE»VliL« 
THIS SUeROOTlNfe PEAIFORMS THfi HSAT ThAWSFER COMPUTATIONS ON THE 
ELEMENTS AND THE REOlSTAiaUTION OF LIQUID* 


C 


C 


SOO FURMAT(2FIE*8I 

801 FURMAT<//*8X»* INITIAL NITRATE CONCENTRATION •»F10*5I 

80«t FQRMATC8X** INITIAC STORAGE TEMPERATURE »K •*FI0*3> 

DIMENSION RSCl011«BEU)li01l»BLC10Ift*AS(101l *ALllOn •VVf lOi I »VEC 101 
l}*HC10il*CU01) *NAX( lOUvRAXSClOll *0011011 
COMMON TECl011*TFf 10n*VSC101l*VLtl011 •ANSUOlltANLUOl) 

COMMON/ OAT A/RHOL*RHOS* TTR *TEU*H«0 •HEUT*HA8T *CPEU*CPB *Clli *TKL 
COMMON/ VALS/H8TR • HLTR*HMT R *HATR*HAEU *HLEU *H«CU *CTR* CEU *C At U • TKATR • 
1TK8TR*TKAEU*TKEEU*CATR 
COMMON/ VC0N/R«AX*R1NS*R1NEN*TAM8*MIT 
IF(38*E0*99I GO TO 9 
INITIAL AND CONSTANT VALUES 
CACU«CPA(T£U) 

CATR«CPACTTR| 

CTR«CLIQCTTRI 

CEU»CLIO(TEU) 

HeTR«HB(TTR) 

HLTRBHL1Q4TTRI 
HMTRsHMO-fCHlATTR 
HATR«HA|TTR) 

HAEU«hACTEU> 

HLEU«HLI0CT£U1 
HKEUsNaO^CSATEU 
TKATRbTHCACTTRI 
TKBTRaTHCBCTTRI 
TXAEU>ThCA(TEU) 

TKEEU«ThCE4TEUJ 
XMaM 

AOVLsVEL«fXM/hITIA*E 
K1«M 
VVL>0* 

REACCSvSOOl CO*TO 
J7»X 

INITIAL CONDITION OF 
DO 8 1«1*M 
TEC I>«TO 
ASiI)*0* 

ALC II«VEL*RHOL 
ANAlflN 
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STUM AGE UNIT 


U*£0*1> AM««li«MB 
' I*£0*M| A9««M49B 

. I 1*0 • 

^XSCI)«0* 

VlC I l«VEL 
VSf ll«0* 

EELOl I|«Ali«(H804C»«Tai 
TL«TLIQICOJ 

ELI II«ALCll*(hLI01TL|4CPTCCUl»tTO-TLl> 
ANSCIIaO* 

8 ANLCIIaCOAALI II 
MRITEUaSOll CO 
MRlr£(6*902) TO 





c 


J9»99 

COMPUTE AXIAC CONDUCTION fffiSl STANCES AND MEAT PEONS 
9 1FCK1«£0«II CO TO 32 
K0«K1<-1 
OC 30 tslvKS 

1F(RAXSC1I«E0«0«I GO TO 34 
MTK«VCC/CA0VE*CVEU|4>TKI.4JILTE442/ftAXSCl I II 
60 TO 30 

34 RTK«I«/CA0VC4TIU.I 

30 RAXCII«l*y(l*/NTK41«/RNAX| 

32 DO 31 IbK 1*M 
IPCAAXSCII.EO.O.I 60 TO 35 
RAXC 1 1« 1 •/( 1 •/ROAX-»i •/MAXSi 1 1 1 
GO TO 31 

35 RAXCIIsftOAX 

31 CONTINUE 
¥ELE«VEL-VSC K1 1 
RAS0V«RAXCK1 14 WL/VEU. 

1PCRAXS€K1I«EQ«0*I GO TO 36 

RBEEOsVEL/C A0VL4 1 VEI.i4TKI.4XI. TE*42yRAXSC K1 1 1 1 
GO TO 37 

36 R0CEOsa«yCAOVL*TKI.I 

37 ReELOsCVELE^VVEI/CVELLOCl •/R«AX41«/RBiU.OII 
RAX CM 1 l«RABOV4ROEEO 
OO 33 t»2»N 

33 OCI ll«2«4€TCtl|oT£l3*lll/CRAXCll4RAXCl*i|| 

QCt llsO« 

OCCN41|sO. 

GO TOI1»101I*K3 

C HEATING SECTION 

1 OO 10 1«1*M 
J«M«141 
AHsNM 

1PCJ«£Q«1I AVs«h«.m8 
1FCJ«£Q*MI AVSHV4WB 

RAMesRlNS 

1FCJ.EQ..1I RAMB«R1I«N 
IFf J»EO«Ml RAMBsRlNEN 
C FIEN RESISTANCE* INSIDE TUBE 

TK«THCFCTFC34U I 
CP«CPFITFC341II 
VF*VlSFtTFC34ill 
REV«G401/VF 

IF I RET •ET •2000* I GO TO 2 

MlaO«0224fREV«60«8lVf(VF»CF/TK|440*3l«TR/01 
GO TO 11 

2 Hla3«664TK/0t 
11 RF«l*/fMl«AlETI 

IFCJ-Kll 21*22*23 

C FUEL EEEMENT* ASSUMED ISOTHERMAE 
21 CAEE RECALCCREI 
R0V«RV4RS f J 1 4RF4RE 
EFF» 1 «-EXPC-l •/ IROV4MT4CF I I 
TF( JI»TFCJ41l4EFF«fTECJI«TFC34t3l 
00F««T4C ENTHFITF C J I l-ENTHF ( TFC3413 I I 40T1 ME 
DQANBbC TEI3 l—TAMB l/RAMB 

ETOT«EE< J l4EEL0i J l-DQF«CDQANB40C( 3l*0Ct 3411 1 6DT 1 ME 
AT0T>ASC3|4AEC3I 
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€ 


C 


C 

c 


CALI.. FULLeL|ATOT«EYOT»OOaAII«XLVE»ALC«ll»ASCJI*IISC JI»TEC Jl •ELCJ>*K3» 
AEEL0C4I • ANL< A !• ANSC4 I •AAXSC 4 I } 

CO TO 29 

CCmPLCTELy SOLID* VARIASLE REFERENCE TEMPERATURE USED 
Z3 RaV«Rtt-»RSCJl'»RF 

EFP* 1 •-EJCPC- 1 •/ IROV*«T«CF I ) 

ALPH-AT«CF«CFF^I •/RARE 

SET AhQCC 4«I l>-QC C 4 14>TAMB/RANS« UTA^AEFF* TF C U 

GANMAbSETAFALPM 

1FICAMMA.CC.TEC4II CO TO «S 

CALL VRTCS«TEC4«*ECLO<4l*AS(41*AM*«T*CF*EFF*TPC4«>l»*XLTE»OTlAe*ANS 
1 C 4 1 •DO*TF C 4 I *RSf 4 I •RAXSC 41 *0014 I • OCC 4A1 I •RAMS •TAMBI 

GO TO 29 

4S CALL VRTM$CANSC4l*ASC4)*TCC4l*ECLO€4l»Afe*TFC4*>ll*CFF*CF»aT«OC*XLTE 
1»DT1ME»ANLC4|*ALC41*TFC4I*ELC4I*RSC4»*RAJISC4> •QCC4> •QCC441 1 •RAMa»T 
2AM8I 
GO TO 29 

THE ELEMENT INCLUOINC THE SURFACE* PARTIALLY FULL 
22 CALL RLCALCCRLI 
ROV»RII4ASC4 1 ARF4RL 
EFFsl •-EXPC tROVAHTACFI I 

CHECK FILLING RATIO. IF *GT. 0.2S* USE REOUCEO TIME INTERVALS 
IF *LT. 0*28* NEGLECT LIQUID AND USE V.R.T* 

VRAT 10* i .-VVL/VEL 

IFC VRAT10.lt .O.EI GO TO 70 

OTMOO*OTlNE/2. 


NX* 2 
GO TO 50 

70 1FCVRATX0.LT.0.343X GO TO 7i 
OTMOO*OTIME/3. 

Ml*3 
GO TO 50 

71 XF<VRATIO.LT.0.25» 60 TO 72 
OTMOOsOTlME/4. 

NX*4 
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so OQFT*0. 

OO 55 1SUB*1*N1 

TFO*TF(441l4EFF«CTEt4l-TFC44XI) 

0QF*MT4CENTHFCTF01-ENT»^CTF144UXXA0TM00 

OQFTsOQFTAOQF 

OQANe*f TEC4l-TAMei/RAMB 

ETOTSELC4 l♦CELO( 41>0QF-C OQAMB4QCC4 I~QC1 441 J ) ADTMOD 
AT0T*ASC4I4ALC4I 

CALL FULLELCATOT.ETOT.OO.AH.XLTE.ALC JI*ASCJI*RSC4I .TEC J1 .LLC JI.K3. 
1CEL0CJ)*ANLCJI*aNSCJ».RAXSC J l> 

55 CONTINUE 

HF03CNTHFCTFC441IX40QFT/I MT«0T1M£> 

TFI 4 l*FLTEMPC HFO I 
GO TO 29 

72 ALPHSMT4CFACFF41./RAMB 

BET AsQCC 441 l-QC 1 4 |4T AMa/RAM84MT«CF 4EFF« TF ( 44 1 1 

6AMMA*B£TA/ALPH 

IFC GAMMA.GE.TEI4II GO TO 75 

CALL VRTCSCTEI4 1 *EELO<4> • AS( 41 .AM* OT.CF *£FF *TF C441 1 • XLTE .DTI Mt • ANS 
1C4I*00*TFC4I*RS<4I*RAXSC41*QCC4».OCC4411.RAMB.TAMBI 
GO TO 29 
75 ELR*0. 

ALR*0. 
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ANLRsQ* 

CALL VRTNS4ANS<JI*ASCJI*TEC4l*EEL0CJl*Ali«TFCJ«l>*£FF«CF»«T»00*XLTE 
A «0T IME. AMLH* ALAtTFC J l•£LR .RSC J)*RAXS< J> .OCC J> •QCC JFU »fiAMB»TAMai 
ANi. C J |3 ANL ( 4 1 ♦ AMUR 
ALC JI*ALCJI«ALR 
ELC4|sCI.CJI<»ELR 
29 VSC J|«ASC^I/RHOS 
VL4 J>«ALC JI/RROL 
VV4Jf«VCL>VS4JI 
AO CONTINUE 
60 TO 200 

C COOLING SECTION 

101 Du AAO AsA«M 
AtIsliM 

AF4I«£a«AI AHsliW-t-lMI 
AF4I*E0«MI AMsHIi««i 8 
RAMBsRlNS 

1F4J.E0.A) RAMBsRANEN 
AF4J«E0*MJ RAHesRlNEN 
C FILM RESISTANCE* INSIDE TUBE 

TK»THCF4TF41>J 
CF«CPF4TF4AAI 
VFsVlSF4TF41»l 
REVsGAOI/VF 

IF4REV*LT*200€» GC TO 102 

Hl30*022*4REV«<l>C.eAF4OH<4>CF/TK|««0.Sl«TIC/0I 
GO TO All 

102 NlsS*6cATK/01 
AAA RFsl./CNJ«>A2£TA 

AF4I->KAI A2A*A22*A2S 
C FULL ELEMENT* ASSUMED ISOTHERMAL 

A2A CALL RLCALC4RL) 

HOVsRO-t-RS4 I l■»RF♦RL 

EFFsi A *>EXP4-1 */4ROV«HT*CF 1 1 

TF4 I«A|sTF4Il4EFF«(T£ClI-TF41JI 

OQF*MT*4£NTHFCTF4l4AI>-ENTHF4TF41l>l*DTlME 

DOAMBsC TEA 1)>TAMBI/RAMB 

ETOTsELC 1 14EEL04 1 »»OQF«4OQAMS40C4 1J-QC4 A4AA >*DT1ME 
AT0T«AS4IA4AL41I 

CALL FULL£L4ATOT*CTOT*00*AH*XLTE*AL41>*AS41 A*RS41A*TE41I •£L41A*K3* 
ACCL04 1I*ANL41»*ANS4II*RAXS41I| 

GO TO A29 

C COMPLETELY SOLID* VARIABLE REFERENCE TEMPERATURE USED 

A2S ROV3Rh«^RSCI}4RF 

EFF«l*«EXPC-A*/4ROV«in*CFIA 
ALPNSUTACFFEFF4A */RAMB 

BETA«0C4 14 AI>0C4 Ii4TAMa/RAMB-»ilTACFA£FF*TF4AA 

GAMMAsBETA/ALPH 

lF4GAMMA*LE*T£4ill GO TO 130 

CALL VRTHS4ANS41l*AS42)*Tc4ll*ECL04II*Ali*TF4U *EFF*CF*MT »D0*XLTE*0 
1T1MC*ANL41I*AL4 11*TF4A4A)*EL4 JI*RS4A»*RAXS4I>*('C4II*0C4I4A>*RAMB*T 
2AMSI 
GO TO 129 

130 CALL VRTCS4TE( 1I*EEL0C1|*A5CII*AM*BT*CF*EFF*TFI1|*XLTE*0T1ME*ANSCI 
AI*DO*TFC 141 J*RS4 I |*RAXS4 I 1 *QC4I I *QC4 |4A | *RAMB*TANB> 

GO TO 129 

C THE ELEMENT INCLUDING THE SURFACE* PARTIALLY FULL 

122 CALL RLCALC4RLA 
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AOV«ilM«RSC 1 l♦RF«m. 
eFF«t«-exPC-l •/CR0V4>IIT«CF» I 

C CHECK riU.IN6 RATIO* IF •€!?• 0*25* USE REOUCEO TIME IHTERVAi^ 

C IP .CT. 0*2«* NCCUECT LlQUlO ANO OS£ V.R.T* 

VRAT IQs 1 ••WL/VEC 
1F4VRAT10*LT«0«SI 60 TO 170 
OTMOOaOTlME/E* 

Nl«2 

60 TO 160 

170 1FC¥RAT10«I.T*0«333» 60 TO 171 
OTHOOsOTlNE/3* 

S5*fo 150 PAGE n 

171 1FCVRAT10«I.T*0«2S1 GO TO 172 OF POOR QUAUIV 

OTNOOsOTlME/4* ^ 

Nl*«« 

ISO OOFTsO* 

OO 15S 1$UB«1*N1 

TFOsTFl 1 |♦EFF«CTEC ll-TFI 111 

OQF^HTFC ENT I^CTFOI<>£NTHFCTFC 111140 TNOD 

OOFTaDOFT-*OOF 

OOANBffCTEC 1 l•TAil6 1/RAMB 

ETOT«Efc.C 1I4EE1.DC 1 loOQFofOQAMS-l-QCC II•0CC1♦1I I40TMOO 
ATOT»ASCll«AI.Cll 

CALI. FULLELCAT0T*ET0T«00»A«*XLTE»ALC11»ASC11*RSC1I •TECH cELC 11*K3« 
lEELOC 1 I • ANLI 1 I • ANSI 1 I •RAXSC 1 I 1 
ISS CONTINUE 

HFOsENTI^CTF 111 l♦OQFT/C MT40T1ME1 
TFC 1411SFLTEMPCHFO) 

GO TO 12S 

172 ALPHs|iT4CF4EFF41»/RAHB 

RET AaQCf 1-»1 1 •fX. C 1 l♦TAMB/RAMB♦HT«CF «£FF« TF C 1 1 

GANNAsBETA/ALPH 

IFC GANNA •LE.TEC 111 GO TO 175 

ELRsO* 

ALR«0« 

ANLRaO* 

CALL VRTHS1ANSC1I«ASC1}*TEC1I*££LOUI*AM»TFUI»EFF*CF*IIT»00»KLTE*0 
lTlN£*ANLR*ALR»TFf J«1 I•ELR•ASC1I•RAXSC II •QCCl l•OCCl♦l I •RANe*TAMSI 
ANLC tIsANLC 1 I^ANLO 
ALC 1I»ALC1|4ALR 
ELC 1I«ELCII4ELR 
GO TO 129 

a7S call ¥RTCSCT£Cll*£ELO(ll»ASCil»AH»NT*OF*EPF»TFCII*XLTE«OTlME*ANSCI 
1I*00«TFC141I«AS<1I«RAXSC1I*QCC1I.QCC141I*RANB*TAI»I 
129 tfSC 1I«ASCII/RH0S 
VlC l|aAL(ll/RHOL 
VVC IlsVEL'-’VSf II 
110 CONTINUE 

C LET LIQUID REDISTRIBUTE DUE TO GRAVITY ASSOCIATED HITH CHANGING 
C VOID SPACE AS SOLID MELTS OR LIQUID SOLIDIFIES 

200 VTOTLsO* 

3«1 

OO 210 lsl*M 

IFCVLCll«£a«O.I 60 TO 210 
VEC JI«VLCII 
CCJIsANLUI/ALCII 
HI2I«ELC1I/ALC1I 
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VTOTi.«VTCTL*VeCJl 

210 CUMTINUE 

C HAVE A K-CON6 ARRAY OF LIQUID VOLUMES AliO ASSOCIATED C AND H 
K*J-1 

lFCK«£0«0i GO TO 231 
VTOTV«0« 

OO 211 l«l»M 

VTOTVsVTOTVFVVC X I 

XF( VTOTV^VTOTLI 211»213»213 

213 Kl«l 
VVL«VTOTV-VTOTL 
GO TO 214 

211 CONTINUE 

C HAVE A Kl-LONG ARRAY OF VOIDS THAT IS 3UST LARGER THAN LIQUID VOLUME 

C REOlSTRlfiUTE THE LIQUID TO FILL SUCCESSIVE VOIDS 

214 3«1 

OO 219 Isl.Kl 
ALCaisO. 

ELI 11«0« 

ANLf l|s0« 

210 IFCVVIXI»VE4JII 21E*216«217 

216 VEC J|«VEC3I»VVC1I 

ALt 1|3ALCXI6VVC1|VRH0L 
ELI llsCLIIl4VVCll«AMOL«HlJi 
ANL € 1 leANLf I |4VV C 1 1 «RHOL«C( 3 J 
Ll« 1 

GO TO 21S 

217 ALC llsAL(Xl4VEC3|4RHQL 

ELC X|sEL(Xl4VEC3|4RH0L4HC3l 
ANL C 1 l»ANL ( 1 1 4 VE C 3 l«IMtOL4C C 3 1 
VVCI13VV4X1-VEI3I 
Ll«2 

IF43.EO.Kl GO TO 219 

33341 

GO TO 218 

219 VL4 1|3AL41I/RH0L 

220 KSSK141 

1F4KS-MI 221.240.232 

221 DO 225 IsKS.M 
ALC U»0. 

EL4 1I«0. 

ANL4I|sO. 

225 VL4 1|30. 

232 CONTINUE 
RETURN 
231 Kiel 
1 

VVL3VEL-VS41I 
GO TO 221 
240 ALCKS>«0. 

EL4KSI«0. 

ANL4KSI«0. 

VL4KS)«0. 

VV4KS1«VEL>VS4KE) 

RETURN 

END 
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orvnrf> 


c 


c 


c 


SUBROUT INIs PUIXei.CAT0T»ET0T»00«Ali»XLie*AC«AS*RS«T«EI.*K3«EEU>*ANlL*A 

'this* SUBROUTINE CONPUTES THE FINAL CONDITION OF AN IN|T I ALLY .? yMr- 
IMCI-IJOINC TOTAL NASS* MASS OF NITRATE* AND ENERCV LEVEL OP 
BOTH THE LlOUtC AND SOLID PHASE* AS BELL AS THE FINAL TEMPERATURE* 
AND THE SOLID H.T* RESISTANCE 

COMMON/DATAyRHDL*RH0S*TTR*TSU*MBO*NEUT*HABT*CPEU*CP8*CM«TKL 

COIMM)N/VALS/HBTR*M.TR*HMTR*MATR*HAEO*HLEU*Hlll£U*CTR*CBU*CAEU*fKATR* 

I TK8TR* TKAEU* TMEEU *CATR 
CISC ANL «ANS 1/ ATOT 

ETESTliATOTAMLIOCTlLI^AMFCHMOACBBTILl 
IFlCTOT*LT*CTESTl> CO TO I 
I^IQUID 

T^ ILAC ETOT-ETESTI 1/ ICNBA««ATOT*CPTCCI > J 
ASsO* 

RSsO* 

RAXSsO* 

ANLsANLAANS 

ANSsO* 

e&0sAmacmm04Cmbt> original ^AGE IS 

reTS?'®®’® ^ 

CHECK SOLID TRANSlTXm AT ALL BETA 
1 CEsCTR 

ALEsC ANL-IANS >/CE 
ASEsATOT-ALE 

ETESTSsASE«H0TR«ALE*l«.TR^AIl«MMTR 
IFCETOT*LT*ETESTEt CO TO E 

RE6ULA FALSI BETBEEN LIQUXOUS ANO TRANSITION 

E1*ETESTI-ET0T 

EEsETESTE'ETOT 

COsCE>fiE«CCI-CEi/CEI-EEI 

TOsTLlOCCOI 

ALOsC ANL4ANS I/CO 

aso*»atot*alo 

ETESTOs ASOBHBC TO l«AL0«HL10( TOI-FAB* CHBO»C BBTOl 
EOsETESTO-ETOT 
IFCEOI 20*EI*EE 
El T*TO 
ALSALO 
ASsASO 
ANL«ANL*ANS 
ANSsO 

EELOsABACHBO’FCBBTI-FASVHBC T) 

ELsETOT'EELO 
CO TO 25 
20 DO EE Jsl*a 

COsCO«EOB( C 1-CO 1/ C El^cOl 
TOsTLlOiCOI 
ALO«( ANL«ANS I/CO 
ASOsATOT-^ALO 

ETC STO* ASOAI^C TO I «>ALO«HLl Q ( T O •♦A BA IH BO*C «ATO } 

E6 COsCTESTOoCTOT 
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ao TO ti 
22 C1«C2 
Eisea 
eO TO 20 

C SOLID ASSIST ANCE 

2d VS«AS/RHOS 

OSsSQATf OOd«2^I •27S2SVS/XLTE) 

RSs ALOGC OS/OO l/C 6 •2832«TMCe4 TI«XLTS I 

AAXS«XLTE*«2/C VSATHCB4TI I 

RETURN 

C CHECK SOLID TRANSITION AT ALL ALPHA 

2 ET£STSaET£ST2«AS2AHAST 
XFCET0T*LT*ETEST31 GO TO 3 
ABs CET0T»£T£ST3|yHA8T 
AASAS2-AS 

AS«AS2 

ALSAL2 

T«TTR 

ANL«ANL-»ANS 

ANS»0 

£a.OBAH*H«TR«ASAHATR<«>Aa«HABT 
EL»ETOT«EELO 
GO T0428*29I»K3 

C SOLID RESISTANCE* COOLING AT TRANSITION 

29 VAsAA/RHOS 

OAsSORTf 00*«2-»I *2732«VA/XLTB1 

RAs ALOGf OA/OO >/C 6 •2832ATKATRAM.TE J 

V8SA8/RH0S 

OS^SORT C GA<t>« 2« I •2732d VB/XLTEI 
RSs ALOGI OS/O A 1/1 6«2832ATK8TR AKLTE 1 ♦RA 
RAXSbXLTS««2/C VA8TKATR-FVB8TKBTR1 
RETURN 

C SOLID RESISTANCE* HEATING AT TRANSITION 

28 VB«AB/RHOS 

OB« SORT C DOA* 2 * 1 *2732« VB/XLTE 1 

RBs ALOGf OB/OO l/C 6 •2B32ATKBTRAXLTC 1 

VAsAA/RhOS 

OS«SQRT C 0EAA2A1 *2732AVA/XLTE 1 
RSsALOGC OS/OB l/C 6*2e32«TKATR AXLTEl ARB 
RAXSVXLTEAA2/C VAATKATRAVBATKB-mi 
RETURN 

C CHECK LIQUID EUTECTIC 

3 CAsCEU 

AL4« C ANLAANS 1/C4 
AS4sATOT~AL4 

ETEST4«AS44HAEU4AL44HLEU4AM«HIBEU 
IFC ETOT.lt *ETEST41 GO TO 4 
C REGC8.A FALSI BETNEEN TRANSITION AND EUTECTIC 

E3SETEST3-ET0T 
E4SETEST4-ET0T 
CO« C4-E4 4 C C2-C4 l/C E3-E4 1 
TO=>TLX0CC01 
ALOeC ANL4ANS l/CC 
ASO«ATOT-ALO 

ETE STOs AS04HA C TO 1 4AL0AHL 1 0 C TO 14 A84 C H 804 C 8«T0 1 
EC2sETESTO>fiTOT 
IFfEO) 30*31*32 
31 TsTO 


C-30 



c 

c 

c 

c 

c 


A4.VAL0 

A$*ASO 

ANL«AML4ANS 

ANS«0 

£Q.OaA«A CMIlO«C««T )♦ ASOHAC Tl 

kl.sETOT-eEl.0 < 

OU TO 35 ( 

30 00 35 3«l«3 

CObCO-EOACC2<-CO 1/CE3-EO) 

T0»TLlQfC01 
AI.Oa C AM.4AMS I/CO 
ASO«ATOT-ALO 

ETESTOs ASOAMAC TO lAACOAHLlOf TOJAAHA IN MOACMATOI 
3e EO«ETESTO-£TOT 
GO TO 31 
32 C2»C4 
E3 bE4 
60 TO 30 

SOLID RESISTANCE 
3S VS^AS/RhOS 

OSb sort C 004424 1 «27324V S/XLTE I 

RSb AL06C OS/00 1/ C 6 •2S324TNCA C T I4XL TE I 

RAXSBXLTE442/CVS4TNCACTI I 

RETURN 

CHECK SOLID EUTECTIC 
4 ETESTSBETEST40AL44HEUT 

IFCETOT«LT«£TESTSI GO TO 5 
TKOoRHASE EUTECTIC 
AL5 bCETOT>ETEST5I/HEUT 
T»TEU 

ASbATOT-ALS 

ALbALS 

ELBAL4HLEU 

EELObETOT-EL 

ANT0TBANL4ANS 

ANL«C44AL 


ANSbANTOT-ANL 
SOLID RESISTANCE 
VABAS4/RH0S 

0 As SORT C 0044241 •27324VA/XLTE I 

RAb AL061 OA/00 1/ C « •2S324TK AEU4XLTE I 

V£B C AL4*AL6 l/RHOS 

OSBSORTC OA44241 .27324VE/XLTE1 

RSb AL06 < OS/ O A 1/ C 6 •ESSEATK EEU4KL TE > 4R A 

R AX SBXLTE442/ C VA4TKAEU4VE4TKEEUI 

RETURN 

ALL SOLIO BEL04 TEU 

S TbT EU-1 ETEST6*fiTOT 1/ 1 CM4Ali4C AEU4 AS44 ALA 4CPEUI 
ALbO* 

ASBATCr 

EELObETOT 

ELbO* 

ANSBANL4ANS 

ANLbO. 

SOLID RESISTANCE 
VABAS4/RH0S 

dab SORT ( D044241 .aTSaAVA/XLTEI 
RAbALOGC 0A/D0l/<6«2e3S4THCAf TI4XLTEI 
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Vfe«AC4/RH0S 

OSaSORT f 0A««a41 «a73a«VE/X4.TE > 
RS»ALOG<OS/OA l/4««a832*THCE€ TI«XI.T£I«RA 
RAX3»M.T£«*3/ ( VA«THCA( Tl« V£«TI^C T > > 
RETURN 
£NO 
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c 

c 

c 


SUBMUTIMfi VRTCSCT£«EEL0*AS«A«*liT*CF«EFP«TFl»XLT£»0TlME»ANS*00«TP0 
i 0J9S0MMS04H «4»^#jeAN0#7AM0> 

THIS SUBROUTINS USES VARIABLE REFERENCE TEMPERATURE TO DETERMINE 
THE COOLINC PROCESS IN A PURE SOLID S^MENT 


C 


* ^00# Qtuuix 


CONMON/OATA/RHOLvRHOS* TTR •TEUvMBO • HEUTvMABT »CPEU »CPB • CM • TKL 
COMMON/ VALS/HBTR • HLTRvHMTR •HATRvNAEU •HLEU*HM£U tCTR*CEU»CAEU tTKATHt 
i TK8TR*TKAEU« TKEEUvCATR 
A«HT«CF«fiFF»l«/RAM8 
a>02-Ql «T AMB/AAMS^MTACFAEFFM TF| 
o*a/A 

1FCAS.£0«0«J GO TO 10 
IFCTE-TTRI 1»S«S 

SOLID A60VE TRANSITION AT START* ALL BETA 
3 SMCaAS<l>CPe«-AM«CM 
TAU«SNC/A 

T0TsC<»< TE-C ISEXPC -OT IME/TAUl 
IFCTDT-TTRI 41 • 42 #42 
10 SMC»AM4CM 
TAUsSMC/A 

TOT«C4l T£>C}4EXPC«OTlME/TAUl 
0«AM4*CM4(TE-T0T t 
RAXS«0# 

GO TO 43 

42 Q«AS4CHB(T£I>HBCT0T1|4AM4CH*(TE«T0TI 

43 EGLO^EELO-Q 

QFsHT*CF 4EFF4(0T1ME«CC>TF1|4TAU4CTE-CI4 U#*EXP<-0TIME/TAU> )i 
TEsTOT 

HFO«ENTHFCTF 1 I40F/CMT40T1 NEl 
TFOsFLT EMP C UFO I 
RETURN 

41 DTI MIST AU4ALOGUTE*C)/CTTR-Cli 

OFs MT4CF4EFF4C DT 1MI4C C-TF 1 14 TAU* C TE-C1« < 1 •-EXPI^OTl Mt /TAUI 1 > 
OTlM2sOTIME<->DTlMi 
TFOTsTF I4EFF4 ITTR-TFI ) 

Qb I MT4CF4 C TF0T«>7F 1 |40Xo024 C T Tfi- T AMB ) /RAMBI40 T1 M2 
OQL«AS*HABT 
IFIO-OQLI S2»S2#83 

52 OTOTBAS4C HBf TEI-HBTR 14AH4CH4 4 T£« TTR 140 
EELOsEELD^OTOT 
OFs MT4CF44 TFOT-TF 1 1 40T 1 N24QF 
TEsTTR 

HFOsENTHFCTFi l♦QF/f MT«OTlMEl 
TFOsFLTEMPIHFOl 
AAsQ/HABT 
ABsAS-AA 
VA*AA/RHOS 

DA* SORT 4 0044 £4| •27324VA/XL TE I 
RAs ALOGI OA/00 1/40 •28324TKATA4XLTE i 
VSsAB/RHOS 

OSs SORT 4 O A44 24| • 27324V8/XL TE ) 

RSs AL064 OS/0 A 1/40 •2B324TKBTR4XLTE > 4R A 
RAXSSXLTE442/ 4 VA4TKATR4V84TKBTRI 
RETURN 

53 DT 1M3«AS4HABT/4MT4CF4 4TF0T-TF1|40I»Q244TTR*TAMBI/RAMB) 

OT IM4«0T |M2»0TIM3 
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QF« HT*CF«4 TFOT^TF 1 ) *OT 1M34QF 
8MC« AS«CATR« AMC« 

TAU«SMC/A 

TOT«C«< TTR-C I AEXP(-OTlM4/TAUi 
QaAS«CH8C Tei*HA<TOTII4>AH«CII«C TS-TOTJ 
fin.o«e£Lo-o 

QF»«T«CFACFFACOTlM4«<C«TFli<»TAU«€TTR-C>FCI*»SXPI*OTlM4/TAUll|40F 

F ACT«TnCe< TC I/THCA4 TOT I 

RSaRSAF ACT 

RAXS«RAX8«FACT 

TfisTOT 

HF0»6NTHFCTF1I«0F/(«T«CTIMEI 

TFO«Fi.TSMFIIiFO} 

RETURN 

C SOLID AT TRANSITION AT START* ALPHA ♦ BETA 
2 TFOT«TFI^EFF«fTE-TFJI 

Q«( MTFCFFCTFOT-TFl )<»01-a2«(TTR-TAMB>/RANai«OTlNE 
ObRTFCF* CTFOT'TF I lAOTlNE 
OEA«ETOT->AM«HMTR>AS«MATR 
IFCQ-OEAI S4*S4*55 

54 TFO«TFOT 
EELOsEELO-O 
TE«TTR 
RETURN 

55 0T1MI«0EA/CHT4CF4CTF0T-TFI|401>-Q24C TTR->TAHS}/RAMBI 
QF»«T4CF4CTF0T«TF J »40T1N1 

OTIMEbOT INE-CTIMI 
SNCBAS4CATR4A44CH 
TAU*SMC/A 

T0TbC4CTTR~CI«EXPC>0T1N2/TAUI 

0«AS4C HATR-HAf TOTI l<t-A««CH4f TTR«T0TI40£A 

Ea.DBEELD*Q 

aFsHT4CF«£FF4COTIN24(C*TFl>4TAU»CTTR-C)4U*>EXPC-OTlH2/TAUI))4QF 

TEbTOT 

HFObENTHFC TF 1 »4QF/< HT«OTl MEI 

TFObFLTEMP f HFO J 

RETURN 

C SOLID BELOW TRANSiTXON AT START 

I IFCANSI 61*61*62 
C PURE ALPHA 

61 SMCbASFCPACTEI^AWWCW 
TAU<*8MC/A 

T0T»C44 TE-C }*fiXP4-OT IHfi/TAUl 
OBAS44 HA4 TEI-HAC TOT I l♦AW•Ctll4 4 TE- TO T> 

SfiLOa>EELO-Q 

QPBWT4CF4EFF440T1NE44C-TF 1 14TAU44 Tfi'CI* I l.*EXP4-OTIME/TAUn J 
TE«TDT 

HFO«ENTHF4TFX |40FyillT40riME> 

TFO«FLTENP4 HFOl 
RETURN 

C ALPHA PLUS EUTECTIC 

62 AE«ANS/CEU 

AA«AS*AE 

SMC^AAWCAEU^AEWCPEU^AW^CW 

TAU«8MC/A 

TOT«C44TE-CI«fiXP4^TXNE/TAU} 

0«AA44HA| TEI-HA4TOT) l♦4Afi«CPEU«'A««CW)44 TE-TOTI 
EEL0«CEL0>0 


af^WT4CF4EFF«40TiME«CC-TFl |-»TAU*C TE-Cl* 4 t*>EXPC*>OTlME/TAU) >» 
TE;bTOT 

HFOSENTHF4TF1 t40F/4WT«0TlME) 

TFObFLTEMPCHFOI 

RETURN 

END 
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c 

c 

c 


10 


11 


14 


OUaROUTlNE VRTNS4ANS»AS»TetEELO*AM«TFlaEEF»CF*MT«OOaXI.TC«DTlMe»AML 
1 • AL • TFO • R8 • fl AXS • ai a 02 • R AM e« TAN6 1 

THIS SUeOOUTINE USES VARIABLE REFERENCE TEMPERATURE TO OETERMINE 
THE HEAT1N6 PROCESS IN AN INITIALLY PURE SOLID ELEMENT 

COMMON/OATA/RHOL aRHOSa TTR a TEUaHHO aHEUTaMAST aCPEUaCPBaCN a TKL 
COMMON/VALS/HBTRaNLTRaHMTRaHATRaHAfiUaHLEUaNMfiUaCTRaCSUaCAEUaTKArRa 
1 TKSTRaTKAEUaTKEEUaCATR 
A«VT«CF*SPF4I a/RANB 
8«02»01 ♦TAMB/RAMB4MT4CF4EFF* TFI 
OB/A 

IFCASaEOaO.I GO TO 30 paft£ |8 

1FCAN8I laSal .^TiTV 

PURS HVMlOXlOEa NO NITRATE ^ POOR OUALn* 

* IFITE-TTRI AaSaS OF PUU*' 

i 8NC«AS4CPB4AM*C« 

TAU«SMC/A 

T0T«C4f TE-C IFEXPI -OTIME/TAUl 

0«A S4 IHBI TOT l-HB ITS ) !♦ AHA CWA 4 TOT* TE I 

EELObEELOAO 

0F«IITACFAEFF«4 4C-TFI IAOTINE4TAUA4 TE*C|« 4 la*EXPC*OTIME/TAUI 1 1 
TE«TOT 

NF0«ENTNF4TF I l♦QFylHTAOTlMEl 

TFO«FLTEMP4HFOI 

RETURN 

> TFOT«TF14EFFA4TE*TFl) 

0« 4 HTACFA4 TF I-TFOT 140S*0 1*4 TiS-TANB 1/RANB l»OT IMG 
0EA«AM4HliTR4AS4HBTR*EEL0 
1F4 0*0EA| lOalOall 
EELO«£EL040 

0F«ST*CFA4TF0T*TF1I«0TIME 
tea TTR 

HFO«ENTNF4TFI H0F/4ilT«0TlNE) 

TF0<*FLTEMP4NF01 
RETURN 

OT lNl«0EA/4 MT4CF«4TF1-TF0T|402*01* 4T£*TAMBI /RAM8I 
QFas WT ACF* 4 TFOT*TF I ) AOT IM 1 
0T1M2«0TINE*0T1M1 
SMC* ASACPB4AH ACM 
TAU«SMC/A 

T0T»Ca4 TTR-C 1 AEXP4*0T1M2/TAUI 
0« ASANB4 TOT 1 4AH A 4 HMOACHATOT |*EELO 
EELO-EELOAO 

QF«MT ACF AEFFA40T IMSA 4 C* TF 1 1 A TAUA 4 T£*C1A 4 1 a*EXP4 *0T1NE/TAU1 1 ) AOF 
T£«TOT 

HF0*ENTHF4 TP 1 lAOF/4 MTAOTlME > 

TFr. .4.TEMP4HF0I 
RETURN 

SNC*ASACPA4 TEIaAMACW 
TAU*SMC/A 

T0T*CA4 T£*C > AEXP 4 *DT IME/TAUI 
1F(T0T*TTR) lAalAalS 
0«ASAHA4T0T|4AIIA4HH0ACMAT0T|*EGL0 
EELO-EELOAQ 

0F«MTACFAEFFA4 4C*TFI IA0TIMEATAUA4 TE*C1A 4 la*fiXP4*0TlME/TAU) I > 

TE*T0T 
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HFU«eNTtir c TP 1 140P/ C MT«OTt HE I 
TF0«P1.T£MP(MP0) 

RtTURN 

15 OTlMl«TAU5ALO«fCTE*CI/(TTR-C>> 

QF»IIT5CF5EFF5f IC'TFn^OTlNi^TAUFCTE^OFCl^-EXPC-OTlMX/TAUlll 

OTlMfl«OTlME-OTIMl 

TFO T«TP 1 4EFF •« TTR- TF I > 

0>C HT5CF5C TF t-TFOT >402-Ql-CTTR-TAMe>/RAMB| «0T1M« 

0EA«AS5HAaT 
iFfO-OEAl i7»17»aa 

17 a«a4AS«CHATR-HACTE} )4AM»Cll4tCTTR-TE > 

EEL0«fi£L04Q 

QF«MT5CF5(TF0T*TFX»*0T1M24QF 

Tt*TTR 

HFO«£NTHFCTFX »4QF/CliT*0TXMEI 
TFO«FI.TEMP<HFO» 

RETURN 

15 0TXM3«0EA/OHT«CF«CTFI>TF0T|402-01-iTYR-TA)/RAMBl 
QF« 5T5CF* ( TFOT'TF 1 ) «07 IM J4QF 
0TXM4S0T 1M2-0T1N3 
SMC«AS«CPe4AM*C5 
TAU»SMC/A 

TOT«C4C TTR-CMEAPI-OTIMA/TAUI 
Qb AS5HBC TOTl 4A5*CH«04CM*T0T 1-SEI.O 
E£1.0«£5LD4Q 

QF« MT 5CF «EFF5 < OT IMA* C C-TF t ) 4 TAU* < T TR<-C J * C 1 • -EXP ( -OTl N4/T AU > Il4 OF 
TE«TOT 

HFOsENTKFC TF 1 14QF/I «T*0T1 Mtl 

T FOBF1.T EMP i HF O I 

RETURN 

SOLID MITM N ITRATE (EUTECTIC I • MUST BE BEL05 TEU 
1 CI»CEU 
AEbANS/CI 
AAbAS**AE 

SMC«AA*CAEU4A£*CPEU4AH*CH 

TAU«SMC/A 

TUTsC4(T£-CJ*£XP(-OTlME/TAU) 

IFCTOT-TEUi 20»20*2I 

£0 Q>AA*( HA( TOTI'HAI TEl >4( AE*CPEU4AH*C« J*( TOT-TE) 

EEL0BEELD4O 

QFb 5T«CF*EFF*((C-TF1>*0TXME4TAUA(TE>CI*(1«-EXP(>DTIME/TAU} Jl 
Tt»T0T 

HFa«ENTNF(TFl )4QF/(MT*0T1M£1 
TF0«FLTEMP(NF0> 

RETURN 

£1 OTlMlsTAU*ALOC((TE-C|/(TEU-C»l 

QFbMTACFAEFF*! (C-TFI l»0TiM14TAUA( TE'O* ( 1 •-EXPf-OTlMl/TAUl J I 

01 1M2«0T1ME-0T1M1 

TFOT«TF14EFF«4TEU-TFII 

Qb ( «T*CF« 4 TF J'TFOT >4 02-01 - 4 TEU-TANB I /RAMBI *0T IM2 
OEAaAE*NrUT 
1F40«0EA> 23»23»24 
23 ALbQ/HEUT 

ObAA*(HA£U-HA4TE )>4AE*CPEU*( TEU-TE|4AL*MEUT4AM*CH*4 TEU-TEI 

£L«AL*HLEU 

EEL0bEELC40-EL 

QFsMT*CF«(TF0T-TFU*0TIM240F 

TE^bTEU 



c 


c 


2S 


^4 
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HfQmeHJHHJe 1 14QF/4 MT«0T1M£ ) 

TFO«FLTEMP(HFO) 

A^AS-M. 

ANL«C1*AL 
ANS«ANS~ANL 

SO. 10 resistance 

¥A«AA/RNOS 

D A« SORT 1 00«« 24 1 • E732«V A/XL TE I 
RA« AL06C OA/00 1/ C 6 •263E4TKAEU4XLTE i 
¥64 C AE«AL l/RHOS 
06«SQRT«0A4424l.2732*VE/Xi.TE> 

RSa At. 061 DS/00 |y ( 6 •283E4TKEEU4M. TE I4RA 
1F( VA«£0.0««ANO.VE«EQ*0.} GO TO 25 
RAXSBXLTE442/C VA4TKAEU4V54TKEEUI 
RETURN 
RAXS«0« 

RETURN 

BcYONO ElOUiO EUTECTIC* CESS ACCURATE AND LESS LIKELY 

DTi M3«0EAy ( MT4CF4 < TF 1-TFOT ) 4a2<-01> C T£U> TA} /RAMS) 

aF«MT«CF«(TF0T*TFII«0TIM340F 

0T1M4«0T1M2-0TIM3 

SMC»AA4CA£U4A£4CPT(C1)4AH«C« 

TAU4SMC/A 

0«>SMC4C TEU-C )* 1 1 *-EXP( -O T1M4/TAUI ) 

0VQ4AA4C HAEU-HA(TE) I4AE4HEUT4C AM4CS4AE4CFEU) 4>< TEU-TO 
0F»WT4CF4EPF«C0T 1M4* f C>- TF 1 }4 TAU« < TEU<>CI 4 C 1 • >EXF C -OT 1 M4/T AU 1 1 > 4 QF 
HFO«ENTHF( TF 1 )4QF/f ■T*OTI ME > 

TF0«FLTEMH(HF0| 

£TOT«E£L04Q 


K3«l 

CALL FULLELl AS*ETOT *00 *AM *XLT£*AL* AS*MS • TE •EL*K3*EEL0* ANL* ANS* RAXS 
il 

RETURN 

30 SMC>AW«CH 
TAU»SMC/A 

T0T«C4CTE-C«4EXPC<'0T1M£/TAU» 

QBAW4CH41T0T-TE I 
EEL0«EEL04CI 

QF« r 4CF4EFF41 C C- TM I*0T1 ME4 TAU» € T£«C>4 1 1 • -EXP 1 -DT 1 ME/T AU J 1 1 
TE«TOT 

HF0«ENThF4TFl)4QF/(MT»0TlMEl 

TFO«FLTEMPCHFQI 

RETURN 

END 
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SUBROUTINE COU.CCT •QC*T1» T2» T3.T4»QB»Ul 
C THIS SUBROUTINE READS A TABLE OF COLLECTOR OUTPUT VS* CYCLE TIME* 

C AND COMPUTES THIS OUTPUT AT A GIVEN CYCLE TIME 

C 

100 FORMAT! I3*4FI2*0*FI0*4*12J 

101 F0RMAT<F10*4*FI2*0I 
DIMENSION Q420i*TC20) 

IF(J*NE*99I 60 TO 10 
R£AC(5*I00I K*TI*T2*T3*T4*QB*J 
DO 1 lsi,K 

1 READ<S*I0II OCll*Ttll 
RETURN 

10 DO i ls£*K 

1F<CT*GE*TCII> GO TO 2 

QOQC 1-1}4(Q( 1>-QCI-II)«CCT«TC1-1II/CTC1I-T41>III 

GO TO 3 

2 CONTINUE 

3 RETURN 
END 


SUBROUTINE RLCALCCRLl 

THIS SUBROUTINE COMPUTES A LIQUID PHASE HEAT TRANSFER RESISTANCE 
IN THERMKEEP* 

FOR NON IT IS A DUMMY ROUTINE* SETTING RLsO 

RLsO* 

RETURN 

END 


FUNCTION FFCR> 

C THIS FUNCTION COMPUTES THE CIRCULAR TUBE FRiCTlCN FACTOR FOR 

C ROUGHNESS RATIO OF 0*0002 

C 

1FCR*GT *2000*1 GO TO 1 

FF»64 */R 

RETURN 

1 IF(R*€T*IOOOOO*I GO TO £ 

FF=0 • 1923«R*«!-0 .1772} 

RETURN 

2 FFs0«O2£ 

return 

END 


-- V-] 





c 

c 

c 


FUNCTION FLTEMPCNl 

THIS FUNCTION CONFUTES FLUID TENPCRATURC IN K FHQN ENTHALPY 
IN KJ/KC USINS H«0 AT AOOK 


FLT ENP«- III* 1«SQNT C 0*261 2E«0 6^S30 • SAHl 

NET URN 

END 


FUNCTION THCF4TI 

THIS FUNCTION GIVES FLUID THERMAL CONDUCTIVITY IN KJ/M-SEC«K 
FROM TEMPERATURE IN K 


THCF«0 .ISTOE^S-O^SaOE-TAT 

RETURN 

END 


ORIGINAL PAGE 18 
OF POOR QUALNY 


FUNCTION CPFITJ 

C THIS FUNCTION COMPUTES FLUIC SPECIFIC HEAT IN KJ/Kfa-K FROM 

C TEMPERATURE IN K 

C 

CPFsO*4ie9«0*277E-02»T 

RETURN 

END 


C 

C 

c 

c 


FUNCTION VISFCTI 

THIS FUNCTION COMPUTES ABSOLUTE VISCOSITY IN Ati/M-oEC FROM 
TEMPERATURE IN K 

FOR T>66 IT IS A FAIR APPROIMATION FROM 367K TO 569K 


X»l079./T-5*34 

VISF»10*V4X 

RETURN 

END 
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FUMCTION eWTHFCTI 

C THIS FUNCTION COMPUTES FLUID ENTHALPY IN KJ/K6 FRCM TEMPcRATUAE 

C IN K US IMS HsO AT 400K 

C 

£NTHF«C I • I7S40. ieS8E>C2« TI«C T-400* } 

RETURN 

END 


FUNCTION CPACTI 

C THIS FUNCTION COMPUTES THE SPECIFIC HEAT OF ALPHA-NAOH IN 

C FROM TEMPERATURE IN K* IT IS ASSUMED CONSTANT. 

C 

CPA32.01S 

RETURN 

END 


FUNCTION HAITI 

C THIS FUNCTION COMPUTES THE ENTHALPY OF ALPHA-NACH IN KJ/KO FROM 

C TEMPERATURE IN K. IT SILL HE ASSUMED THAI HlALPHAlsO AT 500K AND 

C ALL enthalpies MILL BE RELATIVE TO THIS. CP IS ASSUMED CONSTANT 

C 

HA>2.0ia4|T«S00.1 

RETURN 

END 


FUNCTION HBITI 

C THIS FUNCTION COMPUTES THE ENTHALPY OF SETA>NAQH IN KJ/KG FROM 

C TEMPERATURE IN K. CP OF BETA IS ASSUMED CONSTANT. 

C 

HBb 292.0842. 154|T'>S66. I 

RETURN 

END 
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FUNCTION TNCA4TI 

THIS FUNCTION COMPUTES THE THENMAl. CONDUCTIVITY OF ALPHA^NAOH* 
SOLID* IN KJ/SEC>M>« FROM TEMPENATURE IN K. 

THCA«0*001471 

RETURN 

END 


FUNCTION THCeCTI 

THIS FUNCTION COMPUTES THE THERMAL CONOUCTIVITV OF BETA-NACH* 
S4M.10* IN KU/SEC-M-K FROM TEMPERATURE IN K. 

THCBsO.001471 

RETURN 

END 
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FUNCTION THCECTI 

C THIS FUNCTION COMPUTES THE THERMAL CONOUCTIVITV OF EUTECTIC NA-OH* 

C NAN03* SOLID* IN KJ/SEC~M~K FROM TEMPERATURc IN K. 

C 

THCE30*001471 

RETURN 

END 


FUNCTION CLI04TI 

C THIS FUNCTION FINOS THE CONCENTRATION OF NANOS ON THE LIOUIOUS 

C LINE VS* TEMPERATURE IN K* C>MASSCNAN03I /TOTAL MASS.THERMKEEP. 

C THIS IS A SIMPLE VERSION UNTIL PHASE OlASRAM IS FOUND* 

C 

1F(T*LT«£6C*I SC TO I 
CLlQ«£*27eS>0.0C3E46«T 

i CLIOsj«ei84-0*0CE£164T 
RETURN 
END 
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FUNCTION TI.1Q<CI 

THIS FUNCTION COMPUTES TEMPEMATUNE IN K ON Tl*;> LIOUIOUS ClNc VS* 
CUNCENTRAT ION OF NAN03* C«MASSCNANC3|yTO TAi. NASS*THERMK£EP* 

THIS IS A SIMPCE VERSION UNTIC PHASE OIASRAM iS FOUND* 

1FCC*CT *0*101 60 TO 1 

TCl0«99«*•;^i00•«C 

RETURN 

1 TLIOsSaE*OS>ltO*87»C 
RETURN 
END 


FUNCTION CPTCCl 

THIS FUNCTION COMPUTES THE OIQUIO 7HEHMKEEP SPECIFIC HEAT AS A 
FUNCTION OF C0NC£NTKAT10N(KJ/K6-K| 

IT IS ASSUMEC TO EE A FUNCTION OF CONCENTRATION ONEV 
IT IS NOW A OOMMV FUNCTION WITH A CONSTANT VALUE 

CPTs£*14 

RETURN 

END 


FUNCTION HUIOITI 

This function computes the enthalpy of saturated solution! on the 

LIOUIOUS LlNEl AS A FUNCTION OF SOLUTION TcMPERATURE* 

IT IS NOW A SIMPLIFIED FUNCTION Nc6LECTlN6 PRIMARILY HEAT OF 
SOLUTION AS WELL AS AOSOLUTE VARIATION FROM NOS TO OH* 

HlIO»£*14«T>7CO* 

RETURN 

END 
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APPENDIX D 


LIST OP SYMBOLS FOR COMPUTER ANALYSIS 




'steel 


'TK 


shell 


shroud 


d 


solid 


m 


f 

h. 


h 


shell 


h 


shroud 


h, 


TK 


‘ins 


‘liquid 


k 


solid 


K 


steel 


kiK 


total T-66 flow area (in*) 
tube inside surface area per element (m*) 
specific heat of heat transfer fluid (kJ/kg-K) 
specific heat of steel (kJ/kg-K) 

Thermkeep minimum specific heat (kJ/kg-K) 
inside diameter of tubes (m) 
outside diameter of tubes (m) 
inside diameter of vessel (m) 
inside diameter of shroud (m) 
solid layer diameter around a tube (r> 
energy content of an end element (kJ) 
energy content of an interior element (kJ) 

T-66 friction factor 

T-66 side heat transfer coefficient (kW/m*-k) 

inside height of vessel (m) 

inside height of shroud (m) 

enthalpy of Thermkeep (kJ/kg) 

T-66 thermal conductivity (kW/m-K) 

thermal conductivity of insulation fkW/m-K) 

thermal conductivity of Thermkeep, liquid 
(kW/m-K) 

thermal conductivity of Thermkeep, solid 
(kW/m-K) 

thermal conductivity of steel (kV//m-K) 

average axial thermal conductivity of 
Thermkeep in an element (kW/m-K) 


ORIGINAL PAGE IS 
OF POOR QUALITY 


D-1 



L 


T 


N„ 


N 


PR 


N 


RE 

NTU 

P 

P 

*^ax 


ex 


^ins ,e 


‘TOT 


Q 

Q 

Q 

Q 

Q 

Q 

R 

R 

R 


abs 

burn 

c 

del 

env 

loss 

aspect 

f 

Ins , t 


length of each tube (m) 
total mass of Thermkeep (kg) 
number of spatial elements 
T-66 Prandtl number 

OF POOR QUALITY 

T-66 Reynolds number 

total number of heat exchanger tubes 

number of transfer uni' * in an element 
ideal pumping powei’ (k.\> 

axial heat transfer rate between elements (kW) 
heat rate required for boiler (kW) 
collector output (kW) 

excess heating provided by collectors (kW) 

heat transfer rate to fluid per element (kV/) 

heat transfer rate to environment per element 
(kW) 

total heat transfer rate to an element (kW) 

total daily heat absorption by storage (kJ) 

total daily make-up heating required (kJ) 

total daily output of collectors (kJ) 

total daily heat delivered from storage (kJ) 

total daily heat lost to environment (kJ) 

total dally collector heat not used (kJ) 

aspect ratio, vessel height/vessel diameter 

resistance to heat transfer, fluid side (K/kW) 

lateral resistance to heat transfer of insula- 
tion per element (K/kW) 
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^shell 

^solid 

^T,a 

^TK,a 

^TK,a 

\oT,a 

^TOTAL 


^end 

^ins 

^shell 

^shroud, 

bottom 

^shroud, 

side 

^shroud, 

top 

T V. 

amb 

Te 



^o 


vessel axial heat transfer resistance per 
element (K/kW) 

resistance of solid to radial heat transfer 
(K/kW) 

tube axial heat transfer resistance per 
element (K/kW) 

Thermkeep axial heat transfer resistance per 
element (K/kW) 

average Thermkeep axial resistance between 
elements (K/kW) 

tube lateral heat transfer resistance per 
element (K/kW) 

total axial resistance between elements (K/kW) 

total elemental heat transfer resistance (T-66 
and Thermkeep) (K/kW) 

vessel end plate thickness (m) 
insulation thickness (m) 
vessel side wall thickness (m 
shroud bottom thickness (m) 

shroud side wall thickness (m) 

shroud top thickness (m) 

ambient temperature (K) 

storage outlet temperature-upper (K) 
collector outlet temperature (K) 

fluid temperature into an element (K) 

boiler outlet temperature of T-66 (K) 

storage outlet temperature-lower (K) 




T fluid temperature out of an element (K) 

C-/ LI 

TjTp, element temperature (K) 
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^solld 

V 

TK 

\0TAL 

V 

tuDes 

''coll 

''min 


W. 


W 


heat 

exchanger 

ins 


W, 


w 


shroud 


w, 


vessel 


z 

AE 

AM 

AM 

AM, 

AP 

At 

C 


St ,e 
st ,m 


TK 


pf 

^ins 

^steel 

^TK 


solid layer volume per tube per element (m®) 
total Thermkeep volume (m^) 
total internal volume of vessel (m®) 
total volume occupied by heat exchanger (mM 
T-66 mass flow rate from collectors (kg/sec) 
fluid (T-66) mass flow rate to TES unit (kg/sec) 
T-66 mass flow rate to boiler (kg/sec) 
heat exchanger weight, tubes only (kg) 

insulation weight (kg) 

total pump work (ideal) per dally cycle (kW-hr) 

shroud weight (kg) 

total vessel weight (kg) 

factor for time step calculation (kg) 
energy change in an element per time step (kJ) 
mass of steel per end element (kg) 

mass of steel per interior element (kg) 

mass of Thermkeep per element (kg) 

heat exchanger pressure drop (kg/m^) 
time step (sec) 

effectiveness of heat transfer in an element 
T-66 viscosity (kg/m-sec) 

fraction of Thermkeep in solid phase in an 
element 

T-66 density (kg/m^) 
density of insulation (kg/m^) 

density of steel (kg/m^) 

Thermkeep density, liquid phase (kg/m^) 
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APPENDIX E 


LIST OP SYMBOLS FOR SCALING PROCEDURE 

A heat exchange surface area 

Aj^ heat exchange surface area per unit length 

a cross sectional area of storage medium 

a^ tube cross-sectional area 

c fluid sptclflc heat 

c storage medium specific heat 

s 

d . diameter of shroud 

sh 

d^ diameter of vessel 

h^ height of vessel 

k thermal conductivity of storage medium 

k^ thermal conductivity of insulation 

k^ thermal conductivity of tube 

L^ tube length 

M Storage medium total mass 

storage medium mass per unit length 

m fluid mass flow rate 

P reference time period 

s distance along tube 

t storage medium temperature 

fluid temperature 
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reference temperature 


shroud surface temperature 

overall heat transfer coefficient 
fluid volume 

fluid volume per unit length 

distance along flow direction in storage unit 

dimensionless storage medium temperature, t/tj^ 

dimensionless time, t/P 

dimensionless distance, x/h^ 

fluid density 

time 
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